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Préface 


Moustiques. Un mot générique bien pratique pour simplifier une grande diversité 
biologique, écologique, sanitaire, sociale, économique et historique. Les moustiques, 
Culicidae de leur nom scientifique, regroupent les 3600 espèces décrites sur terre, 
plus toutes celles non encore décrites. Si la nature a horreur du vide, les moustiques 
aussi. On les trouve sur tous les continents, dans tous les écosystèmes, et ce, depuis 
bien plus longtemps que les humains. 


Mais, lecteur humain, tu nen as pas conscience. Seul te préoccupe le moustique qui 
empêche de dormir, de profiter de la soirée en plein air et, selon le lieu où tu vis, le 
moustique responsable de maladies dans ta famille ou chez tes animaux domestiques. 


Émergence. Ce mot est devenu à la mode ces dernières années. Les humains ont pris 
conscience que des maladies pouvaient émerger. Charles Nicolle dans Destin des mala- 
dies infectieuses écrivait déjà en 1933 : «il y aura des maladies [infectieuses] nouvelles. 
C'est un fait fatal». Certaines de ces maladies, qui émergent de foyers sauvages, à 
l'issue de changements environnementaux, climatiques, démographiques, sociétaux, 
culturels, sanitaires, économiques, etc. sont à transmission vectorielle, et parfois les 
moustiques en sont les responsables en inoculant virus et parasites. 


Les moustiques sont des insectes, mais leur étude et leur contrôle vont bien au-delà de 
Yentomologie (entoma, l'insecte en grec ancien). De nombreuses disciplines complé- 
mentaires sont mobilisées, de la taxonomie à la santé publique. La télédétection et la 
modélisation spatiale en font partie et sont maintenant des composantes indispensables 
de lentomologie médicale et vétérinaire, mais aussi agricole. 


Hippocrate, philosophe et médecin grec, fait déjà, au v® siècle avant notre ère, le 
lien entre environnement et maladies. Il décrit des fièvres ayant la symptomato- 
logie du paludisme et observe une association entre les marais et ces fièvres dans 
son traité Des airs, des eaux et des lieux. Évidemment, à l’époque, si on se plaignait 
probablement des moustiques, le lien de causalité avec le paludisme était impos- 
sible à faire. Plus proche de nous, sur l'actuel territoire hexagonal, la Vendée, la 
Sologne, les Dombes, la Camargue étaient envahies de moustiques et les fièvres 
y étaient fréquentes jusqu'au début du xx® siècle. La construction du château de 
Versailles, avant canalisation des eaux de surfaces, a été le théâtre de nombreux 
décès probablement dus au paludisme. 


C'est seulement à l'aire pastorienne, donc très récemment, que les liens de causalité 
entre milieux, climat, moustiques, microbes et maladies ont été précisés. Depuis 
20 ans nous sommes maintenant capables d'aller beaucoup plus loin dans la compré- 
hension de ces associations, grâce aux nouvelles techniques de génomique, mais aussi 
grâce à la télédétection, à l'analyse spatiale des phénomènes biologiques, moustiques 
compris, et à la modélisation des risques sanitaires. 


Télédétection et modélisation spatiale 


La diversité biologique des moustiques est extraordinaire. Les 3 600 espèces sont parti- 
culièrement bien adaptées à des milieux et biotopes spécifiques. Certaines espèces ne 
se trouvent, à l'état larvaire, que dans des gîtes très particuliers, comme des phyto- 
telmes arborés ou des urnes de plantes carnivores comme les népenthès. D’autres sont 
moins exigeantes et se développent dans des lacs, des marécages, des bords de rivières; 
d’autres enfin se trouvent presque uniquement dans des collections d’eau artificielles 
créées par les humains. Certaines espèces ne sont connues que dans une seule zone 
(Aedes pia à Mayotte) et d’autres, profitant de l'environnement urbain, se retrouvent sur 
tous les continents (Aedes albopictus). Certaines se nourrissent du sang de nombreux 
animaux, dont l'humain (Anopheles arabiensis), alors que d’autres ont des régimes 
alimentaires très stricts (miellat de fourmis pour Malaya sp.). Plusieurs espèces sont 
capables de passer la saison sèche ou froide sous forme d'œufs en diapause (Aedes) ou 
sous forme adulte au repos dans des abris, dont les maisons ou les étables. Mais tous 
ont besoin d’eau pour pondre leurs œufs et permettre le développement des larves 
et des nymphes. L'eau, par sa présence, sa qualité, ses caractères physico-chimiques 
et biotiques (plantes, nourriture, prédateurs), est donc un paramètre essentiel de la 
biologie des moustiques. Toute approche visant à décrire, analyser et corréler les para- 
mètres liés à leau (pluviométrie, aménagements, végétation, etc.) permet de mieux 
estimer, voire prédire, la présence et l'abondance des différentes espèces et populations 
de moustiques, et les risques afférents. 


Ces risques ne sont pas anodins. L'histoire est riche de destins contrariés par les mous- 
tiques, de la mort d'Alexandre le Grand attribuée au paludisme (Anopheles) ou à la 
fièvre West Nile (Culex), au creusement du canal de Panama freiné par le paludisme 
et la fièvre jaune (Aedes), jusqu’à l'abandon récent, par leurs habitants, d'immeubles 
végétalisés envahis de moustiques-tigres en Chine. La liste des agents infectieux 
transmis aux humains et aux animaux par les différentes espèces de moustiques est 
impressionnante. Près de 100 maladies humaines sont attribuables aux moustiques. 
Certaines sont encore rares, comme la fièvre de Mayaro en Amérique du Sud. D’autres 
sont trop fréquentes comme le paludisme, qui tue près de 400000 enfants par an en 
Afrique, ou la dengue, qui touche, d’après l'OMS, plus de 300 millions de personnes 
par an sur tous les continents. 


Les enjeux sanitaires, sociaux et économiques liés aux moustiques sont donc colos- 
saux, mais les enjeux écologiques le sont également. Si les moustiques ont bien 
leur place dans les chaînes trophiques et sont une composante de la biodiversité, il 
nen reste pas moins que le contrôle des quelques espèces responsables de grandes 
maladies humaines (et animales) est nécessaire. Ce contrôle doit être rationalisé, 
intégré, ajusté, durable, accepté et le moins impactant possible sur l'environnement. 
L'ère de l'usage intensif d’insecticides se termine. D’autres méthodes plus ciblées, y 
compris géographiquement, sont en développement. Les approches de télédétection, 
d'analyse spatiale et de modélisation sont devenues indispensables, bien qu’encore 
insuffisamment utilisées dans l’aide à la décision pour atteindre ces objectifs. 


Les exemples pris dans cet ouvrage — Anopheles et risques de paludisme en 
Camargue, en Guyane française, en Asie ou à Madagascar ; Aedes et risques de dengue 
en Thaïlande, au Brésil ou dans l’océan Indien — démontrent que la télédétection et la 
modélisation spatiale, appliquées aux moustiques et aux maladies dont ils peuvent étre 
responsables, sont indispensables. Ils démontrent également que l’interdisciplinarité 
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est nécessaire. Des modèles bâtis sur des données biologiques mal documentées, non 
seulement sont sans intérêt, mais peuvent donner de faux espoirs à ceux qui les utili- 
seront. Inversement, un travail d’échantillonnage de terrain rigoureux restera mal 
exploité sans une bonne modélisation spatiale. 


Chaque communauté scientifique a ses concepts et son langage. Se tromper de sympo- 
sium spécialisé peut devenir un calvaire si on nen maîtrise pas les codes. Lâchez un 
entomologiste dans un congrès sur la télédétection, ou un géomaticien au meeting 
annuel de la Société d'écologie des vecteurs, il n'est pas certain que les mots ou acro- 
nymes comme réflectance, fauchée, exophile, sternite, univoltin, mode raster, diapause, 
signature spectrale, gonotrophique, MODIS, NDVI et autre NDW/, soient perçus 
comme ils le méritent. Seul le mot vecteur sera peut-être compris par tous, mais avec 
deux sens bien différents : biologique et géomatique. Le mérite de ce livre écrit par des 
spécialistes ayant pu se confronter, voire appartenant, aux deux communautés est de 
rendre accessibles ces concepts à l’aide d'exemples concrets bien documentés. Merci et 
félicitations aux auteurs. Cet ouvrage sera une référence très utile pour tous ceux qui, 
conscients de la nécessité d’une approche globale, spatiale et environnementale, pour 
étudier les moustiques (et au-delà les autres vecteurs) et documenter leur biologie, 
leurs distributions, leurs impacts et leur contrôle, sont à la recherche d'exemples pour 
comprendre et utiliser la télédétection et la modélisation spatiale. 


Ce livre est un pont entre les communautés invitant les entomologistes à « prendre de 
la hauteur », et les télédétecteurs et géomaticiens à découvrir le monde fascinant des 
moustiques. 


Didier Fontenille 
Entomologiste médical, directeur de recherche 
IRD, UMR MIVEGEC (université de Montpellier, IRD, CNRS) 
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THIBAULT CATRY, ÉRIC DAUDÉ, NADINE DESSAY, ANNELISE TRAN 


La télédétection fournit des données d'observation de la Terre qui peuvent être parti- 
culièrement utiles pour la cartographie et la modélisation en santé. Pour l'Organisation 
mondiale de la santé (OMS), l'identification, le suivi et le contrôle des populations 
d’arthropodes vecteurs constituent la priorité en matière de surveillance des maladies 
à transmission vectorielle. À ce titre, depuis plus de deux décennies, de nombreux 
travaux ont montré que les images satellite, et plus généralement les données spatia- 
lisées, permettent l'identification des variables environnementales et climatiques qui 
influencent les dynamiques spatio-temporelles des maladies vectorielles, et notamment 
celles dont agent pathogène est transmis par des moustiques. La diversité actuelle des 
capteurs satellite permet d'accéder à des données à une résolution spatiale et tempo- 
relle suffisamment élevée pour (i) caractériser différentes variables environnementales 
et climatiques (occupation du sol, précipitations, température, humidité, etc.) asso- 
ciées à la présence d’habitats favorables, à l'occurrence et à l'abondance des vecteurs, 
(ii) développer des méthodes et outils de prévision des risques associés à ces vecteurs et 
aux agents pathogènes qu'ils transmettent, et ce, à différentes échelles et (iii) contribuer 
à la surveillance de l’évolution de ce risque. Ces travaux reposent sur des méthodes avan- 
cées de traitement d'images satellite (classifications pixel ou objet, utilisation de séries 
temporelles, algorithmes d'intelligence artificielle, voir partie 1) multicapteurs (optique, 
radar, lidar, etc.) et multirésolutions (moyenne, haute, très haute résolution spatiale), 
ainsi que sur l'intégration de ces variables issues de la télédétection couplées à d'autres 
types d’information spatiale dans des approches de modélisation (basées sur les connais- 
sances, les données, les processus ou les comportements, voir partie 2). Ces modèles 
permettent la prise en compte d’un grand nombre de variables (notamment environ- 
nementales et climatiques) dans des systèmes complexes dynamiques et améliorent la 
compréhension de l'épidémiologie des maladies vectorielles dues aux moustiques et des 
mécanismes de leur transmission, qui représentent un enjeu majeur de santé publique. 


Issues de ces travaux développés dans le cadre de programmes de recherche, des 
méthodes opérationnelles reposant sur la télédétection et la modélisation se sont 
multipliées ces dernières années dans le domaine de la santé. Des produits (cartes de 
risque, chaînes de traitement) ont ainsi pu être mis à disposition par des infrastructures 
comme le pôle de données et de services surfaces continentales Theia! et en particulier 


1. https://www.theia-land.fr/ 
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par son centre d'expertise scientifique (CES) «Risques maladies infectieuses »?. De 
telles initiatives ont permis de fédérer des communautés de chercheurs, venus de la 
géomatique et d’autres disciplines comme l’entomologie et l’'épidémiologie, ainsi que 
des acteurs de la santé publique autour de l'objectif commun d'améliorer les connais- 
sances et les outils de lutte contre les maladies vectorielles dues aux moustiques. En 
particulier, le projet Anisette? (Analyse inter-site : évaluation de la télédétection comme 
outil prédictif pour la surveillance et le contrôle de maladies causées par des mous- 
tiques), financé par le Centre national d'études spatiales (Cnes) entre 2018 et 2022, a 
eu pour objectif de mesurer l’interopérabilité de méthodes couplant télédétection et 
modélisation spatiale pour prédire la dynamique des moustiques vecteurs et des mala- 
dies associées. Il s'appuie sur les résultats de différents projets de recherche, clôturés 
ou en cours, menés par des équipes issues de différentes unités mixtes de recherche 
(ASTRE, Espace-Dev, IDEES et TETIS) qui abordent la modélisation entomologique 
en étroite collaboration avec d’autres organismes comme le réseau des instituts Pasteur, 
et en particulier l’Institut Pasteur de Madagascar. L’interopérabilité des méthodes 
développées dans le cadre de ces projets a été testée sur plusieurs sites en Amérique 
du Sud (Brésil, Antilles, Guyane), en Europe (France), dans l'océan Indien (Madagascar, 
Maurice, La Réunion), en Asie du Sud et du Sud-Est (Inde, Thaïlande, Cambodge) et en 
Océanie (Nouvelle-Calédonie). 


Le présent ouvrage est le fruit des travaux et réflexions menés dans le cadre du projet 
Anisette, et présente de manière synthétique les concepts théoriques, les approches 
méthodologiques, les outils et les principaux résultats obtenus par l’équipe du projet, 
composée principalement de géographes, géomaticiens et modélisateurs. Avec pour 
objectif d'introduire les concepts de télédétection et de modélisation spatiale appli- 
quées à l'étude des maladies impliquant des moustiques vecteurs, il s'adresse à un 
public non expert souhaitant s'initier à ces notions et à leurs applications en santé 
publique*. L'ouvrage s'organise en deux parties distinctes : la première partie aborde 
les méthodes de télédétection pour l'identification et la caractérisation des détermi- 
nants environnementaux et climatiques des populations de moustiques vecteurs. La 
seconde partie se focalise sur l'intégration de ces variables dans différentes approches de 
modélisation pour la mise en place d'outils opérationnels de surveillance des maladies 
vectorielles dont sont responsables certaines espèces de moustiques. Afin de faciliter la 
lecture, plusieurs termes techniques — liés aux différents domaines concernés : ento- 
mologie, épidémiologie, télédétection, géomatique, mathématiques — sont définis au 
fil du texte et regroupés dans un glossaire final. Une liste des sigles utilisés et de leurs 
déclinaisons est aussi proposée en fin d'ouvrage. 


» Notions de télédétection 


La télédétection est définie comme l’ensemble des techniques utilisées pour collecter 
des informations sur des objets à distance; en particulier, l'observation de la Terre 
utilise un instrument (le capteur), embarqué à bord d'une plateforme (satellite, avion, 


2. https://www.theia-land.fr/ceslist/ces-risques-maladies-infectieuses/ 

3. https://anisette.cirad.fr/ 

4. Pour des informations plus complètes sur la biologie et l'écologie des moustiques vecteurs, le lecteur 
pourra se référer à l'ouvrage Le moustique, ennemi public n° 1 ?, coordonné par S. Lecollinet, D. Fontenille, 
N. Pagès et A.-B. Failloux, et récemment publié aux Éditions Que. 
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drone, etc.), pour caractériser la surface de la Terre (les surfaces continentales, les 
océans ou l'atmosphère). Les exemples typiques de télédétection incluent l’utilisation 
de l'imagerie satellitaire ou de la photographie aérienne. 


Principales caractéristiques des capteurs de télédétection 


Il existe différents types de capteurs en télédétection. Les capteurs passifs mesurent le 
rayonnement naturel émis ou réfléchi par la surface observée, ce sont par exemple les 
capteurs optiques, qui dépendent d'une source d'énergie externe ` la lumière du soleil. 
Les capteurs actifs mesurent, eux, la réflexion du rayonnement qu'ils émettent, c'est 
le cas des radars, qui émettent leur propre source d'énergie, et mesurent la rugosité et 
l'humidité des surfaces. 


Les signaux mesurés par les capteurs de télédétection sont appelés «rayonnement 
électromagnétique » et possèdent des propriétés que l’on peut mesurer et qui servent 
à les décrire : la longueur d'onde (représentant la périodicité spatiale des oscillations, 
c'est-à-dire la distance entre deux maximas de oscillation), l'amplitude (ou intensité, 
correspondant à la valeur maximale de l’oscillation) et la polarisation (relation entre 
l'amplitude et la direction de propagation de l'onde). Les capteurs mesurent la quantité 
d'énergie transportée par le rayonnement électromagnétique émis ou réfléchi par les 
surfaces et, en particulier, l’albédo ou réflectance hémisphérique directionnelle, défini 
comme le rapport entre l'énergie réfléchie et l'énergie reçue. Les images panchroma- 
tiques, en noir et blanc, sont obtenues à partir de l'enregistrement de ce rayonnement 
dans un unique intervalle de longueur d'onde. Lorsque cette mesure s'effectue dans 
différentes longueurs d'onde, on parle de mesure « multispectrale». Les capteurs de 
télédétection peuvent mesurer des signaux dans le spectre du visible (télédétection 
optique), de l'infrarouge ou des micro-ondes (télédétection radar), apportant ainsi des 
informations complémentaires. Les capteurs peuvent être situés au sol, sur des avions 
ou des drones (on parle de capteurs aéroportés) ou sur des satellites. La plupart des 
capteurs d'observation de la Terre organisent l'information acquise sous forme d'images 
numériques, caractérisées par la taille de pixel et la largeur de la prise de vue (la fauchée). 


Les capteurs de télédétection sont principalement caractérisés par trois résolutions : 
— la résolution spatiale est la taille minimale des objets observables, étroitement liée à 
la taille élémentaire de chaque pixel de l’image. On parle de basse résolution pour des 
images dont la résolution spatiale est supérieure à 1km, de moyenne résolution quand 
elle est comprise entre 100 m et 1 km, de haute résolution entre 10 et 100m et de très 
haute résolution pour des valeurs métriques ; 

— la résolution temporelle qui correspond au temps de revisite d’un satellite 4 un 
endroit donné, c’est-à-dire l'intervalle de temps entre deux prises de vue identiques; 
— la résolution spectrale qui caractérise la capacité d’un capteur a distinguer des 
signaux de longueurs d'onde différentes. 


La notion de signature spectrale 


Chaque type de surface est caractérisé par sa signature spectrale, définie comme la 
variation de la réflectance en fonction de la longueur d'onde (figure 1.1). La signa- 
ture spectrale dépend de la nature de la surface, de ses propriétés physiques et de 
l'interaction qu'elle produit avec l'onde électromagnétique qu'elle reçoit. 
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Figure I.1. Les signatures spectrales des surfaces naturelles dans différentes gammes de 
longueurs d'onde. Adapté de https://e-cours.univ-paris1.fr/. 


Extraction de l'information des images satellite 


Il existe différentes façons d'exploiter les informations spectrales contenues dans les 
images satellites. La méthode la plus simple consiste en une interprétation visuelle 
de l’image ou photo-interprétation. Des méthodes d'analyse plus complexes sont 
utilisées pour réaliser des classifications du contenu spectral des images, en se basant 
sur l'information contenue dans les pixels (approches orientées pixels) ou dans des 
objets qui sont des regroupements de pixels en régions homogènes (approches 
orientées objet). Ces approches de classification peuvent être non supervisées (sans 
connaissance a priori sur l’image à classer) ou supervisées (quand des connaissances 
a priori sont disponibles) [figure 1.2]. Les algorithmes de classifications (K-means, 


Image optique Sentinel-2 
Classification non supervisée 


Algorithme de clustering = Labellisation 
(p. ex. K-means) des classes 


Étapes de production d’une classification = 
Évaluation 


de la 
classification 


d'images satellite 


Classification supervisée 


Création du jeu Algorithme de classification 
d'entraînement (p. ex. Random Forest) 


Image radar Sentinel-1 


Figure I.2. Principes généraux des approches de classification supervisée et non supervisée 
d'images satellite optiques et radar. 
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Random Forest, Support-Vector Machine, etc.) traduisent l'information contenue 
dans chaque pixel ou objet en une classe thématique qui décrit l’image étudiée. 


Diffusion des images satellite 


Il existe actuellement une très grande diversité de portails qui mettent à disposition 
des images satellite optiques et radar. C'est le cas pour les capteurs Sentinel-1 et 2 de 
l'Agence spatiale européenne dans le cadre du programme Copernicus’. 


On trouve également un large éventail de logiciels, outils et algorithmes de traitement 
des images satellite qui sont libres et qui participent a la vulgarisation de cette techno- 
logie. Les informations extraites des images de télédétection peuvent par la suite 
être combinées et analysées avec d’autres données spatialisées à l’aide d'un système 
d'information géographique (SIG). 


» Introduction aux systèmes d’information géographique 


Les systèmes d’information géographique (SIG) sont des outils informatiques permet- 
tant l'acquisition, le stockage, la mise à jour, l'intégration, l'analyse, la visualisation et la 
restitution de données numériques géoréférencées (c'est-à-dire que l’on peut localiser 
dans l’espace par leurs coordonnées géographiques). Ils permettent de manipuler et de 
traiter des données spatialisées de natures et de sources différentes. 


L'organisation des données géoréférencées (ou à référence spatiale) dans un SIG repose 
sur le principe suivant : chaque type d'objet (la végétation, les cours d’eau, les villes ou 
les résultats de piégeage des moustiques) est représenté par une couche d’information 
différente (figure L.3). La superposition de ces couches par correspondance spatiale 
permet de visualiser et d'analyser chaque couche d’information séparément (requête 
horizontale telle que : Quelles sont les espèces de moustiques observées et en quelles 


Lieu de collecte entomologique 
(format vecteur — points) 


Réseau hydrographique 
(format vecteur — lignes) 


Occupation du sol 
(format vecteur — polygones) 


Altitude 
(format raster) 


Figure I.3. Principe d'organisation d'un système d’information géographique — couches d'infor- 
mation thématiques et modes de représentation. 


5. https://www.copernicus.eu/fr/propos-de-copernicus/infrastructure/decouvrez-nos-satellites 
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abondances?) et les relations entre des couches d’information différentes (requête 
verticale telle que : Dans quel type d'occupation du sol observe-t-on les abondances de 
moustiques les plus élevées ?). 


Dans une même couche d’information spatiale, les objets sont de même nature et 
représentés selon deux modes possibles : 

— le mode vectoriel ou «vecteur» : dans ce mode, chaque objet est représenté sous 
forme de polygones (p. ex. une parcelle de végétation), de lignes (p. ex. une route ou 
une rivière) ou de points (p. ex. localisation d’un site de piégeage) [figure I.3]. Le format 
de fichier vecteur le plus répandu est le format «.shp » (shapefile) ; 

— le mode matriciel ou «raster » : dans ce mode, l'information spatiale est représentée 
sous forme d’une image (ou grille) composée de cellules de même taille appelées pixels 
(comme sur une image satellite). Le format de fichier raster le plus répandu est le 
format «tif». 


Dans les deux cas, l'information géographique est couplée avec une information 
thématique renseignant sur les propriétés de l’objet. En mode vecteur, cette informa- 
tion thématique est stockée dans une table attributaire associée (par exemple, pour 
les résultats d’un piégeage entomologique, représenté sous forme de point, la table 
associée listera : date du relevé, espèces capturées, abondances, etc.). En mode raster, 
la valeur du pixel contient l'information représentée (pour une image satellite multi- 
spectrale, la valeur des pixels sera la valeur de réflectance mesurée par le capteur). 
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Informations spatiales pour 
la surveillance des moustiques vecteurs 
et des maladies associées 


La première partie de cet ouvrage s'intéresse à l'identification des différentes variables 
environnementales, climatiques et démographiques qui influencent la présence et la 
dynamique des populations de moustiques, avec un focus particulier sur les images 
satellite et leur apport pour la thématique des maladies à transmission vectorielle. 


Le premier chapitre de cette partie recense ces différentes variables ainsi que les données 
et méthodes de télédétection satellite qui permettent d’y avoir accès. Les chapitres 
suivants présentent quant à eux différentes approches reposant sur l'imagerie satellite 
pour l'extraction de ces variables : l’utilisation d'indices spectraux pour l’eau et la végé- 
tation (chapitre 2), l'étude des températures de l'air (chapitre 3), la caractérisation des 
populations humaines (chapitre 4) et enfin l'apport des approches basées sur la texture 
dans les images pour la caractérisation de l'environnement urbain (chapitre 5). 
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Chapitre 1 


Liens entre moustiques vecteurs 
et environnement : apport des méthodes 
de télédétection satellite 


RENAUD MARTI, CLAIRE TEILLET, 
HOBINIAINA ANTHONIO RAKOTOARISON, FLORENCE FOURNET 


Certains traits biologiques sont communs aux près de 3600 espèces de moustiques 
(Diptera : Culicidae) répertoriées dans le monde et conditionnent une grande partie 
des relations avec l’environnement et différents facteurs météorologiques. En tant 
qu'ectothermes, les moustiques ont une température interne contrôlée par la tempé- 
rature de leur environnement, qui conditionne ainsi leur physiologie, comportement, 
écologie et in fine leur survie (Reinhold et al., 2018). Les moustiques ont quatre prin- 
cipaux stades de vie (figure 1.1). Les trois premiers ont une forme aquatique : l'œuf, la 
larve (avec 4 stades larvaires entrecoupés chacun d’une mue) et la nymphe. Le dernier 
stade, l'adulte, associé à la période de reproduction et de dispersion, a une forme 
aérienne. Les stades aquatiques se développent dans des collections d'eau appelées 
«gîtes larvaires», dont les caractéristiques varient selon les espèces de moustiques. 
En influant sur le développement de chacun des stades, de leurs transitions et de leurs 
mortalités associées, les facteurs météorologiques, tels que la température, les précipi- 
tations ou l'humidité, sont des déterminants clefs de la distribution et de la dynamique 
des populations de moustiques. 


Seule la femelle des moustiques est hématophage et c’est au cours d'un repas de sang 
qu'elle peut être infectée par un agent pathogène (virus, bactérie ou parasite) et, après 
multiplication de ce dernier, le transmettre à un nouvel hôte : le moustique est alors 
appelé «vecteur ». En plus de leur influence sur l'écologie du moustique, les variations 
météorologiques impactent l'ensemble du cycle de transmission des agents pathogènes 
et de leur développement, et conditionnent les interactions entre populations humaines 
et moustiques vecteurs (Morin et al., 2013; Stresman, 2010). Un cycle de développement 
plus rapide des moustiques et une durée de survie accrue accélèrent leur taux de repro- 
duction et favorisent la prise d'un nombre croissant de repas sanguins. L'augmentation 
de la température ambiante est également associée à un taux plus rapide de réplica- 
tion virale dans le vecteur ou à une accélération du cycle parasitaire, et donc à une 
période d’incubation extrinsèque plus courte. Si l'on considère une augmentation de 
quelques degrés, compatible avec la biologie actuelle du moustique, l'ensemble de ces 
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phénomènes tend à favoriser la transmission des agents pathogènes. Les effets du chan- 
gement climatique en cours restent cependant difficiles à évaluer, car ils impactent de 
multiples facteurs et processus à différentes échelles et s’accompagnent d’autres modi- 
fications majeures conjointes, comme des processus de coévolution induits entre les 
populations d’arbovirus ou de parasites, d’arthropodes vecteurs et d'hôtes vertébrés. 
Par ailleurs, les microclimats locaux, auxquels sont exposés les moustiques, peuvent 
différer considérablement des mesures macroclimatiques, en particulier dans les 
environnements urbains hétérogènes (Wimberly et al., 2020). 


Si la plupart des arbovirus ou parasites transmis par des moustiques sont respon- 
sables de fardeaux bien plus importants dans les milieux tropicaux et subtropicaux, 
les pays des latitudes tempérées sont aussi affectés par certaines maladies vectorielles. 
Avec l’intensification des échanges internationaux et les effets des changements envi- 
ronnementaux et climatiques, les risques potentiels ou avérés de transmission ou de 
réémergence imposent d’adapter les stratégies de surveillance. Les espéces apparte- 
nant aux trois genres de moustiques Anopheles, Aedes et Culex présentent un intérét 
majeur dans l'étude des maladies à transmission vectorielle (MTV). Les anophèles 
ont une répartition très vaste, quasiment mondiale à l'exception des zones polaires, 
et seule la distinction entre espèces anophéliennes vectrices et non vectrices permet 
de comprendre finement la géographie de la transmission du paludisme (Duvallet 
et al., 2017). Cette distinction est parfois rendue délicate par la présence, au sein 
d'un même complexe, d'espèces dites jumelles, qui ne présentent pas de différences 
morphologiques (impossibilité de les distinguer à l'œil nu), mais des caractéristiques 
biologiques, écologiques et génétiques bien distinctes. Considéré comme le vecteur 
principal du paludisme à l'échelle mondiale, Anopheles gambiae désigne en réalité un 
complexe d'espèces, dont la distinction impose l’utilisation de méthodes moléculaires 
(Coetzee et al., 2013; Davidson, 1964). Sur les quelque 540 espèces du genre Anopheles 
recensées (Duvallet et al, 2017), une soixantaine ont été identifiées comme étant 
vectrices du paludisme, c’est-à-dire capables de transmettre des parasites du genre 
Plasmodium (Manguin et al., 2008). L'Afrique subsaharienne abrite 85 à 95% des cas 
de paludisme. Les principales espèces de moustiques vecteurs responsables de cette 
pathologie qu'on y trouve sont An. gambiae, Anopheles arabiensis et Anopheles coluzzii 
du complexe d'espèces An. gambiae et Anopheles funestus, du groupe d'espèces 
An. funestus (Dahan-Moss et al., 2020). En dehors de la transmission des parasites 
de genre Plasmodium, les anophéles sont également vecteurs de filaires et du virus 
O'nyong-nyong, un arbovirus d'intérêt en santé publique, transmis notamment par 
An. gambiae et An. funestus, et décrit pour la première fois en Ouganda en 1959. 


Les espèces Aedes aegypti et Aedes albopictus, du genre Aedes dans la classifica- 
tion historique et Stegomyia dans la nouvelle classification! proposée par Harbach et 
Horward (2007), sont originellement associées aux climats de la région intertropicale, 
chauds et humides. Ces espèces sont impliquées dans la transmission de plusieurs flavi- 
virus, responsables de maladies comme la dengue, la maladie à virus Zika et la fièvre 
jaune (Afrique et Amérique intertropicales), et de l’alphavirus causant le chikungunya. 
Si Ae. aegypti, caractérisé par une meilleure compétence vectorielle, est le moustique 
vecteur principal de ces arbovirus à l'échelle mondiale, Ae. albopictus est un vecteur 


1. https://mosquito-taxonomic-inventory.myspecies.info/ 
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secondaire important, notamment en Asie et dans les îles de l'océan Indien. Son aire de 
distribution est en extension ces dernières années dans des régions aux latitudes plus 
élevées, comme en Amérique du Nord ou dans le Sud de l’Europe (Kraemer et al., 2019). 


Parmi les espèces de l'assemblage Culex pipiens, les espèces Cx. pipiens pipiens et 
Cx. quinquefasciatus présentent la répartition la plus étendue à l'échelle mondiale et 
sont observées dans toutes les régions tempérées et tropicales urbaines et suburbaines, 
où elles sont souvent identifiées comme les principaux vecteurs d’arbovirus d'intérêt 
en santé vétérinaire et en santé publique. Le genre Culex est notamment associé à 
la transmission de deux flavivirus, l’un responsable de l’encéphalite japonaise (avec 
comme principal vecteur Culex tritaeniorhynchus) observée en Asie, et l'autre de la 
fièvre à virus West Nile (avec comme vecteur en particulier Culex pipiens) présente en 
Europe, Afrique, Asie et Amérique du Nord et du Sud. Les genres Aedes (notamment 
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Figure 1.1. Cycle de vie des moustiques femelles illustré avec les différents stades aériens et 
aquatiques de l'espèce Aedes albopictus. Photographies des stades larvaires : © Nicolas Henon — 
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l'espèce Aedes [Aedimorphus] vexans) et Culex (notamment l'espèce Culex poicilipes) 
peuvent être associés à la transmission de la fièvre de la vallée du Rift, une zoonose 
virale observée en Afrique et causée par un phlebovirus. Enfin, les trois genres 
Anopheles, Aedes et Culex peuvent étre impliqués dans la transmission de la filariose 
lymphatique, une maladie tropicale parasitaire causée par des nématodes de la famille 
des Filariidae. 


» Liens entre moustiques vecteurs et environnement 


Le système vectoriel, associant agents infectieux, moustiques vecteurs et hôtes verté- 
brés, s'inscrit dans un environnement. Défini à la fois par des facteurs biotiques 
(couvert végétal, présence de gîtes naturels, de sources de nourriture, de prédateurs, 
etc.) et abiotiques (pluie, température, humidité, rayonnement, etc.), l'environnement 
influe sur les densités des différentes populations composant le système vectoriel et 
sur la fréquence de leurs interactions. 


Influence des facteurs météorologiques sur Le cycle de vie 


Considérations générales 


La température influence de manière non linéaire les différents stades de dévelop- 
pement des moustiques, avec l'observation d'effets de seuil. Une température trop 
faible de l'eau des gîtes larvaires limite voire stoppe le développement des stades aqua- 
tiques, tandis qu'une augmentation des températures jusqu’à une autre valeur seuil a 
tendance à l’accélérer. L'augmentation de la température de Pair, à condition que celle- 
ci reste dans un intervalle optimal, favorise également la survie et l’activité du stade 
adulte, dont la prise de repas de sang. 


Chaque espèce de moustique pond ses œufs dans des gîtes larvaires spécifiques, qui 
peuvent être naturels (flaques, trous dans les arbres, rochers, coquilles, plans ou cours 
d’eau, etc.) ou artificiels (contenants de tailles et matières diverses, vides sanitaires, etc.), 
remplis d'eau douce, saumâtre ou salée et plus ou moins riches en matière organique. 
D'une manière générale, la disponibilité en gîtes larvaires est déterminée par les niveaux 
de précipitations (pluie), qui peuvent suivre un cycle saisonnier selon la latitude et le 
climat. La mise en eau des gîtes peut également être artificielle, selon le contexte et les 
pratiques humaines associées : arrosage des plantes en intérieur, des jardins en extérieur, 
irrigation et mise en eau des cultures agricoles (p. ex. rizière). La température de lair 
influe sur le maintien en eau des gîtes via le phénomène d'évapotranspiration. Enfin, les 
précipitations favorisent le développement de la végétation, source de sucre et de lieux 
de repos pour les moustiques adultes (Sallam et al., 2017). 


Focus sur des espèces d'intérêt médical 


Les œufs des moustiques du genre Anopheles sont pondus isolément à la surface de 
l'eau. Chez An. gambiae, une température optimale pour les œufs semble comprise 
entre 24 et 30°C, avec des durées de développement inférieures à 7 jours, tandis que 
des températures inférieures (12°C) et supérieures (48°C) réduisent nettement la 
viabilité des œufs (Impoinvil et al., 2007). Les œufs d’anophéles sont pondus isolément 
à la surface de l’eau et résistent très peu à la dessiccation, avec une très nette réduction 
du taux d’éclosion sur substrats secs (Duvallet ef al., 2017). A contrario, en conditions 
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humides, des œufs d'An. gambiae et d'An. funestus ont pu éclore jusqu’à 12 et 10 jours 
respectivement après la ponte. Le développement des stades larvaires d'An. gambiae 
est également favorisé par l'augmentation de la température avec l'observation d'effets 
de seuils et de non-linéarités : développement non linéaire à partir d’un seuil initial à 
16°C, développement proportionnel entre 22 et 28°C (optimum), déclin non linéaire 
jusqu’à un seuil maximal de 34°C, avec l'observation d'effets délétères au-delà (Bayoh 
et Lindsay, 2003). La plupart des espèces d’anophéles subissent une saison défavo- 
rable — froide en zone tempérée ou sèche en zone tropicale —, qui peut conduire à 
différentes stratégies de survie de la femelle adulte, selon l'espèce et le contexte biogéo- 
graphique : diapause, dormance, migration, estivation ou discordance gonotrophique 
(repas de sang sans production d'œufs) [Duvallet et al., 2017]. 


Les œufs des moustiques du genre Aedes sont pondus isolément sur un support à 
proximité de l'eau. Durant les périodes prolongées sans précipitation, ils résistent remar- 
quablement bien à la dessiccation (Duvallet et al., 2017). Par ailleurs, des conditions de 
sécheresse prolongées peuvent, dans certains contextes, conduire à un accroissement 
du stockage en eau sous la forme de récipients ou conteneurs domestiques, susceptibles 
d'offrir de nouveaux gîtes de ponte (Pontes et al., 2000). A contrario, de trop fortes préci- 
pitations sur une courte période peuvent engendrer un phénomène de lessivage des 
gîtes, conduisant à une surmortalité des stades aquatiques (Dieng et al., 2012). Le couple 
température et pourcentage d'humidité détermine fortement l’activité reproductrice et 
la survie des Aedes, particulièrement sensibles à de faibles variations de ces deux para- 
mètres, mutuellement influents. Pour une humidité fixée à 80%, les femelles Ae. aegypti 
semblent ainsi survivre plus longtemps et produire davantage d'œufs à 25°C qu'à 35°C, 
tandis qu'une humidité à 60% semble réduire drastiquement voire inhiber le mécanisme 
doviposition (Costa et al., 2010). En contexte tropical, marqué par des alternances de 
fortes périodes pluvieuses et de périodes plus sèches, les variations de température et 
d'humidité, même d’amplitudes limitées, peuvent conduire à une augmentation ou à une 
réduction saisonnière de la densité de population d’Aedes (Nasir et al., 2017). En contexte 
tempéré, à la suite du développement d'un mécanisme adaptatif de certaines souches, 
la femelle Ae albopictus a la capacité de pondre des œufs diapausants, qui résistent à 
des conditions climatiques défavorables (faibles températures), puis reprennent leur 
développement jusqu’à l'éclosion dès que les conditions plus favorables sont à nouveau 
réunies. Si les œufs d'Ae, aegypti ont montré une bonne tolérance au froid en labora- 
toire, voire sur le terrain (Kramer et al., 2020), le mécanisme de diapause n’a pas encore 
été observé en milieu naturel, ce qui pourrait expliquer la difficulté de son implantation 
pérenne dans le continent européen. De récents résultats expérimentaux obtenus pour 
certaines souches d’Ae. aegypti associées à une région tempérée (Buenos Aires, Argen- 
tine) semblent cependant montrer une sensibilité des formes adultes à la réduction de la 
photopériode conduisant à une diminution de la fréquence d’éclosion des œufs (Fischer 
et al., 2019). Si elle était confirmée, cette capacité d’Ae. aegypti à suspendre l'éclosion 
des œufs en réponse à une photopériode de jours plus courts pourrait favoriser son 
expansion vers des régions aux hivers plus froids. 


Culex pipiens et Cx. quinquefasciatus sont présents dans la plupart des zones habi- 
tées du monde et sont souvent étroitement associés aux humains, ce qui leur vaut les 
noms de moustiques domestiques du Nord et du Sud (Farajollahi et al., 2011). Les 
femelles pondent leurs œufs de façon regroupée en formant une barquette (ou radeau 
ou nacelle) flottant sur l'eau, au sein de latrines, puisards ou caniveaux. La corrélation 
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positive entre la température et le taux de développement de ces espèces apparaît 
moins marquée à partir de 24°C, avec généralement une augmentation concomi- 
tante de la mortalité adulte pour des températures plus élevées (Ciota et al., 2014). 
La distribution mondiale des moustiques du complexe Cx. pipiens s'explique par leur 
adaptation aux environnements modifiés, voire altérés par les humains, ces derniers 
favorisant leur dispersion, ainsi que par leur mode d'alimentation mixte ciblant les 
oiseaux et les mammifères (y compris les humains). Les pluies abondantes, l'irrigation 
des cultures, le mauvais entretien des réseaux d'assainissement des villes ainsi que 
des températures élevées favorisent la pullulation des Culex : ces moustiques affec- 
tionnent les collections d’eau souillée riche en matière organique (Darriet, 2014). Sous 
les climats tempérés, les adultes femelles de Cx. p. pipiens passent l’hiver dans les 
caves et les égouts (Darriet, 2014). 


Interactions des moustiques vecteurs avec leur biotope 


Considérations générales 


Dans un environnement favorable, les moustiques présentent une stratégie d’exploi- 
tation maximale des ressources disponibles, accompagnée d’un taux de natalité élevé, 
d'un temps de génération court et d'une bonne capacité de dispersion (Duvallet et al., 
2017). La distribution spatiale est alors conditionnée par la structuration de l'environ- 
nement, avec en premier lieu la disponibilité de collections d’eau favorables à la ponte 
des femelles gravides et au développement des stades larvaires, et par la présence 
d'hôtes (animaux et humains) pour les repas de sang. Du fait de leur remarquable 
plasticité écologique, certaines espèces de moustiques des genres Anopheles, Aedes et 
Culex se sont parfaitement adaptées aux environnements anthropisés et présentent 
un comportement anthropophile qui favorise d'autant plus la transmission d'agents 
pathogènes (Cohuet et al., 2010). Ces espèces profitent notamment des aménagements 
urbains ou périurbains, et des activités humaines associées, pour adapter leurs sites 
de reproduction, de ponte et de repos. Elles peuvent cibler les humains, soit comme 
principale source de repas de sang, soit de manière plus opportuniste. 


Environnement des stades aquatiques 


Si les trois stades juvéniles — œuf, larve et nymphe — sont aquatiques, les larves 
sont les seules à se nourrir et la majorité des effets qui affectent leur développement 
et/ou leur taux de mortalité est liée à leur densité dans le milieu (Beck-Johnson et al., 
2013). Le développement des larves passe par quatre stades (L1 à L4), qui sont les seuls 
moments où les moustiques croissent de manière continue et progressive en taille. 
La durée du développement larvaire est très variable (de quelques jours à quelques 
mois voire plus si les larves passent en diapause larvaire) et dépend fortement des 
conditions climatiques, écologiques et environnementales. Les conditions du déve- 
loppement larvaire sont définies par la qualité et la quantité de nourriture (levures, 
bactéries, microplancton, microalgues, grains de pollen, etc.), la température, les 
compétitions intra- et interspécifique et la présence de sources de stress biotiques ou 
abiotiques. Selon l'espèce et la période considérées, la disponibilité en gîtes larvaires 
devient donc un facteur clef de la dynamique de population. La durée de vie de la 
nymphe, déterminée par la température, est comparativement courte, de l’ordre de 
quelques jours, au cours desquels celle-ci ne s’alimente pas. 
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La majorité des espèces d’anophéles nécessite des gîtes d'eau douce, bien que 
certaines présentent une forte tolérance à la salinité. Il est à noter quun phéno- 
mène d'adaptation à l’eau polluée est observé dans les contextes urbains pour 
certaines espèces, comme An. coluzzii qui a développé une tolérance à l’ammoniac, 
l’un des principaux polluants des gîtes aquatiques urbains (Duvallet et al., 2017). 
Pour Ae. albopictus et Ae. aegypti, tous les types de récipients artificiels ou natu- 
rels de petites dimensions et susceptibles d'être mis en eau constituent des gîtes de 
ponte potentiels (Estallo et al., 2008). Les œufs des moustiques du genre Aedes sont 
pondus isolément sur un support à proximité de l’eau. Les gîtes larvaires des Culex 
sont nombreux et hétéroclites (puisards, caniveaux, etc.), largement disséminés au 
sein des eaux polluées et des effluents de la ville, où les larves se nourrissent de la 
matière organique abondante (Darriet, 2014). 


Environnement des stades aériens et du cycle gonotrophique 


L'émergence de l'adulte représente une courte période (dizaine ou quinzaine de 
minutes) de très grande vulnérabilité à la prédation et aux modifications environ- 
nementales. Certaines espèces, dites « eurygames », ont besoin de vastes espaces pour 
s'accoupler (p. ex. Ae. caspius, Ae. detritus, Cx. pipiens pipiens), tandis que d’autres, 
sténogames, peuvent s’accoupler au sein de petits espaces (Ae. albopictus, Cx. pipiens 
molestus). Après l’accouplement, la femelle néonate (nullipare) cherche un hôte 
pour un premier repas de sang, nécessaire a la production des ceufs dans le cas des 
espèces dites « anautogénes». Au cours de leurs quelques semaines d’existence (entre 
3 et 5 semaines pour la majorité des espéces en zone tropicale), le comportement 
des femelles est dominé par la réalisation d’une dizaine de cycles gonotrophiques 
(tropho- sang, -gonique développement ovarien). Chaque cycle peut étre schémati- 
quement décomposé en trois étapes principales, d'une durée totale de 3 à 5 jours selon 
les espèces et les conditions climatiques (Duvallet et al., 2017; figures 1.1 et 1.2): 

— recherche et choix de l'hôte vertébré par la femelle à jeun; 

— repos pour la digestion du sang et la maturation des ceufs; 

— recherche du site de ponte et oviposition par la femelle gravide. 


Pour Anopheles, le repas de sang se déroule au crépuscule ou de nuit et humain peut 
être identifié comme une source prioritaire ou opportune, selon l’espéce. Le repas 
peut être pris à l’intérieur (endo-) ou à l'extérieur (exophagie) des habitations. Pour 
Aedes, le repas de sang est diurne. Il s’effectue sur l'humain en priorité et à l’intérieur 
pour Ae. aegypti et se répartit selon les mammifères présents, plus généralement à 
l'extérieur des habitations pour Ae. albopictus. Culex pique la nuit, les mammifères 
en général (y compris l’homme), ciblant davantage les oiseaux (bird biting mosquito) 
pour Cx. pipiens pipiens en milieu rural (tempéré) et l'humain en milieu urbain, 
avec une endophilie (préférence marquée pour l’intérieur des maisons) prononcée 
pour Cx. pipiens molestus (milieu tempéré) et Cx. quinquefasciatus (milieu tropical) 
[Duvallet et al., 2017; Farajollahi et al., 2011]. 


Sauf à l'émergence, les adultes des deux sexes ne boivent pas d'eau, mais prennent 
régulièrement des repas de jus sucrés sur divers végétaux, incluant des nectars 
floraux. Les moustiques mâles dépendent exclusivement de ces sources d’alimenta- 
tion. La préférence pour les plantes varie selon les espèces de moustiques, les habitats 
géographiques et la disponibilité saisonnière. Lorsque les sources de sucre sont rares, 
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les moustiques femelles de certaines espèces, comme An. gambiae ou Ae. aegypti, 
peuvent compenser en prenant des repas sanguins plus importants et plus fréquents 
(Barredo et Degennaro, 2020). Par ailleurs, certaines espèces de moustiques peuvent 
également se nourrir de suspensions riches en sucre et de fruits en décomposition, 
potentiellement abondants dans les déchets produits par les populations humaines. 
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Figure 1.2. Inscription du cycle de vie des moustiques Aedes albopictus dans un biotope de 
type urbain ou périurbain. 


Outre leur capacité de dispersion active leur permettant d'assurer leur cycle de développement, une 
dispersion passive, à l’aide des moyens de transport humains (véhicules individuels, transports publics, 
etc.), peut assurer leur dissémination sur de plus grandes distances. Lors de leurs déplacements (recherche 
d'hôtes pour le repas de sang ou de gîtes de ponte), les femelles moustiques peuvent être exposées 
à une surmortalité (symbolisée par un petit crâne rouge). Les principales sources de caractérisation de 
l'environnement par télédétection sont indiquées dans la partie supérieure de la figure. 
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En ville, les ressources nécessaires au cycle de vie des moustiques adaptés aux 
habitats artificiels sont disponibles sur de faibles étendues (partenaires de repro- 
duction, sources d'alimentation, sites de ponte et de repos, et hôtes pour les repas 
de sang). Les capacités de dispersion active de Ae. aegypti ou Ae. albopictus sont 
par ailleurs considérées comme limitées, généralement inférieures 4 1km. Leur 
transport passif peut être en revanche associé à de longues distances soit sous la 
forme d'œufs résistants à la dessiccation (transport par bateaux de pneus usés ou de 
certains végétaux comme les lucky bamboo), soit sous leur forme adulte au travers de 
moyens de transport par route ou rail, dont l'importance du trafic favorise la proba- 
bilité d'implantation sur de nouveaux territoires (Eritja et al., 2017). La distance 
de dispersion active pour An. gambiae est considérée comme plus importante, de 
l'ordre de 2 à 3km par génération (Duvallet et al., 2017). Sa dispersion passive est en 
général assurée par des moyens de transport humains (avion, bateau, voiture, etc.), 
véhiculant sa forme adulte. 


» Description de l'environnement par des approches 
de télédétection satellite 


Intérêt de la télédétection 


La télédétection satellite a une couverture mondiale et les approches qui y sont 
liées sont reproductibles de manière dynamique, en différents lieux et à différentes 
résolutions spatiales ou temporelles. Elles peuvent ainsi s'avérer efficaces pour estimer 
des variables environnementales essentielles comme certains paramètres météo- 
rologiques (température, précipitation, répartition en vapeur d’eau, vent, etc.) ou bien 
pour établir un premier niveau de description de l'environnement biophysique au 
travers d'indices spectraux ou de produits plus élaborés comme l'occupation ou l'usage 
du sol (Parselia et al., 2019; voir Introduction générale et chapitres 2 à 5). L'intérêt 
pour ce type de techniques en écologie de la santé est régulièrement renouvelé au 
gré des innovations technologiques et de la disponibilité de nouveaux capteurs et des 
produits associés (Goetz et al., 2000; Herbreteau et al., 2007). L'avènement récent 
du programme Copernicus? et des satellites Sentinel est à ce titre remarquable, avec 
la mise à disposition de séries temporelles d'images (revisite hebdomadaire) et de 
produits dérivés à une résolution spatiale décamétrique. 


Le croisement des données de télédétection avec d’autres paramètres socio- 
économiques ou d'aménagement au sein des systèmes d’information géographique 
(SIG) permet d'enrichir la représentation de l’environnement sur un territoire, en 
lien avec l'écologie des moustiques (stade larvaire ou adulte), en prenant par exemple 
en compte la gestion de l’eau ou des déchets au sein de la ville (Kolimenakis et al., 
2021). Les produits issus de la télédétection peuvent servir de données d'entrée aux 
différentes approches de modélisation, en permettant d'intégrer une structura- 
tion spatiale de l’environnement dans les modèles fondés sur les observations et la 
recherche de relations statistiques (chapitre 6), les modèles fondés sur les connais- 
sances (chapitre 7), ou encore dans les modèles explicitant les processus et les traits 
biologiques du moustique (chapitre 8 et 9). 


2. https://www.copernicus.eu/ 
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Dans le contexte des maladies à transmission vectorielle humaines comme la dengue, 
la détection et la caractérisation des milieux urbains apparaissent prioritaires 
(chapitres 4 et 5) en raison de l'adaptation de certains vecteurs d'importance médicale 
aux environnements urbains et de l'accroissement global de la population humaine en 
ville, plus nombreuse qu'en milieu rural depuis 2008 (Prasad et al., 2016). Entre 1950 
et 2050, le taux de personnes habitant en ville devrait ainsi passer de deux à sept sur 
dix, avec une dynamique urbaine particulièrement forte attendue en Asie et en Afrique 
pour les prochaines décennies. 


Estimation des variables météorologiques par télédétection 


À la surface de la Terre, les stations météorologiques permettent notamment de 
mesurer la température et l'humidité de lair, la durée d’ensoleillement et le vent, 
variables susceptibles d’influencer le développement et la dispersion des moustiques 
(cf. section précédente). Dans le cas de sites non instrumentés, certaines de ces 
variables peuvent étre estimées indirectement par télédétection satellite (tableaux 1.1 
et 1.2). Contrairement aux satellites géostationnaires (orbite distante de 36000km), 
les satellites héliosynchrones (qui revisitent un point du globe à la même heure solaire 
locale) ont des orbites beaucoup plus proches de la surface terrestre (le plus souvent 
entre 500 et 1000km), qui offrent des résolutions spatiales davantage compatibles 
avec l'étude des phénomènes écologiques. 


Si les capteurs embarqués des satellites météorologiques permettent de reconstruire 
des profils verticaux de température, les satellites destinés à l'étude de l’environ- 
nement permettent de réaliser des mesures directes de la température de brillance 
du sol (tableau 1.1). Les instruments embarqués mesurent le spectre d'émission dans 
la gamme de l’infrarouge thermique (3-14um). A l’aide d’algorithmes spécifiques 
(p. ex. Split Window), les valeurs d'émissivité et de température du sol (LST pour 
Land Surface Temperature) peuvent être démélées et estimées comme des variables 
distinctes. Les produits de température de surface permettent de bien identifier et de 
délimiter la structure spatiale des îlots de chaleur qui se forment au sein des villes. 


Tableau 1.1. Principaux instruments de télédétection satellite actuellement en service, 
associés à la variable «température de surface ». 


Caractéristiques d’acquisition 


Plateformes satellites Emprise Résolution Fréquence Exemples d’utilisation 
et instruments associés d’une image spatiale de revisite dans le contexte des MTV 


ECOSTRESS 384km 70 m 1 jour (Boser et al., 2021) 
(Ecosystem Spaceborne 

Thermal Radiometer 

Experiment) 


Landsat 8/9 TIRS 185km 100 m 16 jours (Ogashawara et al., 2019) 
(Thermal Infrared 
Sensor) 


MODIS MOD11A2 1100km 1km 8 jours (Yue et al., 2018) 
(Moderate 

Resolution Imaging 

Spectroradiometer) 
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Ces îlots sont susceptibles d’influencer l'écologie du moustique, en créant des condi- 
tions locales plus favorables au sein des villes (Huraux et al., 2017). L'association 
d'une température de surface à une température de l’air reste cependant délicate par 
une approche physique et se modélise généralement au travers d’une relation statis- 
tique avec des observations in situ (Boser et al., 2021; Weiss et al., 2014; chapitre 3). 
La résolution spatiale des capteurs permettant d’estimer la température de surface 
est actuellement, au mieux, de l’ordre de la centaine de mètres voire du kilomètre, 
ce qui ne permet pas de travailler à une échelle fine, comme les quartiers d'une ville. 
La future mission franco-indienne Trishna devrait offrir une résolution spatiale à 
ce jour inédite (57m) et permettre de dépasser certaines contraintes d'exploitation 
thématique actuelles (Lagouarde et al., 2018). Dans certains contextes régionaux, 
le produit Sentinel-3 SLSTR LST L2 pourrait s'avérer également intéressant avec 
son emprise de 290 km, sa résolution spatiale de 1 km et une fréquence de revisite 
journalière (Shumilo et al., 2019). 


Concernant la mesure des précipitations, les produits dérivés de la mission Global 
Precipitation Measurement (GPM, depuis 2010) offrent une continuité aux données 
historiques de la mission Tropical Rainfall Measuring Mission (TRMM, 1997-2015). 
Les produits dérivés GSMaP (Global Satellite Mapping of Precipitation) estiment les 
précipitations tous les 3 jours avec une résolution de 10 km (Guilloteau et al., 2014) et 
le produit Integrated Multi-satellite Retrievals for GPM (IMERG) combine les données 
TRMM et GPM pour estimer les précipitations entre 2010 et aujourd’hui a une résolu- 
tion de 10km et une résolution temporelle allant du temps quasi réel au pas journalier 
et mensuel (Tsantalidou et al., 2021). 


Tableau 1.2. Principaux produits de précipitations a l'échelle mondiale intégrant des 
données de télédétection satellite. 


Caractéristiques d’acquisition 


Exemples d'utilisation 


Produits Période Résolution Fréquence 

météorologiques d'acquisition spatiale d'acquisition dansle contexte 
sq q P aq des MTV 

Product 3B43 v6 1997 -2015 30km Moyenne (Scavuzzo et al., 2018) 

Tropical Rainfall mensuelle 

Measuring Mission 

TRMM 

GSMaP/MVK v6 2014 - 10km Moyenne (Guilloteau et al., 2014) 

Global Precipitation en cours à 3 jours 

Measurement Mission 

GPM IMERG 2010 - 10km Moyenne quasi- (Tsantalidou et al., 

Fusion de TRMM en cours temps réel, 2021) 

(2010-2015) et GPM quotidienne 

(2014-present) et mensuelle 

Worldclim Historique, 1km Moyenne (Tsheten et al., 2021) 

Variables climatiques présent, futur mensuelle 


et bioclimatiques 


(T, P, V) ou selon 
la nature de la 
variable générée 
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Par ailleurs, de nombreuses bases de données nationales, régionales ou mondiales offrent 
des données climatiques (température, humidité, précipitation, vent, etc.) sous la forme 
de réanalyses générées par des modèles climatiques intégrant des données d’observa- 
tions satellites et in situ. Des variables bioclimatiques peuvent également être dérivées 
des valeurs mensuelles de température et de précipitations afin de produire des variables 
plus significatives sur le plan biologique (plage diurne moyenne, isothermes, moyennes 
annuelles de températures et de précipitations, etc.). À une résolution proche de 1 km, 
la base de données « WorldClim » fournit des données météorologiques et climatiques 
maillées pour les conditions historiques (proches des conditions actuelles) et futures. 
De telles données ont déjà été exploitées dans le cadre d'études sur l'aire de distribution 
mondiale de certains moustiques vecteurs (Kraemer et al., 2015) ou de la transmission 
de certaines maladies vectorielles (Tsheten et al., 2021). 


Description de l'environnement biophysique par télédétection 


Capteurs optiques et radar utiles dans le contexte des maladies 
à transmission vectorielle 


Selon leurs caractéristiques, les images de télédétection permettent d'identifier des 
objets de quelques centimètres (p. ex. images de drone) ou d'obtenir des images de la 
surface terrestre avec une fréquence d'acquisition quotidienne (comme avec le satel- 
lite MODIS). Dans le cadre du programme Copernicus, l'Agence spatiale européenne 
(ESA) a lancé la constellation de satellites « Sentinel». Mis en orbite respectivement 
le 3 avril 2014 et le 25 avril 2016, Sentinel-1A et Sentinel-1B fournissent des images 
radar (en bande C) a haute résolution (10m), par tous les temps, de jour comme de 
nuit. Lancés respectivement le 22 juin 2015 et le 7 mars 2017, Sentinel-2A et Sentinel- 
2B produisent des images optiques de haute résolution (10 et 20m principalement) 
particulièrement adaptées au suivi des surfaces végétalisées, aquatiques, agricoles 
ou encore des milieux urbains (tableau 1.3). Plusieurs articles de synthèse listent de 
manière assez complète les capteurs mobilisables dans le cadre d’études sur l'écologie 
des moustiques vecteurs et l’épidémiologie des maladies vectorielles, en fonction 
des surfaces à identifier (Catry et al., 2018b; Herbreteau et al., 2018). De nouvelles 
missions devraient permettre d’atteindre des niveaux de description plus fins de l'envi- 
ronnement, avec des résolutions proches de la télédétection aéroportée : 30 cm pour 
Pléiades Neo par exemple (Soubirane, 2019). 


Calculs d’indices spectraux 


Dans le domaine des applications de la télédétection, la communauté scientifique a 
développé de nombreux indices exploitant la combinaison de différentes bandes spec- 
trales pour caractériser les propriétés biophysiques des surfaces imagées (tableau 1.4). 
Associé à l’activité chlorophyllienne, l’indice de végétation par différence normalisée 
(NDVI pour Normalized Difference Vegetation Index) repose sur le comportement 
spectral de la végétation en mesurant la différence entre les valeurs radiométriques 
enregistrées par le canal «proche infrarouge» (forte réflexion de la végétation) et 
le canal «rouge» (absorption). Le NDVI est un indice avec des valeurs continues 
comprises entre —1 (absence de végétation) et +1 (végétation dense avec une forte acti- 
vité chlorophyllienne). D’autres indices de végétation peuvent également étre extraits, 
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Tableau 1.3. Exemples de capteurs de télédétection satellite optiques et radar, mobilisables 
dans le cadre de l'étude de l'écologie des moustiques vecteurs et de la transmission des MTV. 


Caractéristiques d'acquisition 


Mode d'acquisition Emprise Résolution Fréquence Exemples d'utilisation 
et plateforme satellite d’une image spatiale de revisite dans le contexte 
(au nadir) des MTV 

(Série temporelle) 80km 10m 5 jours (Hardy et al., 2019) 
Sentinel-1 
(Série temporelle) 290km 10 ou 20m 5a 10 jours (Tran et al., 2019) 
Sentinel-2 
(Série temporelle) 185km 15m (P) 16 jours (Kofidou et al., 2021) 
Landsat 7 ETM+ 30m (XS) 
Landsat 8 
(Série temporelle) 2230 km 250 m 1 jour (Troyo et al., 2009) 
MODIS 500 m Synthèse a 7 j 
(Programmation) 60 km 1,5 m (P) Program- (Orta-Pineda et al., 2021) 
SPOT 6/7 8m (XS) mation à la 

demande 
(Programmation) 20km 0,7 m (P) Program- (Georganos et al., 2020) 
Pléiades 2,8m (XS) mation ala 

demande 


Les lettres (P) et (XS) de la colonne « Résolution spatiale» désignent respectivement les caractéristiques 
panchromatique ou multispectrale des images (un seul canal d’acquisition à large spectre ou des canaux 
spécifiques à certaines longueurs d'onde, p. ex. rouge, vert, bleu, proche infrarouge, etc.). 


Tableau 1.4. Principaux indices dérivés des images de télédétection pour caractériser la 
présence de végétation, d'eau, de bâtiments. 


Caractéristiques de production 


Indice spectral Images Surfaces Méthode Exemples d'utilisation 
satellite caractérisées de calcul* dans le contexte 
sources des MTV 

NDVI Multi- Surfaces Différence (Richman et al., 2018) 

(Normalized Difference spectrales  végétalisées normalisée 

Vegetation Index) des canaux PIR 

et Rouge 

MNDWI Multi- Surfaces Différence (Malahlela et al., 2018) 

(Modified Normalized ` spectrales en eau normalisée 

Difference Water Index) entre canaux 

(Vert ou PIR) 
et (SWIR 
ou MIR) 

NDBI Multi- Surfaces Différence (Demets et al., 2020) 

(Normalized Difference  spectrales construites normalisée 

Built-up Index) (batiments) entre canaux 

SWIR et PIR 


* PIR: proche infrarouge ; SWIR : Short Wave Infrared / infrarouge à ondes courtes ; MIR : moyen infrarouge. 
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comme lEVI (Enhanced Vegetation Index), qui permet d'estimer les variations de 
structures de la canopée, ou le SAVI (Soil Adjusted Vegetation Index) qui permet 
de prendre en compte la contribution due au sol. Les indices MNDWI (Modified 
Normalized Difference Water Index) ou NDWI (Normalized Difference Water Index) 
permettent d'évaluer si la surface est en eau, tandis que le NDBI (Normalized Diffe- 
rence Built-up Index) ou le BI (Brightness Index) permettent de détecter si les surfaces 
sont de type «bâti». Un seuillage de chaque indice permet d'extraire les surfaces consi- 
dérées comme végétalisées (ND VD), en eau (MNDW/), ou bâties (NDBI) [chapitre 2]. 


Produits cartographiques mondiaux et régionaux disponibles 


Des agences spatiales nationales et des grands centres scientifiques (NASA, ESA, DLR, 
JRC ou le pôle Theia?) mettent régulièrement à disposition des produits réalisés à l'échelle 
mondiale ou régionale (tableau 1.5). En dépit de certaines limites inhérentes à ce type 
de production (qualités inégales, dates parfois anciennes), ces produits peuvent s’avérer 
particulièrement utiles dans des contextes pauvres en données d'observation in situ. 


Tableau 1.5. Produits cartographiques pertinents au regard de l'écologie de certains 
vecteurs : emprises de bâtiments, topographie, occupation du sol. 


Caractéristiques de production 


Mot produit Types Plateformes Méthodes Publications 
denvironnement d'acquisition associées 
GHSL Emprise Sentinel-1 et2 Chaîne (Pesaresi et al., 2013) 
des bâtiments MASADA 
GUF Emprise TerraSAR-X Texture et (Esch et al., 2017) 
des batiments Support 
Vector Data 
Description 
SRTM Topographie Navette Interférométrie (Slater et al., 2006) 
Endeavour radar 
OSO Occupation Sentinel-2 Chaine iota? (Inglada et al., 2017) 
(France du sol 


métropolitaine) 


Production de cartographies et de classifications ad hoc 


Dans de nombreux cas, les produits mondiaux ou régionaux de classification d’occu- 
pation du sol peuvent s'avérer non adaptés pour représenter l'écologie des vecteurs et 
l'environnement de transmission des maladies vectorielles : sémantiques des classes 
incomplètes ou peu appropriées, résolutions spatiales inadaptées, produits devenus 
obsolètes. L'utilisateur disposant des compétences techniques adéquates peut alors 
générer ses propres produits cartographiques soit à l’aide de chaînes de traitement 
open source (tableau 1.6), soit en construisant ses propres outils à l’aide de briques 
logicielles spécifiques au traitement d'image (comme Orfeo Toolbox* ou SNAP’). 


3. https://www.theia-land.fr/pole-theia-2/ 
4. https://www.orfeo-toolbox.org/ 
5. https://step.esa.int/main/download/snap-download/ 
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Tableau 1.6. Exemples de chaînes de traitement open source permettant de caractériser 
des composantes environnementales d'intérêt écologique pour l'étude des moustiques 
vecteurs et des MTV, à partir d'images satellite. 


Caractéristiques de production 


Nom de la chaîne Images satellite Types r Publications 
j SC Méthodes Ze 
de traitements sources d'environnement associées 
Iota? Sentinel-1 et 2 Généraliste Approche (Inglada et al., 
orientée pixel 2017) 
Moringa Sentinel-2 Généraliste Approche (Dupuy et al., 
et THRS orientée objet 2020) 
(p. ex. SPOT 6/7, 
Pléiades) 
Urba-Opt Sentinel-2 Urbain Approche (Puissant et al., 
orientée objet 2019) 
Fototex Pléiades, Principalement Calcul de (Teillet et al., 
Sentinel-2 urbain, mais transformée 2021) 
exploitable de Fourier par 
sur la végétation fenêtre glissante 
(indice de texture) 
WaterDetect Sentinel-2 Surfaces en eau Technique (Cordeiro et al., 
de clustering 2021) 
Classificateur 
Naive Bayes 


En lien avec l'écologie du vecteur et son interaction avec l’environnement (cf. sections 
précédentes), l'utilisateur restitue l'environnement biophysique sous un format numé- 
rique, compatible avec les SIG. Cette représentation se fait sous la forme d'objets 
géographiques extraits des images ou bien de classes d'occupation du sol associées à la 
couverture biophysique (forêt, surface de type agricole, etc.) et/ou de classes d'usage 
du sol associées à la fonctionnalité de l’espace (p. ex. zone d'exploitation forestière, 
culture rizicole). Ces classifications ou « extractions » d'objets doivent être réalisées en 
cohérence avec les images disponibles (résolution, emprise, date), la méthode retenue 
(classification orientée objet ou pixel, supervisée ou non supervisée), le choix d’un jeu 
de paramètres appropriés, la disponibilité de données d'apprentissage de qualité pour 
les méthodes supervisées, et la connaissance de l'écologie du vecteur (espèce) dans le 
contexte biogéographique de l'étude (climat, saisonnalité, habitats). 


Quelques grandes typologies géographiques (qu’il convient d'adapter selon l'espèce et 
le cycle de transmission) permettent de structurer et de caractériser le paysage dans 
le contexte générique de l'étude des moustiques vecteurs et des maladies vectorielles : 
— l'emprise des bâtiments, qui renseigne sur la présence humaine et structure fortement 
l'ensemble du territoire (Troyo et al., 2009). La superficie des bâtiments et leur hauteur 
(nombre d'étages) peuvent permettre de construire un indicateur de la densité humaine 
dans le cadre de l'exposition à la piqûre (cf. chapitre 4; Georganos et al., 2020); 

— la typologie fonctionnelle associée aux emprises des bâtiments, qui caractérise 
parfois les habitudes et les comportements des populations humaines et module 
différents compartiments du risque de transmission : zones de résidences indivi- 
duelles avec jardins, immeubles d’habitation à plusieurs étages, zones industrielles 
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ou commerciales (p. ex. marchés) ou encore habitats spontanés ou en construction 
(Flamand, 2015). Le type d'activité humaine associé à une construction permet d’iden- 
tifier des caractéristiques écologiques pertinentes pour le moustique vecteur (p. ex. 
densité de gîtes larvaires pour la phase d’oviposition) ou de potentielles interactions 
entre moustiques vecteurs et humains (repas de sang selon le taux d'occupation du 
lieu, en fonction des heures de la journée); 

— les voies de communication : routes de différents gabarits, qui marquent la mobilité 
humaine et peuvent être synonymes de barrières écologiques. Selon leur revêtement, 
les routes et les trottoirs bétonnés peuvent cependant générer des conditions propices 
à la ponte en cas d’accumulation d'eau et de présence de plaques et de bordures 
dégouts (Montalvo et al., 2022) ; 

— les surfaces en eau, qui, selon leurs dimensions et leur contenant (eau claire, sale 
ou saumâtre) et l'écologie larvaire spécifique de l'espèce de moustique considérée, 
peuvent être de potentiels gîtes de ponte (p. ex. An. gambiae et flaques d’eau claire) ou 
des barrières écologiques (p. ex. cours d’eau, lacs pour Ae. albopictus); 

— les surfaces de forêt, et leur évolution au cours du temps, avec la prise en compte de 
la pression anthropique exercée sous la forme d'actions de déforestation, d'urbanisation 
et d'extension des terres agricoles, qui viennent perturber les écosystèmes à l'équilibre et 
favoriser les contacts humains-vecteurs (Orta-Pineda et al., 2021; Stefani et al., 2013); 
— les zones agricoles, qui sont susceptibles d'être exploitées comme gîtes larvaires, avec 
notamment l'importance des riziéres pour les anophèles (Diuk- Wasser et al., 2006), et qui 
constituent des zones d'exposition à la piqûre pour les personnes résidant à proximité; 
— le sol nu, qui peut jouer un rôle important en début et en fin de saison humide 
(formation de flaques d'eau stagnante), ainsi que les zones humides (wetlands), dans le 
cas de Anopheles notamment (Giardina et al., 2015); 

— la végétation, qui joue un rôle important en ville, en offrant un microclimat favo- 
rable aux moustiques vecteurs adaptés aux zones urbaines, sous la forme de gîtes de 
repos (zones d'ombre, de fraîcheur et d'humidité ou zone plus tempérée en hiver), et en 
favorisant l'existence de gîtes potentiels de ponte ou de sources d'alimentation. Selon 
la végétation du lieu, les moustiques peuvent se nourrir au travers de sucs et de nectar 
sur les arbres ou les fleurs (Honorio et al., 2009). La structuration, la hauteur, le type 
de végétation et la distance à une surface végétalisée sont des paramètres explicatifs 
importants de l'abondance des moustiques en milieu urbain, souvent considérés sous 
une forme agrégée (Manica et al., 2016). Diverses classes peuvent ainsi être consti- 
tuées, comme «arbres», «canne à sucre», «chaume», «pelouse», «sol à végétation 
éparse » (Machault et al., 2014); 

— la topographie de l’environnement, qui peut être associée au gradient de température 
(altitude, exposition), à une typologie de végétation ou à des considérations hydro- 
logiques, au travers d'une variable dérivée (pente) ou du calcul d'un indice hydrologique 
(p. ex. Topographic Wetness Index calculé par Homan et al. [2016] dans le contexte 
de l'étude du paludisme). La donnée de topographie peut être disponible, à moyenne 
résolution spatiale, sous la forme d'une couche globale (comme SRTM, tableau 1.5) ou 
générée à haute résolution par l'utilisateur (p. ex. modèle numérique de terrain dérivé 
d'une interférométrie radar, de données lidar ou d'un couple d'images stéréoscopiques). 


Différentes métriques paysagères disponibles dans la littérature permettent d'exploiter 
les cartes d'indices spectraux et les classifications d'occupation et d'usage du sol 
produites, en caractérisant l’hétérogénéité, la fragmentation, la composition (comme 
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le degré d’anthropisation, voir Orta-Pineda et al. [2021]) ou encore louverture d'un 
paysage (présence de clairières en forêt pour Anopheles, voir Stefani et al. [2013]). 
Linteraction humain-vecteur est parfois appréhendée simplement par la distance à un 
élément du paysage lié à l'écologie du vecteur, comme la distance à un dépôt de pneus 
ou à un cimetière dans le cas d’Ae. albopictus (Estallo et al., 2008), par le périmètre 
défini par la lisière de la forêt amazonienne pour An. darlingi (Li et al., 2016) ou bien au 
travers de statistiques géospatiales plus élaborées (Giardina et al., 2015). La présence 
d’habitats favorables à d’autres espèces susceptibles de servir d'hôtes intermédiaires 
doit également être prise en compte dans le cas d’un cycle zoonotique. 


En lien avec l'écologie du vecteur impliqué dans la transmission de la maladie, de 
nombreuses études ont ainsi exploité les données de télédétection issues des différents 
capteurs cités dans ce chapitre (optiques et, plus rarement, radar) dans des contextes 
biogéographiques décrits à différentes échelles et emprises (figure 1.3). Les études sur 
la dengue (Ae. albopictus et Ae. aegypti) et sur la fièvre à virus West Nile (Cx. pipiens, 
Cx. tarsalis) s’attachent surtout à décrire l'organisation spatiale des milieux urbains 
(figure 1.3a), tandis que les études sur le paludisme (An. gambiae notamment) 
incluent également les environnements agroforestiers aux interfaces des habitats 
humains (figure 1.3b). Les cartes de risque dérivées s'éloignent souvent de la réso- 
lution native des pixels pour être agrégées à des secteurs administratifs ou à d’autres 
types d’emprises incluant des sources de données complémentaires, qui caractérisent 
des aspects socio-économiques et/ou des pratiques des populations. Dans de rares 
cas, les données satellites sont désagrégées pour étre croisées avec des données plus 
finement résolues (Richman et al., 2018) ou sont intégrées dans des modèles de risque 
a différentes résolutions (Demets et al., 2020; Giardina et al., 2015). 
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Chapitre 2 


Indices spectraux et classifications 
d'images multispectrales pour 
la cartographie du risque vectoriel 


ANNELISE TRAN, RENAUD MARTI, VINCENT HERBRETEAU 


Les méthodes classiques d'analyse d'images de télédétection multispectrales! 
incluent la classification d'images et le calcul d’indices spectraux. Ces deux approches 
exploitent le lien existant entre les caractéristiques des surfaces observées et leur 
signature spectrale. 


Les méthodes de classification permettent de regrouper des zones de l’image ayant 
des réponses spectrales similaires dans la même classe. Deux types de classification 
peuvent être utilisés pour produire des cartes d'occupation du sol à partir d'une image 
multispectrale. Les méthodes de classification supervisée utilisent des algorithmes 
(maximum de vraisemblance, Random Forest, etc.) pour discriminer les différentes 
classes d'occupation du sol présentes dans l’image (eau, zones urbaines, zones de végé- 
tation, sols nus, etc.) et préalablement définies par l'utilisateur (Rodriguez-Galiano 
et al., 2012). Ces méthodes nécessitent donc de créer une base de données d’appren- 
tissage pour chacune des classes à cartographier avant de lancer la classification 
proprement dite. Les méthodes de classification non supervisée (p. ex. algorithmes 
K-means, Isodata) sont des processus plus automatiques qui ne nécessitent pas de 
connaissances a priori des classes (Duda et Canty, 2002). Les deux types de méthodes 
peuvent être appliqués pour classer chaque pixel de l’image (classifications orientées 
pixel) ou des objets (classifications orientées objet) en agrégeant des pixels voisins 
ayant des valeurs de réflectance similaires avant la classification, par un processus de 
segmentation de l’image (Blaschke et al., 2014). Dans les deux cas, des données de vali- 
dation doivent être collectées sur le terrain pour évaluer la qualité de la classification. 


Les indices spectraux exploitent la combinaison de différentes bandes d’une image 
multispectrale pour mettre en évidence certaines propriétés des surfaces observées 
(voir chapitre 1). Par exemple, la signature spectrale de la végétation révèle une forte 
réflectance dans le domaine du proche infrarouge et, au contraire, une faible réflec- 
tance dans le domaine du rouge (figure I.1). Pour mettre en évidence les zones en 
végétation d'une image multispectrale (figure 2.1a), les indices de végétation utilisent 


1. Voir Introduction générale. 
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la différence (comme le NDVI, Normalized Difference Vegetation Index, cf. chapitre 1) 
ou le rapport entre les valeurs de réflectance dans le proche infrarouge (PIR) et le 
rouge (R) [figure 2.1b] : 


NDVI = EE (équation 1) 
PIR+R 
Pour les zones en eau, dont la signature spectrale révèle des valeurs de réflectance 
décroissantes depuis les longueurs d'onde du domaine du bleu vers l’infrarouge, les 
indices d’eau (figure 2.1c) exploitent les différences entre les valeurs de réflectance 
dans le domaine du vert (V) et le proche ou moyen infrarouge (MIR), comme pour le 
MNDWI (Modified Normalized Difference Water Index) : 


V — MIR 
MNDWI = ——— (équation 2) 
V + MIR 
Le seuillage des indices spectraux permet ensuite l'identification et la cartographie des 
zones d'intérêt écologique pour les moustiques (figure 2.1d). 


N 
À 0 0,5 1km 
=) 
Figure 2.1. Exemples de calcul d'indices de végétation et d’eau a partir d'une image multi- 
spectrale. 


(a) Extrait d'une image Sentinel-2, composition colorée en couleurs naturelles. (b) Indice de végétation. 
(c) Indice d’eau — du noir, zones à faible valeur d'indice, au blanc, zones avec les valeurs d'indice les plus 
élevées. (d) Résultat du seuillage de l’image de l'indice d’eau — les surfaces en eau apparaissent en bleu. 
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Dans ce chapitre, différents exemples d'application de ces méthodes classiques d’inter- 
prétation d'images multispectrales sont présentés. Le premier est l’utilisation d'une 
carte d'occupation du sol dérivée d'une image multispectrale à haute résolution spatiale 
pour cartographier la distribution de moustiques Anopheles en Camargue, à une 
échelle locale. Dans le deuxième exemple, un indice d'eau est calculé à partir d'images 
à moyenne résolution spatiale et utilisé dans un modèle prédictif du risque d’occur- 
rence de cas humains de la fièvre du Nil occidental (West Nile fever) en Europe, à une 
échelle continentale. Enfin, nous présentons un exemple d'outil en ligne, Sen2Extract, 
qui permet d'extraire facilement des séries temporelles d'indices spectraux. 


» Cartographie de l'occupation du sol par télédétection pour 
modéliser La distribution des moustiques Anopheles en Camargue 


Contexte 


La Camargue est située dans le delta du Rhône dans le Sud de la France (figure 2.2a), avec 
un climat méditerranéen caractérisé par des étés chauds et secs, des hivers humides 
et doux et des précipitations importantes en automne et au printemps. Peu peuplée, 
cette région est une mosaïque de paysages de marais, sous l'influence des entrées d'eau 
salée par le sud, et de zones agricoles plus au nord, composées principalement de 
riziéres et prairies. 


Les différentes zones humides présentes en Camargue sont particulièrement propices 
au développement des moustiques, qui y sont très abondants, et à certaines maladies 
qu'ils transmettent. Le paludisme y était endémique jusqu'au début du xx° siècle (Jetten, 
Takken, 1994), puis il a été éliminé dans les années 1950 avec l'usage des insecticides 
(essentiellement le DDT) et la chimiothérapie/prophylaxie à base de quinine. Plus récem- 
ment, différents travaux de recherche ont étudié le risque de réémergence du paludisme 
dans cette région dans un contexte de réchauffement climatique (Ponçon et al., 2008). 
Ces études ont identifié l'espèce Anopheles (Anopheles) hyrcanus (Pallas), dont les rizières 
constituent le principal gîte larvaire, comme vecteur potentiel (Ponçon et al., 2007). 


Dans l'exemple présenté ici, les liens entre l'occupation du sol et des relevés entomo- 
logiques sur la présence de larves d'An. hyrcanus ont été analysés, afin de cartographier 
la distribution de cette espèce sur toute la Camargue (Tran et al., 2008). 


Données 


Données entomologiques 


Un total de 80 gîtes larvaires potentiels (différents types de marais, rizières) ont été 
visités mensuellement pendant la période d'activité des moustiques d’avril à octobre 
2006 par les équipes de l’Institut de recherche pour le développement (IRD). Les larves 
de moustiques ont été collectées selon les protocoles entomologiques standards et 
identifiées au niveau de l'espèce (Ponçon et al., 2007). 


Données de télédétection 


Deux images Landsat 7 Enhanced Thematic Mapper (ETM+, résolution spatiale de 
30m) ont été acquises 4 deux périodes contrastées (en juillet et en octobre) pour 
cartographier les principaux types d’occupation du sol dans la zone (figure 2.2a). 
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Méthodes 


Classification de l'occupation du sol 


Une classification supervisée orientée objet (logiciel eCognition?) a été appliquée, à 
partir d'un jeu de données d'entraînement construit sur la base d'observations terrain et 
en concertation avec des entomologistes, afin de cartographier les habitats de plusieurs 
espèces retenues pour leur intérêt, pour cinq principaux écotypes (figure 2.2b). Cette 
typologie inclut l'ensemble des gîtes larvaires potentiels pour les moustiques Anopheles 
dans la région : eau libre, rizières, roselières, marais à scirpes, marais à joncs. 


Analyses statistiques 


Les corrélations entre la présence ou l’absence de larves d'An. hyrcanus (variable à 
expliquer) et le type de gîte larvaire et la distance aux rizières (variables explicatives) 
ont été étudiées à l’aide d'un modèle de régression logistique. 


Principaux résultats 


Les rizières, les marais à scirpes et les roselières ont été identifiés comme significa- 
tivement corrélés à la présence de larves d'An. hyrcanus par le modèle de régression 
logistique. Ces résultats ont permis de confirmer le potentiel d’autres biotopes que les 
rizières comme gîtes larvaires d'An. hyrcanus, ce qui explique la présence de l'espèce 
en fin d'été et en automne lorsque les riziéres sont à sec. D’autre part, une carte de 
probabilité de présence de l'espèce a pu être générée sur toute la zone d'étude à partir 
de la carte d'occupation du sol (figure 2.2c) et cette carte a pu être validée par des 
données de terrain sur la distribution des moustiques An. hyrcanus adultes. Il faut 
noter que la même approche appliquée à des données d’autres espèces de moustiques 
Anopheles et Culex a montré que la distribution des gîtes larvaires en zone humide 
impactait la répartition spatiale des espèces à l'état adulte, lorsque celles-ci étaient 
fortement inféodées à un type de gîte larvaire (Cailly et al., 2011). 


415E WE 445E 


415E #30E A'4SE 


Indice larvaire pour An. hyrcanus E Zones urbaines 
OB: zz 


Zones en eau 


Figure 2.2. Cartographie d’un indice des zones favorables au moustique Anopheles hyrcanus en 
Camargue, à partir d'une image satellite optique. Adapté de Tran et al., 2008. 


(a) Extrait d'une image Landsat 7 ETM+. (b) Résultat de la classification en différents types d'occupation du 
sol. (c) Indice larvaire pour l'espèce An. hyrcanus, variant de 0 (zone très peu favorable) à 1 (très favorable). 


2. https://fr.geospatial.trimble.com/products-and-solutions/trimble-ecognition 
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» Indices spectraux dérivés d'images de télédétection, 
comme facteurs environnementaux des cas humains de fièvre 
du Nil occidental en Europe 


Contexte 


Le virus du Nil occidental (West Nile virus, WNV), responsable de la fièvre du Nil 
occidental (West Nile fever, WNF), est transmis par piqûre de moustiques du genre 
Culex, avec comme hôtes principaux des espèces d'oiseaux domestiques et sauvages, 
et parfois une transmission à l'humain ou au cheval, considérés comme des hôtes acci- 
dentels et des impasses épidémiologiques (dead-end host). En Europe, la circulation du 
WNYV a été confirmée épisodiquement ces dernières décennies dans plusieurs pays du 
bassin méditerranéen, mais durant l'été 2010, particulièrement chaud, un très grand 
nombre de cas humains de WNF a été rapporté dans des zones jusqu'alors indemnes 
(Paz et al., 2013). L'étude résumée ici visait à identifier les principaux déterminants 
environnementaux du risque d’occurrence de cas humains de WNF en Europe, en 
testant en particulier des indices spectraux issus de la télédétection (Tran ef al., 2014). 


Données 

Données épidémiologiques 

Les données sur les cas humains confirmés de WNF ont été répertoriées à l'échelle 
du district en Europe sur la période 2002-2013. Pour chaque année, chaque district 
(n = 1 113) était considéré comme «infecté » si au moins un cas humain de WNF avait 
été rapporté dans ce district et «non infecté » dans le cas contraire. 

Données de télédétection 

Une série temporelle d'images MODIS a été téléchargée afin de couvrir toute l’Europe, 
et ce, avec une fréquence de 8 jours pour toute la période considérée. 


Autres données spatialisées 


D'autres sources de données à référence spatiale ont été utilisées pour caractériser, à 
l'échelle du district, des facteurs de risque potentiels d’occurrence de cas humains de 
WNF : températures, population, routes de migration aviaire et zones humides. 


Méthodes 


Calculs d'indices de végétation et d'eau 

Les indices de végétation (NDVI) et d'eau (MNDWI) présentés précédemment (équa- 
tions 1 et 2) ont été calculés à partir de la série temporelle d'images MODIS. 
Analyses statistiques 


Les corrélations entre le statut d'infection par le WNV (variable à expliquer) et la popu- 
lation, la présence de zones humides, la présence de routes de migration aviaire, les 
valeurs d'anomalies de température, de NDVI et de MNDWI (variables explicatives) 
ont été étudiées par une analyse univariée (en testant les variables une à une), puis par 
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une analyse multivariée, combinant les variables les plus significatives à l'aide d'un 
modèle de régression logistique. Les données épidémiologiques de la période 2001- 
2011 ont été utilisées pour construire le modèle et les données épidémiologiques de la 
période 2012-2013, pour tester sa capacité prédictive. 


Résultats 


Le meilleur modèle permettant d'expliquer le statut infecté/indemne de cas humains 
de WNF par district inclut comme variables explicatives les anomalies de température 
au mois de juillet et d'indice d'eau MNDW/I début juin (9-16 juin), la circulation du 
WNV l’année précédente, la présence de zones humides, le type de route migratoire 
aviaire et la population humaine (figure 2.3a) ; toutes ces variables étant positivement 
et significativement corrélées avec la probabilité d'infection (Tran et al., 2014). A 
partir de ce modèle de régression logistique, des cartes de la probabilité d'infection par 
le WNV par district peuvent être produites annuellement (figure 2.3b). Ces résultats 
mettent en évidence le rôle des zones en eau dans le risque de transmission du WNV 
en Europe et la pertinence d’un indice d’eau simple (le MNDWI dérivé des images 
MODIS) pour détecter des surfaces en eau au-dessus des moyennes saisonnières au 
mois de juin, susceptibles de favoriser la prolifération des moustiques durant les mois 
d'été et la transmission du virus à l'humain en fin d'été-début d'automne. 
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Figure 2.3. Cartographie des zones à risque de transmission du virus du Nil occidental en Europe 
basée sur des indicateurs environnementaux et météorologiques. Adapté de Tran et al. (2014). 


(a) Analyse d'un indice d’eau et de différents prédicteurs environnementaux du virus du Nil occidental par 
un modèle de régression logistique. (b) Cartographie du risque d'infection. 
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» Production automatisée d'indices spectraux : 
exemple de l'outil Sen2Extract 


La télédétection reste un domaine technique demandant une bonne connais- 
sance des méthodes, des logiciels et algorithmes de traitement d'images et du large 
éventail des données spatiales disponibles. Il n’est donc pas étonnant que l'essentiel de 
l'usage des images satellite dans le domaine de la santé se soit fait indirectement via la 
mise à disposition de produits pré-calculés, tels que les indices spectraux, gratuitement 
et simplement accessibles par téléchargement. Ces produits ont permis d’atteindre 
une communauté d'utilisateurs, essentiellement des épidémiologistes et des statisti- 
ciens, sans compétences en télédétection, qui ont pu les intégrer dans leurs études. 
Bien que ces indices ne soient pas forcément les mieux adaptés aux caractéristiques 
des zones d'étude ou à la problématique, dans bien des cas leur disponibilité constante 
dans l’espace et le temps en a fait de précieux indicateurs des fluctuations saisonnières. 


C'est à partir de ce constat que des outils ont été développés pour automatiser la 
production d'indices spectraux à partir des images des satellites Sentinel-2A et 2B 
du programme Copernicus de l’Agence spatiale européenne (ESA). Comme leur 
nom le suggère, ces deux satellites ont été conçus à des fins de surveillance envi- 
ronnementale et ils peuvent donc être indirectement utiles pour la surveillance de 
maladies liées à l’environnement. En effet, ils fournissent des images optiques ayant 
une résolution spatiale de 10m, disponibles gratuitement tous les 5 jours, depuis 
mi-2015, sur la quasi-totalité du globe. Parmi les outils développés, Sen2Extract est 
une application en ligne? qui permet d'extraire simplement des séries temporelles 
d'indices spectraux sur une période et des sites choisis selon les besoins de l’utili- 
sateur. Une interface simple et intuitive (figure 2.4) invite à indiquer son/ses sites 
d'intérêt en chargeant un fichier au format vectoriel (points ou polygones, format 
shapefile zippé), puis à choisir la période et les indices à calculer (indice de végé- 
tation NDVI et différents indices d’eau). Dès que le calcul est réalisé, un e-mail est 
envoyé pour indiquer le lien de téléchargement du résultat de la requête, qui prend 
la forme d’un tableau de données au format csv. Pour chaque entité spatiale, poly- 
gone ou point, et chaque date, des statistiques descriptives des valeurs de l'indice 
choisi sont calculées (moyenne, minimum, maximum, écart-type, médiane, premier 
et troisième quartiles, nombre de pixels ennuagés ou sans données). Le fichier final 
est donc très léger (quelques kilo-octets), comparé au poids d’une image Sentinel-2 
et de ses traitements qui peuvent représenter près d'un gigaoctet. Ceci facilite l’uti- 
lisation des images dans des pays où l'accès à internet est limité. 


Sen2Extract est l'interface utilisateur d'un ensemble d’algorithmes développés dans le 
cadre du projet S2-Malaria (appel à projets de recherche Tosca, financé par le Cnes, 
de 2017 à 2020) pour produire des données d'observation satellitaire utiles à la surveil- 
lance de maladies. Ces algorithmes ont ensuite été optimisés pour le suivi des dégâts liés 
aux passages des cyclones (Projet Interreg Renovrisk-impact de 2018 à 2020). En amont 
de Sen2Extract, la chaîne Sen2Chain, codée en Python, est au cœur de la production 
des indices. Elle permet de télécharger les images Sentinel-2 sur deux des principaux 
catalogues, le Copernicus Open Access Hub* de PESA et la plateforme d'exploitation 


3. https://web.seas-oi.org/sen2extract/ 
4. https://scihub.copernicus.eu/ 
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des produits Sentinel’ (Peps) du Cnes. Elle intègre l'algorithme Sen2Cor de ESA pour 
effectuer les corrections atmosphériques, de terrain et des cirrus sur les images brutes 
(niveau 1A) afin d'obtenir le niveau 2A, nécessaire à la comparaison des dates et des 
scénes. Elle permet ensuite de produire différents indices spectraux — il est aussi possible 
den intégrer de nouveaux selon les besoins —, puis d’extraire les valeurs des pixels pour 
calculer leurs statistiques descriptives sur des entités spatiales définies, à chaque date. Les 
codes de Sen2Chain sont ouverts® et peuvent donc être réutilisés et améliorés librement. 
Sen2Extract se greffe en aval pour offrir l'interface de requête et lancer le calcul des séries 
temporelles d'indices selon les besoins. Cet outil est codé en R” et l'interface en ligne 
repose sur le package R Shiny®. Les codes sont également ouverts pour pouvoir être réim- 
plémentés’. Ils permettent d'envoyer les requêtes vers la chaîne Sen2Chain, de paramétrer 
les extractions et les calculs statistiques puis d'envoyer le lien par e-mail. L'interface acces- 
sible en ligne est installée sur les serveurs de la station de réception satellitaire SEAS-OI à 
La Réunion. Les requêtes de séries temporelles peuvent se faire sur l'emprise des images 
qui ont été acquises sur la région océan Indien essentiellement. Pour d’autres zones, il 
convient donc d'implémenter ces outils sur un nouvel ordinateur ou serveur. 


L'ensemble des tâches proposées par ces chaînes est entièrement automatisable. Ceci 
permet de produire des séries temporelles en routine pour le besoin d'outils de surveil- 
lance. C’est le cas de l'application Lepto Yangon”? qui peut réaliser des cartographies des 
environnements favorables à la transmission de la leptospirose sur l’agglomération de 
Yangon au Myanmar, dès qu’une nouvelle image est disponible (tous les 5 jours). Une 
telle automatisation de la production d'indices spectraux dans des temps courts serait 
envisageable pour le suivi d’environnements favorables aux espèces de moustiques pour 
lesquelles les liens entre ces indices et leur présence ou abondance ont pu être établis. 


= Dä s.,2c 
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Browse... No file selected Tiles in database 


Chagos Archipetaga 


Choose a spectral index (multiple 
choice possible) : 

NDVI 

C NDWIGAO 

C NDWIMCF 

MNDWI 

C NDRE 

(1 IRECI 


Temporal smoothing 


O Activate 2 


Select the date range : SES Leaflet | © OpenStreetMap contributors 


2022-06-01 | to | 2022-06-30 


Figure 2.4. Interface de l'outil en ligne Sen2Extract (https://web.seas-oi.org/sen2extract/). 


5. https://peps.cnes.fr 

6. https://framagit.org/espace-dev/sen2chain 
7. https://www.R-project.org/ 

8. https://shiny.rstudio.com/ 

9. https://framagit.org/espace-dev/sen2extract 
10. https://leptoyangon.geohealthresearch.org/ 
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Chapitre 3 


Estimation des températures de Pair 
à partir d'images satellite 
et de stations météorologiques 


BARBARA BOUFHAL, ALEXANDRE CEBEILLAC, ÉRIC DAUDÉ 


Les moustiques du genre Aedes sont thermosensibles, c'est-à-dire que leur niveau d’acti- 
vité (recherche de repas sanguin, déplacements, etc.) et la vitesse de développement de 
leurs larves sont conditionnés par les températures. La zone de confort de ces mous- 
tiques se situe entre 21°C et 32°C, avec une probabilité de survie qui se réduit au-delà 
de ces bornes (Brady et al., 2013). A l'échelle du globe, si l'incidence des maladies dues 
à ces moustiques est donc pour l'essentiel concentrée dans les zones tropicales et sub- 
tropicales (Bhatt et al., 2013), la présence des vecteurs s'étend aujourd’hui jusqu'aux 
zones tempérées d'Amérique du Nord et d'Europe du Sud notamment, à la faveur du 
changement climatique (Messina et al., 2019). Les moustiques de l'espèce Aedes aegypti 
sont également synanthropiques, c'est-à-dire que leur cycle de vie est étroitement lié à la 
présence de populations humaines. Les femelles Ae, aegypti privilégient en effet de petits 
gîtes créés par l’homme (pots de fleurs, citernes à eau, détritus, etc.) pour pondre et ont 
besoin de repas de sang pour la maturation de leurs œufs. En concentrant la chaleur et 
les humains, les villes sont donc de parfaits attracteurs pour ces populations de vecteurs. 


Dans ce chapitre, nous nous intéressons plus particulièrement à l'estimation des tempé- 
ratures en milieu urbain en ce qu'elles impactent potentiellement les dynamiques de 
population vectorielle. En effet, les différentiels de température présents à micro-échelle 
pourraient, selon les contextes, favoriser ou limiter le développement biologique des 
moustiques et donc leur propension à se déplacer, à prélever des repas de sang, mais 
aussi agir sur la durée de leur cycle de reproduction et leur taux de survie. Connaître 
précisément les températures à échelle fine est donc un préambule à la détection d’envi- 
ronnements locaux qui leur sont favorables ou défavorables. Nous présentons dans ce 
chapitre différents types de données pour caractériser les températures du sol et de l'air, 
puis un ensemble de méthodes pour estimer les températures de l'air et, pour conclure, 
une comparaison des résultats de ces méthodes appliquées à Bangkok (Thaïlande). 


» Des données pour mesurer les températures 


La connaissance des températures, ou tout au moins de leurs variations relatives à échelle 
fine, doit permettre de caractériser plus précisément les micro-environnements qui sont 
susceptibles de favoriser ou non le maintien et le développement de populations de 
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moustiques. Les îlots de chaleur urbains sont caractéristiques de ce différentiel de tempé- 
rature qui peut être mesuré en zone intra-urbaine. La présence d’un îlot de chaleur est due 
à la concentration de zones bâties qui restituent la nuit la chaleur emmagasinée le jour, 
contribuant ainsi à l'augmentation des températures minimales des zones densément 
peuplées. Ceci peut donc favoriser d'autant plus le maintien d’Aedes, notamment lors des 
périodes les plus froides pour les zones urbaines des régions intertropicales — comme à 
Bangkok où la dengue est endémique — et réduire sa probabilité de survie lors des périodes 
les plus chaudes. Classiquement, deux classes de températures sont mobilisables pour 
identifier la présence d’ilots de chaleur : les températures au sol et les températures de l'air. 


La mesure des températures au sol (LST pour Land Surface Temperature) passe par 
l'imagerie satellitaire via différents capteurs. Les images sont proposées à différentes 
résolutions spatiales et temporelles. Nous pouvons citer le capteur TIRS embarqué 
sur Landsat 8, avec lequel il est théoriquement possible d'estimer les températures de 
surface à une résolution spatiale de 100m pour une périodicité bimensuelle. Pour une 
périodicité plus courte et donc plus adaptée au cycle de développement du moustique, 
le capteur MODIS embarqué sur les satellites Aqua et Terra fournit des tempéra- 
tures du sol à une résolution de 1km, mais quatre fois par jour, lorsque les conditions 
météorologiques sont optimales. Lune des contraintes fortes en effet pour les images 
satellite, notamment dans les régions tropicales et subtropicales, est la présence d'un 
plafond nuageux très fréquent qui les rend peu exploitables. 


Si les températures varient spatialement à la surface du sol, elles varient également 
lorsqu'on les compare à celles de l'air (figure 3.1). Ces différences s'expliquent par la 
combinaison de nombreux facteurs qui agissent localement : humidité du sol, rugosité 
de surface, vitesse du vent (Stisen et al., 2007), etc., et ces différences de température 
ne sont pas sans effet sur le développement du vecteur (Misslin et al., 2016). 


La mesure des températures de l’air passe par des stations météorologiques. Celles- 
ci sont positionnées selon des normes très précises, à l'ombre notamment, et gérées 
par des organismes officiels locaux. Les relevés sont donnés à 2m d'altitude et les 
contrastes de valeur enregistrés entre les différentes stations d’une même ville rendent 
compte des variations locales du bilan énergétique de surface. Lun des défis dans 
l'utilisation de ces données est généralement la densité du maillage territorial, le 
nombre de stations étant souvent trop faible pour réaliser des interpolations fiables à 
l'échelle de grandes métropoles architecturalement et structurellement hétérogènes. 


Ce chapitre propose donc d'associer ces deux sources de données : des images satellite 
(LST) qui ne sont pas complètes dans le temps et des relevés de stations météoro- 
logiques qui ne sont pas complets dans l’espace, pour comparer différentes méthodes 
d'estimation des températures de lair dans la métropole de Bangkok (figure 3.1). 
Pour cette étude, nous avons fait le choix de privilégier le capteur MODIS pour les 
fréquences élevées de passage et la simplicité d'accès aux images (site de l'USGS!). 


» Estimation des températures de l'air : différentes méthodes 


Il existe différentes méthodes pour estimer les températures de l’air à partir de la 
connaissance des températures de surface. Une première méthode simple repose sur 
des régressions linéaires entre les températures de l'air enregistrées par les stations au 


1. https://Ipdaac.usgs.gov/data/get-started-data/collection-overview/missions/modis-overview/ 
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sol et les températures du sol obtenues à partir d'images satellite. Si le nombre d’enre- 
gistrements est important (réalisés sur une longue période et avec un nombre suffisant 
de stations) et le coefficient de corrélation significatif, la régression permet d'estimer 
les températures de l'air sur l'ensemble de la zone à partir des données de tempé- 
rature du sol (Zaksek et Schroedter-Homscheidt, 2009). Un autre type d'approche 
utilise le bilan énergétique des surfaces observées et mobilise des concepts de ther- 
modynamique (Sun et al., 2005). Enfin, des régressions (linéaires ou non) entre la 
température de surface et d’autres variables annexes, comme l'occupation du sol ou 
l'altitude, peuvent être utilisées (Benali et al., 2012; Oyler et al., 2014). La méthode 
TVX (Temperature- Vegetation Index) repose sur ce principe. 


Une première hypothèse de la méthode TVX postule un équilibre thermique entre les 
températures du sol et de lair au niveau de la canopée (Vancutsem et al., 2010). Elle 
part également du postulat qu'il existe une corrélation négative entre la température 
et l’activité végétative, mesurée par l'indice de végétation par différence normalisée 
(NDVI pour Normalized Differenced Vegetation Index), et que les conditions atmos- 
phériques sont uniformes localement (Stisen et al., 2007). A partir de ces principes, la 
méthode TVX consiste à calculer une régression linéaire entre les valeurs de LST (Ts) 
et celles de NDVI dans une fenêtre mobile de taille donnée : 


Ts~aNDVI + B (équation 1) 


En se basant sur l'hypothèse de l'équilibre thermique de la canopée, on conserve les 
paramètres de l'équation 1 et on pose alors que : 


Ta-aNDVImax + B (équation 2) 


ou NDVImax peut être soit choisi parmi les valeurs de NDVI des stations au sol de sorte 
qu'il minimise l'erreur moyenne absolue de la relation de l'équation 2 (méthode 1 ou 
« TVX classique »), soit le NDVI maximum dans la fenêtre locale (méthode 2 ou « TVX 
local»). L'avantage de la méthode TVX est qu'elle nécessite peu de données auxiliaires 
en entrée et que, contrairement à l’approche statistique, elle ne requiert pas un réseau 
de stations trés dense pour donner des résultats corrects (Misslin et al., 2018). 


Nous appliquons par la suite la méthode TVX à la métropole de Bangkok (méthode 1 
«TVX classique») en proposant différentes variantes, telles que l’utilisation du 
NDVImax local (méthode 2 « TVX local») et la correction des températures en lien 
avec la part de surface bâtie (méthode 3 « TVX local & bâti»). En effet, si les bâtiments 
influencent les températures en générant l'effet îlot de chaleur, leur prise en compte 
peut potentiellement améliorer l'estimation des températures de l'air. 


» Applications à Bangkok 


Pour ce cas d'étude, nous nous appuyons sur la collecte des images de LST et de NDVI 
issues du capteur MODIS de Terra et Aqua pour la période du 30 juin 2015 au 30 avril 
2022. Au total, 1 885 images sont potentiellement exploitables, avec une couverture d’au 
moins 50 % de la zone. Néanmoins, les mois de mai à novembre, du fait d'une couverture 
nuageuse importante, ne représentent que 17 % de l'échantillon (tableau 3.1). 


Nous utilisons également une couche géographique contenant plus de 1,4 million de 
bâtiments à Bangkok (voir chapitre 4 pour le détail de la méthode). Nous appliquons 
pour chacune des images une régression linéaire dans une fenêtre mobile (entre la 
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température du sol et le NDVI pour les méthodes 1 et 2 et avec l'ajout de la couche 
de batiment pour la méthode 3). Nous obtenons alors les coefficients de régression 
locaux, ainsi que les ordonnées à l’origine, comme explicité dans la figure 3.2. 


MODIS LST MODIS NDVI Batiments 


i 


Fenêtre mobile 7x7 


NDVImax Modèle linéaire : LST ~ NDVI + Bâtiments 


onto de régression 


_ Calibration : 
Ta ~ intercep + Coef1*VALEUR* + Coef2*Batiments Meilleure valeur NDVI 
Max global 


Températures des stations 


Validation 


Figure 3.2. Présentation des méthodes employées. 


La méthode 2 (« TVX local») suit le cheminement en noir, auquel on rajoute les éléments en rouge pour la 
méthode 3 (« TVX local & bâti») ou ceux en bleu pour la méthode 1 (approche « TVX classique »). "Valeur" 
minimise la moyenne des erreurs pour l’ensemble des stations pour la méthode 1, tandis qu’elle correspond 
au NDVImax local pour les méthodes 2 et 3. 


Concernant les données météorologiques, sept stations sont présentes sur la zone 
d'étude, mais nous ne disposons d’informations sur les surfaces bâties que pour cinq 
dentre elles (les plus proches des limites administratives de Bangkok, figure 3.1). 
Ces dernières fournissent les températures de l’air avec une fréquence d'une à trois 
heures, mais seules les températures aux heures des passages les plus proches des 
satellites Aqua et Terra sont conservées a des fins de calibration et de validation. Nous 
proposons par la suite trois approches issues de la méthode TVX afin de les comparer. 


La première méthode (méthode 1) reprend l’approche de Misslin et al. (2018). Elle 
consiste à trouver la valeur NDVImax telle que l'écart entre les températures de lair 
enregistrées par les stations et les températures de l'air estimées par l'équation 2 soit 
minimisé. La méthode consiste donc à utiliser l'ensemble des valeurs de NDVI enre- 
gistrées à la position de chacune des stations et de sélectionner celle qui donne l'erreur 
moyenne absolue (MAE pour Mean Absolute Error) la plus faible : 


Ta'-T 
mag - ZF#"-Ta 
n 


avec Ta’ la température de l'air estimée, Ta la température de l'air enregistrée et n le 
nombre d'observations. 


C'est cette valeur de NDVImax, soit une valeur pour l’ensemble de la zone, qui est 
alors utilisée dans l'équation 2 pour estimer les températures de l'air. La mise à jour 
saisonnière de NDVImax permet d'estimer les évolutions des températures de l'air sur 
la zone (tableau 3.1). 
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Tableau 3.1. Valeurs du NDVImax par période (par mois). 


Janvier Février Mars Avril Mala Décembre 
novembre 

NOVa 0,56 0,49 0,47 0,45 0,62 0,58 
(diurne) 

NON 0,08 0,02 0,08 0,06 0 0,05 
(nocturne) 

e 333 388 334 200 328 302 
d'images 


Remarquons que si ces valeurs NDVImax varient entre les périodes, elles varient 
également à différents moments de la journée (tableau 3.1). Les très faibles valeurs 
de NDVImax obtenues pour la nuit impliquent l'importance de lordonnée à l'origine 
(le coefficient DI dans l'estimation des températures nocturnes. En effet, si la valeur 
de NDVImax est proche de zéro, Ta’ devient alors très proche de f (équation 2). Ce 
premier constat devra être pris en compte lors des comparaisons des qualités des 
méthodes, notamment dans l'estimation des températures nocturnes. 


Les deux autres méthodes testées (méthodes 2 et 3) calculent un NDVImax local à 
partir d'une fenêtre mobile. Par rapport au NDVImax global calculé par la méthode 1, 
l'hypothèse est que le bilan énergétique varie fortement en zone urbaine du fait de la 
variété des usages du sol, ce qui justifie de prendre en compte localement cette source 
de variation. La méthode 2 fait ainsi un focus sur le lien entre les températures de 
surface et le NDVImax local, et la méthode 3 (TVX & bâti) prend en compte la part 
de zones bâties dans le calcul des coefficients de régression, en modifiant l'équation 1 : 


Ts~a NDVI + a båti + B (équation 1 modifiée) 


La figure 3.3 montre la distribution des coefficients de corrélation calculés pour chaque 
image MODIS dans une fenêtre mobile (figure 3.3a) et les moyennes locales obtenues 
par la méthode 2 (figure 3.3b) et 3 (figure 3.30). 


Nombre d’images 
a 2004 


Ts ~ NDVI + bâti 


Méthode 
E Ts~ NDVI E Ts ~ NDVI + bâti 


Figure 3.3. Comparaison des coefficients de corrélation entre les températures de surface et le 
NDVI pour la méthode 2, et les températures de surface et le NDVI avec la prise en compte de 
la densité de batiments pour la méthode 3. 


(a) Histogramme des valeurs du coefficient de corrélation R? global calculé pour chaque image MODIS. (b) 
R? local (fenêtre 7 x 7) moyen sur l'ensemble des images pour la relation Ts ~ NDVI. (c) Ts ~ NDVI + bâti. 
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Nous pouvons noter que l'ajout de la couche de bâti dans l'application de régres- 
sions linéaires locales permet d'augmenter les R? globaux (moyenne à 0,18 pour la 
méthode 2, et 0,25 pour la méthode 3, figure 3.3a). Le même constat peut être fait 
localement (figure 3.3b et c), mais si un grand nombre de zones présentent des coef- 
ficients de corrélation élevés, il est délicat de conclure à une quelconque influence de 
la structure urbaine sur la relation Ts - NDVI + bâti. Si le niveau de corrélation n’a 
pas nécessairement besoin d’être élevé pour appliquer la suite de la méthode TVX 
(Misslin et al., 2018), on peut observer malgré tout que l’utilisation d'une couche 
de bâti en plus du NDVI augmente fortement la relation et pourrait améliorer 
l'estimation des températures de l'air. 


La figure 3.4 présente les résultats obtenus pour chaque méthode, avec, à titre de 
comparaison, les LST des images MODIS. La figure 3.4a montre qu'il y a d’une part une 
surestimation des températures de l'air évaluées le jour et une sous-estimation la nuit. 


a Densité Jour Nuit 
0,15 H i 


0,10 


0,05 


b Valeur MAE R2 
4 


Source [_]LST [_] Méthode? ` ` | Méthode2 [_] Méthode 3 


Figure 3.4. Ecart des températures de lair estimées par rapport aux températures observées. 


Pour chaque méthode : (a) le jour (gauche) et la nuit (droite); (b) : erreur absolue moyenne (gauche) et 
coefficient de corrélation (droite) sur 6 périodes de l’année. 
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Cependant, le jour, les modes des températures estimées sont plus proches de ceux des 
températures enregistrées que le mode des LST. À noter que la méthode 1, bien que 
présentant une amplitude des écarts plus importante, permet d’avoir plus de valeurs 
proches des températures des stations. En revanche, la nuit, les méthodes d'estimation 
des températures de l'air ne réduisent pas les écarts avec les températures observées, 
surtout pour la méthode 1. La figure 3.4b montre l'écart des estimations sur différentes 
périodes de l’année, sous le prisme du MAE et du R°. Au regard de ce dernier indicateur, 
les températures de l'air pour les mois de mars et d’avril sont les plus fiables, notam- 
ment avec la méthode 2 qui présente des chiffres légèrement meilleurs que ceux de la 
méthode 3 (TVX + bâti). La méthode 1 réduit effectivement l'erreur absolue moyenne 
par rapport aux LST, mais reste moins efficace que les méthodes 2 et 3, où les écarts 
oscillent entre 3,2 et 3,7°C selon les mois. 


La méthode de validation, qui passe par la comparaison avec les stations météo- 
rologiques au sol, souffre de plusieurs écueils. Tout d’abord, comment comparer 
efficacement des données à une échelle méso à des données à l'échelle micro? Les 
valeurs enregistrées par ces dernières dépendent de leur environnement proche, alors 
que les estimations des températures de l'air sont des moyennes sur des zones nette- 
ment plus importantes. De plus, le nombre de stations est relativement faible (5) et 
leur maillage n'est pas suffisamment homogène pour représenter tous les assemblages 
des structures urbaines de Bangkok, susceptibles d’influencer les températures (zone 
d'habitat dense, éparse, hybride, etc 1 


Si nous comparons les écarts moyens entre les différentes méthodes d'estimation des 
températures de l'air et les LST (figure 3.5), nous pouvons noter la présence de zones 
où les températures de l'air sont supérieures aux températures du sol issues de MODIS, 
notamment dans le centre-ville. L'inverse s'observe également, surtout à l’est de la ville, 


a) Ensemble de Bangkok 
TVX(1) vs LST TVX(2) vs LST TVX + bati(3) vs LST 
vi SE ir : 


EEH im 


b) Zoom au sud de Bangkok 


TVX(1) vs LST TVX(2) vs LST TVX + bâti(3) vs LST 
L_ = ii] ("e 6 km Le "Eu 
Différence de température (°C) We ` O 


0 


Figure 3.5. Différences entre la somme des températures diurnes estimées (de gauche à droite : 
méthodes 1, 2 et 3) et les LST issues de MODIS, pondérées par le nombre de pixels contenant 
une information sur l’ensemble de la période : (a) sur l’ensemble de Bangkok, (b) sur une zone 
plus petite et hétérogène au sud de la ville. 
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dans la zone de l'aéroport (station 48429) [figure 3.1]. Une grande hétérogénéité dans 
l'occupation du sol à une échelle locale peut expliquer ces écarts, notamment dans la 
zone du méandre au sud de Bangkok (figure 3.5b), où des zones essentiellement végé- 
talisées côtoient des zones densément peuplées. Cette hétérogénéité locale influence 
les coefficients de régression locaux et le NDVImax local. 


» Conclusion 


Nous avons comparé ici différentes méthodes permettant d'estimer les températures 
de l'air à partir des températures du sol provenant d'images MODIS. Si les traitements 
proposés permettent de réduire les écarts entre les températures de l'air observées et 
estimées, ces derniers restent relativement importants et varient de 3 à 5°C selon les 
périodes de l’année et les méthodes. Il est toutefois délicat de comparer des tempé- 
ratures estimées à l'échelle méso à des données mesurées à micro-échelle. En effet, il 
est tout à fait possible que la température moyenne de l'air effective sur une maille de 
1km? soit différente de la température enregistrée localement, du fait de l’emplace- 
ment de la station et de son environnement local d'une part, et d'effets de lissage sur 
la maille d'autre part. La qualité de l'estimation dépend aussi de la densité de mail- 
lage des stations météorologiques de référence et de leur environnement local, pas 
nécessairement homogène sur la maille d'analyse. 


Notons également que les méthodes employées entraînent des effets de bordure liés à 
l'utilisation d'une fenêtre mobile dans des zones hétérogènes. Le calcul des coefficients 
de régression locaux sur des zones homogènes d'un point de vue de l'occupation du sol 
permettrait d'affiner la comparaison des différentes méthodes. 
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Chapitre 4 


Du recensement au bâtiment : 
génération de populations synthétiques 


ALEXANDRE CEBEILLAC, OLIVIER GILLET, ÉRIC DAUDÉ 


La connaissance de la répartition des populations humaines à haute résolution spatiale 
est un enjeu majeur pour l'étude des maladies vectorielles. Dans le cas de la dengue, la 
taille des foyers infectieux varie en effet entre 200 et 800 m de rayon (Salje et al., 2012), 
ce qui induit de possibles explosions épidémiques dans les zones les plus peuplées. Ceci 
conduit à une stratégie éprouvée de la lutte antivectorielle avec des interventions ciblées 
autour des domiciles de cas index pour limiter l'apparition et la taille des clusters. En 
cas de présence de plusieurs foyers, cette stratégie doit être affinée pour prioriser les 
interventions dans les zones les plus à risque, donc potentiellement les plus peuplées, 
avec une présence avérée de vecteurs. Le couplage de données épidémiologiques 
(p. ex. via les hôpitaux sentinelles) et entomologiques (comme des cartes de risque de 
présence de moustiques) avec la répartition de la population doit donc permettre de 
mieux anticiper les risques d'apparition et de propagation de foyers infectieux. 


Véritables outils d'intervention sur le terrain, les cartes de population ou de densité 
de population à une échelle spatiale très fine ne sont malheureusement pas dispo- 
nibles dans une grande partie des métropoles des zones intertropicales où la dengue 
sévit tout particulièrement. En effet, les données institutionnelles sont généralement 
agrégées selon un découpage administratif qui ne fait pas sens pour les mécanismes 
épidémiologiques. Par ailleurs, elles sont souvent trop anciennes au regard des taux de 
croissance élevés que connaissent certaines mégapoles et, parfois aussi, inexistantes. 
Si ces deux derniers aspects, liés aux capacités de recensement de chaque État et/ou 
municipalité, ne sont pas du ressort de ce chapitre, il est toutefois possible d'appliquer 
des algorithmes de désagrégation spatiale des données disponibles à une échelle large 
afin de ventiler la population à des résolutions spatiales plus fines (Viel et Tran, 2009). 
Des projets comme Landscan (Dobson et al., 2000), Worldpop (Linard et Tatem, 2012) 
ou Global Human Settlement Layer (GHSL; Florczyk et al., 2019) vont dans cette 
direction en proposant des cartes de répartition de la population à l'échelle mondiale 
avec des résolutions spatiales maximales de l’ordre du kilomètre et jusqu’à 250 m. 


Pour aller plus loin dans la désagrégation spatiale des données de population, nous 
présentons un cas d'étude sur Bangkok (Thaïlande). Pour cette métropole, nous dispo- 
sons des données de recensement de la population à l'échelle de l'arrondissement, ainsi 
que d’une couche représentant les bâtiments en deux dimensions. Nous décrivons 
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ensuite une seconde méthode de ventilation des données de population sur un autre 
site, Rouen (France). Nous montrons l’utilisation d’algorithmes qui permettent 
d'estimer la répartition de la population par bâtiment, ainsi que la description réaliste 
de la composition de chaque ménage. Ce dernier aspect peut être tout particulière- 
ment intéressant dans le contexte de la dengue pour laquelle il existe un lien fort entre 
taux d’immunité de la population et pyramide des âges au sein de la population. Les 
résultats de ces algorithmes peuvent également représenter une source de données 
utilisable dans le cadre de modèles de simulation (chapitre 9). 


» Désagréger et redistribuer la population 


Cartographie dasymétrique classique 


Utilisée dès 1911 par Semenov-Tian-Shansky (Petrov, 2012) et popularisée par Wright 
en 1936 (Li et Weng, 2005), la cartographie dasymétrique (littéralement « mesure de 
densité ») consiste à distribuer la population connue depuis une zone source vers un 
ensemble de zones cibles sur la base de caractéristiques thématiques. Historiquement, 
il s'agissait de redistribuer la population d'une zone administrative dans les zones 
habitées de cette même entité spatiale. Ce principe a été progressivement affiné par 
l'emploi des régressions linéaires multiples entre la population enregistrée par maille 
et l’utilisation/occupation du sol de cette même maille (équation 1). 


P, = Bo + BU, + BU, +... + BU, (équation 1) 


Cette méthode propose de décrire la distribution de la population P, d'une zone source 
à l’aide d'une régression multiple où à chaque type d'utilisation du sol U, est associé un 
coefficient f,,. Ces coefficients, qui représentent les densités de population associées 
aux différentes classes d'utilisation du sol, sont ensuite utilisés comme pondérateurs 
pour désagréger la population source vers chaque zone cible, en fonction de sa surface. 
Idéalement, le coefficient de corrélation doit tendre vers 1. 


Cette méthode a donc été utilisée pour estimer la population de Bangkok a échelle fine 
(Misslin et Daudé, 2017) a partir des images satellite OLI et TIRS (Landsat 8). La part 
de bati lache, dense et trés dense au sein des sous-districts de la capitale thailandaise a 
été estimée a partir d’une classification non supervisée. La cartographie dasymétrique 
dans ce contexte repose sur l’hypothése qu'un pixel de bâti dense accueillera plus de 
population qu'un pixel de bâti épars. Chaque pixel se voit donc attribuer un coefficient 
(£) obtenu à partir de la régression multiple sur l'ensemble de la zone et correspondant 
à sa catégorie (U), que Ton multiplie par la population de la zone géographique dans 
laquelle il se situe (figure 4.1). 


Comme les images satellite Landsat 8 sont disponibles pour l'ensemble de la planète 
et que la méthodologie employée est éprouvée (Hallot et al., 2019), cette approche a 
l'avantage d’être théoriquement reproductible en différents lieux, dès lors qu'il existe des 
données de recensement dans la zone étudiée. D'autres données thématiques peuvent 
également être utilisées pour caractériser le lien avec la distribution de la population, 
telles que le taux d'imperméabilisation des sols et leur surface relative (Hallot et al., 2019) 
ou encore les données de luminosité artificielle du ciel nocturne (Briggs et al., 2007). 


Néanmoins, toute classification non supervisée reste imparfaite et l'apparition de 
nouveaux services proposant des données relatives à l'emprise spatiale des bâtiments 
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et à leur utilisation (Google Maps, ESRI!, etc.) ouvre de nouvelles perspectives. Il s'agit 
alors de dénombrer les bâtiments habités et de prendre en compte leur emprise spatiale 
par zone de recensement pour appliquer les mêmes approches par régression linéaire. 
Les données issues du projet OpenStreetMap (OSM), par exemple, sont disponibles 
gratuitement et renseignent les types d'occupation et d’utilisation du sol d'une zone 
géographique donnée et, parfois, l'empreinte géographique des bâtiments. Mais le niveau 
dexhaustivité et de véracité des données de ce projet participatif dépend grandement du 
lieu étudié. D’une grande précision en France, ces informations ne sont pas aussi précises 
dans d’autres zones du monde, notamment à Bangkok. Il est cependant possible, avec la 
démultiplication des producteurs de données géospatiales, de contourner cet obstacle. 


Nombre d’habitants 


Fréquence (nb de pixels)* 


100 000 


50 000 
10km 


(0) 
110 20 30 40 50 60 70 80 90100 
Nombre d'habitants 


Données : NSO (2010) et Landsat8 OLI et TIRS (04/01/2015) * Les pixels non habités ne sont pas comptabilisés 


Figure 4.1. Répartition de la population à Bangkok sur une maille de 30 x 30m par cartographie 
dasymétrique. Source : Misslin et Daudé (2017). 


Affiner les données thématiques pour améliorer 
la cartographie dasymétrique 


Les méthodes de webscrapping permettent d'accéder à des bases de données en ligne 
(privées ou non), comme celles issues du WMS (Web Map Service) d'ESRI. Ces bases de 
données permettent d'obtenir une carte thématique décrivant l'emprise spatiale des bâti- 
ments (figure 4.2a). À Bangkok, une comparaison visuelle entre les images satellite et les 


1. https://www.esrifrance.fr/ 
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données représentant les bâtiments issues du WMS ESRI montre une bonne représenta- 
tivité qualitative de ces dernières. Un code a donc été développé pour extraire l’ensemble 
de ces informations’, tout d'abord sous forme de raster (figure 4.2b), puis de façon 
vectorisée? (figure 4.2c et d), formant une base de données de 1,4 million de polygones. 


Figure 4.2. Processus d’extraction de batiments a partir d'une image du WMS d’ESRI. 


(a) La base de données brutes, (b) le changement d'intensité des couleurs, (c) l'extraction des polygones bâtis 
et (d) leur projection sur l'image satellite pour comparaison visuelle. 


Cette couche 2D des batiments de Bangkok* permet, pour chaque zone (ici les sous- 
districts), de calculer un nombre de bâtiments, un périmètre cumulé de ces bâtiments 
et une surface de bâti. Bangkok est également une ville où se mêlent des quartiers 
sinueux faits de constructions de un à deux étages et de grands buildings d'habitation. 
Lajout de données sur l'estimation de la hauteur peut donc apporter une informa- 
tion supplémentaire utile pour désagréger la population. Pour ce faire, nous avons 
collecté les données de la mission SRTM (Shuttle Radar Topography Mission) qui 
permettent d'obtenir des modèles numériques d'élévation (MNE) avec une résolution 
au sol d'environ 30m à l'équateur’. Un modèle numérique de terrain (MNT) a ensuite 
été construit en interpolant l’altitude moyenne des canaux et rivières de la ville pour 
définir un niveau de base que nous avons soustrait au MNE, afin d'obtenir l'élévation 
localef. Enfin, les zones habitées ont été distinguées des zones non habitées grâce à 
la couche d'utilisation du sol issue d'OSM. Celle-ci a permis de supprimer les bâti- 
ments situés dans des zones non résidentielles (temples, zones industrielles, centres 
commerciaux, zones de loisirs, etc.) ainsi que les zones végétalisées ou en eau. 


À partir de ces données, sont calculées plusieurs régressions linéaires entre la popula- 
tion à Bangkok par sous-districts et les informations sur l’utilisation du sol : le nombre 
de bâtiments, la superficie de l'emprise spatiale du bâti, le linéaire de bâti et l'altitude 
moyenne dans chacun de ces sous-districts (figure 4.3a). Nous sélectionnons ensuite la 
combinaison de facteurs issue de l'équation 1 qui permet d'obtenir une erreur moyenne 
absolue la plus faible possible et le meilleur coefficient de corrélation. Ce dernier est 
obtenu en utilisant tous les indicateurs de bâti (figure 4.3b), mais la combinaison dau 
moins deux paramètres (peu importe lesquels) donne des coefficients de corrélation 
R’ supérieurs à 0,81. En effet, au regard de la figure 4.3a la plupart des paramètres qui 
caractérisent le bâti sont très corrélés avec la population (R° > 0,75). 


2. À noter que la méthode de webscrapping employée est désormais légale aux yeux de la Cour de justice 
américaine, https://techcrunch.com/2022/04/18/web-scraping-legal-court/. 

3. https://drive.univ-rouen.fr/f/a20e1dc3baad49caadec/?dl=1 

4. https://drive.univ-rouen.fr/f/d71c73ffd3504ae28068/?dl=1 

5. https://www2.jpl.nasa.gov/srtm/ 

6. https://drive.univ-rouen.fr/f/f282bacb588e43818c75/?dl=1 
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a Nombre Bâtiment 
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Figure 4.3. Liens entre les caractéristiques des bâtiments d'une zone géographique et la 
population associée. 
(a) Corrélations (R°) entre paires de paramètres et (b) régression entre la population connue au sous-district 


et la population estimée avec l'ensemble des indicateurs d'occupation du sol (droite noire, la droite bleue 
correspond à la relation y = x). 


Nous procédons ensuite en trois étapes : tout d’abord nous créons une grille de maille 
d'une taille désirée, par exemple 50 m, à laquelle nous inférons les informations sur 
les bâtiments (le nombre, la longueur, la hauteur et la surface de bâti dans chaque 
maille) ; ensuite nous appliquons à chacune des mailles les coefficients de la régression 
linéaire qui lient la population à ces indicateurs à l'échelle du sous-district; et enfin, 
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nous ajustons la population au sein de chaque maille selon le sous-district de sorte à 
conserver la population des arrondissements du recensement. Nous obtenons alors 
une ventilation de la population dans une maille de résolution donnée (figure 4.4b). 
Néanmoins, comme le montre la figure 4.4c, selon les modèles de régression multiple 
utilisés en entrée, des différences locales importantes peuvent survenir. Elles se situent 
essentiellement dans des zones où la hauteur des bâtiments est importante, ce qui 
implique que les régressions qui prennent en compte l'élévation vont tendre à concen- 
trer la population dans ces mailles, à l'inverse des autres approches qui vont répartir la 
population de manière plus homogène. 


La cartographie à haute résolution spatiale des densités de population peut donc aider à 
la décision (contrôle vectoriel accru dans les zones les plus denses) ou être utilisée comme 
donnée d'entrée pour différents modèles (voir chapitre 9). Pour ces derniers, les évolu- 
tions les plus récentes conduisent à aller encore plus loin dans la qualification des données 
de population à haute résolution spatiale. Il s’agit non plus seulement de décompter le 
nombre de personnes présentes dans une zone, mais de caractériser chacun des indi- 
vidus présents dans la zone. Nous illustrons cette démarche sur la commune de Rouen. 


» Population synthétique, une méthodologie au service 
d’une analyse spatiale fine des problématiques de santé 


La construction d'une population synthétique consiste à concevoir un ensemble 
exhaustif et statistiquement représentatif d'individus et de ménages sur un territoire 
donné (Lovelace et Dumont, 2016). Plusieurs méthodes (reconstruction synthé- 
tique, optimisation combinatoire, apprentissage statistique) ont été développées ces 
dernières années pour générer une population d'agents avec des attributs démo- 
graphiques réalistes. La disponibilité des données et le niveau de précision de celles-ci 
vont participer au choix de la méthodologie à employer : avec ou sans échantillon, 
déterministe ou stochastique. 


Principe des générateurs de populations synthétiques 


Lalgorithme le plus utilisé pour générer une population synthétique est ’IPF, pour 
Iterative Proportional Fitting (H6rl et Balac, 2021). Cet algorithme utilise deux types 
de données : un échantillon représentatif d'individus décrit selon plusieurs moda- 
lités et le nombre total d'individus dans la population par modalité. L’algorithme va 
alors générer une population d'individus équivalant à la population globale tout en 
respectant la distribution relative de l'échantillon par modalités croisées. LIPF n'est 
cependant pas adapté pour respecter des contraintes d’emboitement, par exemple 
un nombre de ménages composés de n individus. Cet algorithme vise en effet une 
répartition optimale à l'échelle soit des ménages soit des individus (Khachman et al., 
2021). Cette perte d’information peut être problématique notamment pour la prise 
en compte des mobilités des individus, celles-ci dépendant à la fois des caractéris- 
tiques individuelles et du ménage. D’autres méthodes ont été développées pour pallier 
cette perte comme l’Iterative Proportional Update (Ye et al., 2009) ou le Hierarchical 
Iterative Proportional Fitting (Yameogo et al., 2021). Ces algorithmes prennent donc 
en considération non seulement les informations au niveau individuel, mais aussi au 
niveau des ménages dans leurs processus de reconstruction synthétique. 
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De la population aux ménages, du ménage à l'individu 


La France, via les recensements de la population, dispose d'informations très précises 
sur les caractéristiques de la population : âge, profession, moyens de transport utilisés, 
conditions de logement, etc. Ces données (2017) sont ici mobilisées pour générer 
une population synthétique à partir de l'algorithme IPF sur l'ensemble des IRIS 
(îlots regroupés pour l'information statistique) de la commune de Rouen. Au total, 
62121 ménages et 110 142 individus ont été générés sur l'ensemble des 42 IRIS de la 
zone d'étude. La seconde étape dans ce processus de génération consiste à ventiler les 


e Recensement 
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Figure 4.5. Répartition des ménages sur trois IRIS à Rouen. (a) Répartition spatiale. (b) Structure 
des groupes sociaux : (1) ménage d’une seule personne, (2) famille monoparentale, (3) couple 
sans enfant, (4) couple avec enfant et (5) ménages complexes. 
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individus et les ménages dans les bâtiments de chacun des IRIS. Il faut d’abord estimer 
le nombre d'habitants à partir des caractéristiques des bâtiments résidentiels de la 
BD TOPO’ de l'IGN. Le nombre d'habitants est donc fonction du nombre d'étages, 
de la superficie de chaque étage au prorata de la surface totale des bâtiments résiden- 
tiels. Les ménages sont ensuite ventilés dans les bâtiments en fonction du nombre 
d'individus qui les composent et de la capacité d’accueil théorique des logements. 


La figure 4.5 présente quelques résultats de cette génération de population synthétique 
sur trois IRIS de la commune de Rouen (Vieux-Marché — Palais de Justice, Vieux- 
Marché — Saint-Patrice, Beauvoisine). La carte représente le nombre de ménages 
dans chaque bâtiment et l’histogramme présente la distribution des types de ménage 
dans les trois IRIS. La même méthode peut ensuite être employée pour ventiler la 
population en fonction de la pyramide des âges. 


» Conclusion 


La connaissance de la répartition des populations doit permettre de limiter les risques 
d'exposition des humains aux maladies à transmission vectorielle. Le risque épidé- 
mique résulte en effet de la coprésence spatiale et temporelle des hôtes et des vecteurs. 
Nous avons montré dans ce chapitre quelques-unes des méthodes qui permettent, 
depuis diverses sources de données, de désagréger des données issues de recensements 
qui sont trés souvent disponibles 4 un maillage administratif assez peu pertinent pour 
des enjeux épidémiologiques. Si la spatialisation des individus est à ce jour possible à 
une très haute résolution spatiale, l’un des enjeux à venir est de s'attaquer à l'évolution 
temporelle de ces spatialités. Les données sociodémographiques que nous avons utili- 
sées ici reposent en effet sur des recensements au lieu de résidence, elles ne permettent 
donc pas d'évaluer les risques d'exposition des individus lors de leurs déplacements. 
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Chapitre 5 


Texture des images satellite 
et caractérisation des milieux urbains 
favorables aux moustiques vecteurs 


CLAIRE TEILLET, OPHÉLIE HOARAU, NAUSICAA HABCHI-HANRIOT, 
BENJAMIN PILLOT, THIBAULT CATRY, ANNELISE TRAN 


La texture est une notion importante du traitement d'images, utilisée lors des étapes de 
classification, de segmentation ou de synthèse des images (Busch et al., 2004). Linforma- 
tion texturale est souvent utilisée en complémentarité avec les informations spectrales 
pour décrire une image (Pacifici et al., 2009) en la décomposant en fréquences et orien- 
tations. Haralick (1979) a montré qu’une texture peut être représentée selon deux 
dimensions : la premiére correspond 4 la description du motif en termes de fréquence 
et d'orientation, la seconde correspond à l’organisation spatiale de ces motifs (Haralick, 
1979). La notion de texture, bien que trés dépendante des cas d’application et difficile- 
ment généralisable, renvoie d’abord à une recherche d’homogénéité au sein des variations 
de tons de gris via l'identification de motifs dans l’image (Caloz et Collet, 2001). 


Les principales notions liées à l’analyse de texture sont illustrées dans la figure 5.1. 
La répétition réguliére d’un motif dans une ou plusieurs directions correspond a 
une texture anisotrope, tandis qu'une répartition aléatoire des pixels dans toutes les 
directions est considérée comme une texture isotrope, pour laquelle on observe une 
certaine homogénéité de la structure. Une texture anisotrope est donc caractérisée par 
des propriétés dépendantes de la direction, alors que l’isotropie est indépendante de la 
direction (figure 5.1a). Une micro-texture et une macro-texture représentent respecti- 
vement les détails trés fins et les détails les plus grossiers d’une texture (figure 5.1b). De 
plus, la texture varie avec l'échelle d’observation, une fenêtre d’analyse d'un objet doit 
étre suffisamment grande pour inclure le motif définissant la texture, mais suffisam- 
ment petite pour conserver la dépendance spatiale et permettre d’identifier la texture 
correspondant à l’objet (figure 5.1c). Lors de changements d'échelle et donc de fenêtre 
d'analyse, un motif peut être perdu mais d’autres motifs peuvent être mis en évidence. 
Ces notions d’échelles d'analyse de la texture font écho aux résolutions spatiales des 
images d'observation de la Terre utilisées en entrée des méthodes d'extraction de 
l'information texturale (figure 5.1d). 


L'analyse de texture a été largement utilisée pour le traitement et l'extraction d'infor- 
mation des images satellite, en particulier dans des études sur les milieux forestiers 
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et urbains (Lefebvre, 2011; Pandey et al., 2021; Proisy et al., 2007; Puissant et al., 
2006 ; Ruiz Hernandez et Shi, 2018). En milieu urbain, la densité et la proportion de 
végétation et de bâti, ainsi que les arrangements spatiaux entre bâti et végétation, sont 
les principales composantes de la texture. L'information texturale a jusqu'à présent 
très rarement été exploitée dans le cadre d'études sur les maladies vectorielles dues 
à certaines espèces de moustiques (cf. chapitre 1), alors qu’elle pourrait compléter 
les indicateurs environnementaux, climatiques et démographiques présentés dans 
les chapitres précédents, en particulier pour les moustiques vecteurs de la dengue 
Ae. aegypti et Ae. albopictus, principalement présents en milieu urbain. 


#7 Ke 
anisotrope 


Micro-texture Macro-texture 
isotrope anisotrope isotrope isotrope 


Sol isotrope 


Echell 
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Figure 5.1. Notions clefs pour l'analyse de texture d'images satellite. 


Dans ce chapitre, différentes applications d'analyse texturale sont présentées. En parti- 
culier, l'implémentation de l'approche FOTO (Couteron et al., 2006) sur des images 
Pléiades à très haute résolution spatiale a permis d'étudier les relations entre des 
variables urbaines issues des images, la distribution des cas de dengue et la disponibilité 
en gîtes larvaires pour des applications au Brésil et à La Réunion. 


» Différentes méthodes pour caractériser la texture d'une image 


Il existe de nombreuses méthodes d'analyse de texture, dont l'approche statistique et 
l'approche fréquentielle. 


Un exemple d'analyse statistique : la matrice de cooccurrence 
de niveaux de gris et les indices d’Haralick 


Il s’agit d'une approche basée sur la représentation des niveaux de gris en statistique 
qui consiste à calculer une matrice de cooccurrence de ces niveaux (GLCM pour Grey 
Level Co-occurrence Matrix). Cette matrice est la représentation dans une image du 
nombre d'occurrences de paires de niveaux de gris séparées par une certaine distance 
dans une direction donnée (figure 5.2a). Haralick (1979) propose le calcul d'indices 
à partir de cette matrice, dont les plus pertinents pour les études du milieu urbain 
sont : l'énergie, l’entropie, le contraste, la corrélation et ’homogénéité (Maillard, 
2003; Pacifici et al., 2009). 
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Un exemple d'analyse fréquentielle : la transformée de Fourier 


Fourier a démontré au x1x* siècle que les fonctions périodiques peuvent être décom- 
posées selon une somme de fonctions sinusoïdales et co-sinusoïdales. Un signal 
discret peut ainsi être représenté par une fonction composée d’une fréquence fonda- 
mentale (la plus basse), des harmoniques (les multiples) et des coefficients de la série 
de Fourier. Ce sont ces coefficients qui permettent de comprendre la contribution de 
chaque fréquence à la formation du signal (Caloz et Collet, 2001). 


La représentation fréquentielle d'une transformée de Fourier est un spectre d'amplitude 
correspondant à des valeurs d'énergie pour chaque fréquence (figure 5.2b). Le spectre 
rend donc compte de la distribution énergétique de l’image en respectant la périodi- 
cité et l'orientation de la texture. Plus on s'éloigne du centre, plus la fréquence observée 
est élevée. Une texture lisse aura au centre un spectre correspondant à des valeurs de 
basses fréquences. À l'inverse, une texture rugueuse aura des valeurs plus étalées dans le 
spectre et correspondra à des hautes fréquences (Regniers, 2014). Les transformées de 
Fourier permettent ainsi de passer du domaine spatial au domaine fréquentiel à partir de 
la décomposition fréquentielle des signaux qui constituent une texture. 


Lalgorithme FOTOTEX repose sur cette approche fréquentielle par transformée de 
Fourier. Les détails de la méthode et le principe de sa mise en œuvre sont présentés 
dans la section suivante. 
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Figure 5.2. Exemple d'analyse de texture par (a) une analyse statistique grace à une matrice 
de cooccurrence de niveaux de gris et aux indices d’Haralick et (b) une analyse fréquentielle 
a travers les spectres de la transformée de Fourier en 2D pour : (1-1’) prairies, (2-2’) céréales, 
(3-3’) prairies + céréales. Sources : (a) Brynolfsson et al. (2017); (b) Lefebvre (2011). 
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» Étude des relations entre variables urbaines et distribution 
des cas de dengue à Brasilia par une approche texturale 


Contexte 


Dans le cadre du projet Apureza (Analyse par télédétection des relations entre 
paysages urbains dengue et Zika) financé par le Cnes (2017-2020), un algorithme de 
caractérisation multi-échelle de la texture des zones urbaines a été développé. Il vise 
à évaluer l'apport d'une approche texturale fréquentielle à l'étude des relations entre 
variables urbaines et distribution des cas de dengue, en particulier dans la ville de 
Brasilia, quartier de São Sebastião (Teillet et al., 2021). 


Ces dernières années, une très forte augmentation des cas de dengue a en effet été 
observée à Brasilia : en 2004-2005 il y avait 200 cas confirmés alors qu'en 2015 leur 
nombre atteignait 11722. 


Données 


Pour la caractérisation de l'emprise urbaine, une image Sentinel-2 à haute résolution 
spatiale (HR) du 31 août 2018 a été utilisée, téléchargée depuis la plateforme de mise 
à disposition des produits Sentinel du programme Copernicus! de PESA (Agence 
spatiale européenne). Les images Sentinel-2 comptent treize bandes spectrales pour 
une résolution spatiale comprise entre 10 et 60 m. 


Pour la caractérisation de l’intra-urbain, une image Pléiades à très haute résolution 
spatiale (THR) du 14 janvier 2020 a été fournie par l'équipe GEOSUD de Montpellier), 
Elle est composée de plusieurs bandes : une bande panchromatique P (480-830nm) 
à 0,7m de résolution rééchantillonnée à 0,5m et des bandes multispectrales (MS) 
à 2,8 m, rééchantillonnées à 2m, incluant les canaux bleu (430-550 nm), vert (490- 
610 nm), rouge (600-720 nm) et proche infrarouge (750-950 nm). 


Méthode 


La méthodologie développée utilise des informations de texture extraites de ces 
données d'observation de la Terre acquises à différentes résolutions spatiales pour 
caractériser les paysages urbains des pays du Sud à trois échelles différentes, définies 
dans le cadre du centre d'expertise scientifique « Urbain » de Theia? : 

— l'empreinte urbaine (macro-échelle) ; 

— l'échelle des unités urbaines (méso-échelle) ; 

— l'échelle des objets tels que les bâtiments (micro-échelle). 


Pour cela, l'algorithme FOTOTEX est utilisé et se fonde sur l'algorithme «FOTO » 
(pour Fourier-based Textural Ordination), qui visait à l’origine à caractériser les 
modèles de végétation des écosystèmes tropicaux (Couteron et al., 2006; Proisy et al., 
2007). Afin de disposer d'un outil simple, robuste et efficace, cette version de l’algo- 
rithme a été développée en Python et mise à disposition en ligne‘. Lalgorithme a été 


1. https://scihub.copernicus.eu/dhus/#/home 
2. https://ids-dinamis.data-terra.org/ 

3. https://www.theia-land.fr/ceslist/ces-urbain 
4. https://framagit.org/espace-dev/fototex 
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optimisé en matière de processus et de temps de calcul. Les détails de l’utilisation et de 
la paramétrisation de l’algorithme FOTOTEX sont disponibles dans l’article de Teillet 
et al. (2021). La méthode est présentée ici de façon synthétique. 


Étape 1 : Partitionnement de l'image 


La première étape de la méthode consiste à partitionner l’image (figure 5.3a). Ce parti- 
tionnement correspond à la définition de la règle de parcours d’une fenêtre sur l’image. 
Deux méthodes de découpage sont disponibles : la méthode par bloc (méthode block), 
qui analyse la texture pour chaque fenêtre de l’image bloc par bloc, et la méthode 
par déplacement (méthode moving window), qui analyse la texture pour une fenêtre 
donnée en glissant d'ouest en est et du nord au sud avec un décalage d’un pixel. 


Étape 2 : Analyse spectrale par transformée de Fourier 
et calcul des «r-spectres » 


Une transformée de Fourier rapide (FFT) est calculée pour chaque fenêtre d'analyse 
de l’image originale et permet d'exprimer la variance de l’image sous la forme 
d'une somme pondérée de formes d'ondes cosinus et sinus dont la direction et la 
fréquence spatiale varient (Couteron et al., 2006; Lang et al., 2018; Proisy et al., 
2007; figure 5.3b). Les coefficients de pondération quantifient la contribution de 
chaque fréquence et chaque direction à la formation du signal. Un « périodo- 
gramme », calculé pour chaque fenêtre, exprime la proportion de variance expliquée 
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a) Découpage de l’image -------- + des fréquences spatiales contenues dans chaque --- 
fenêtre et regroupement dans un r-spectre 
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Figure 5.3. Cadre méthodologique simplifié de l'algorithme FOTOTEX. 
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par chaque paire de fréquences (p, q), où p et q sont des fréquences spatiales dans 
les directions des colonnes et des lignes de l’image (Proisy et al., 2007). Finalement, 
un r-spectre isotrope est calculé moyennant les valeurs du périodogramme dans 
toutes les directions (de 0° à 360°) pour obtenir un spectre radial moyen, appelé 
«r-spectre», noté I(r), pour chaque fenêtre d'analyse. Un tableau synthétisant 
l'information des r-spectres est produit; dans ce tableau, chaque ligne correspond 
au r-spectre d'une fenêtre et chaque colonne correspond à une fréquence spatiale 
(Lang et al., 2018). 


Étape 3 : Analyse en composantes principales 


Une analyse en composantes principales (ACP) est appliquée au tableau des r-spectres 
afin de résumer l'information qu'il contient (figure 5.3c et d). Cette méthode statis- 
tique consiste à transformer des variables corrélées entre elles en nouvelles variables, 
nommées «composantes principales» ou «axes principaux», décorrélées les unes 
des autres — l'information est ainsi réduite, car les premiers axes suffisent à repré- 
senter une grande proportion de la variabilité totale. Les trois premiers axes résultant 
de ACP sont des indices texturaux qui peuvent ensuite être visualisés sous la forme 
d'une image composite rouge-vert-bleu (RVB) [Proisy et al., 2007]. 


Brésil 


Brasilia 


San Sebastiad 


H = =3 im 


Figure 5.4. Distribution des clusters de cas de dengue à São Sebastião (Brasília) en fonction des 
types de quartiers identifiés par la texture des images Pléiades. 

(a) Localisation de la zone d'étude, le quartier de São Sebastião, Brasilia (Brésil). (b) Clusters de cas de 
dengue identifiés par SATSCAN pour la période 2007-2017. (c) Composition colorée RVB en sortie de 
FOTOTEX. (d) Organisation du bâti et de la végétation sur une image Pléiades définissant la typologie 
urbaine correspondant à chaque cluster. 
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Résultats 


La figure 5.4 présente les résultats obtenus sur São Sebastião, quartier de Brasilia au 
Brésil (figure 5.4a). L'information de texture issue de FOTOTEX permet d'une part de 
délimiter l'emprise urbaine (échelle macro) à partir de l’image Sentinel-2 et, d'autre 
part, d'identifier des quartiers représentés par des informations de textures différentes 
au sein de cette zone urbaine (figure 5.4b). 


Ces quartiers sont liés à une distribution spatio-temporelle hétérogène des cas de dengue 
(figure 5.4c). En effet, des clusters de cas ont été identifiés à l'aide du logiciel SATSCAN?® 
sur les données de cas de dengue géolocalisés dans la zone pour la période 2007-2017 
issus de la base SINAN*®. Chaque cluster est caractérisé par une configuration spécifique 
du paysage urbain, avec une organisation spatiale particulière entre le bâti (structure, 
densité, hauteur) et la végétation urbaine notamment (figure 5.4d). 


Des analyses statistiques ont montré que la densité du bâti et la proportion de végétation, 
facteurs principaux composant l'information de texture en milieu urbain, sont corré- 
lées positivement avec le nombre de cas de dengue à Säo Sebastiäo. Des corrélations 
positives existent également avec d’autres variables comme la température de surface. 


» Cartographie de La distribution des gîtes larvaires potentiels 
du moustique-tigre à La Réunion 


Contexte 


Dans le cadre du projet Anisette (Analyse inter-site : évaluation de la télédétection 
comme outil prédictif pour la surveillance et le contrôle de maladies causées par des 
moustiques), l’algorithme FOTOTEX décrit précédemment a été appliqué sur l’île de 
La Réunion. Ce département français de l'océan Indien présente un relief accidenté 
(entre 0 et 3070 m) qui contraint à concentrer l’habitat sur les zones littorales. Le 
moustique-tigre, Aedes albopictus, y est le principal vecteur d’arbovirus, comme ceux 
de la dengue (épidémies en 1977, 2004 et depuis 2016) et du chikungunya (épidémie en 
2005-2006). Dans cette étude, les relations entre les indices texturaux obtenus et des 
données de terrain sur le nombre de gites larvaires du moustique-tigre ont été analy- 
sées, ceci afin de cartographier les zones propices a la prolifération d’Ae. albopictus, a 
une échelle régionale (Hoarau, 2021). 


Données 


Données sur les gites larvaires du moustique-tigre 


Dans le cadre de ses actions de lutte antivectorielle (LAV), l’agence régionale de santé 
(ARS) de La Réunion collecte des données sur les gites larvaires en procédant 4 une 
prospection des jardins par porte-a-porte (figure 5.5). Pour chaque zone de surveil- 
lance entomologique (ZSE), constituée d’environ 200 maisons, les données collectées 
sont le nombre et le type de gites potentiels. Pour cette étude nous disposions des 
données de 346 ZSE (sur un total de 1203 suivies par l’ARS, soit un quart) collectées 
entre janvier 2020 et mars 2022. 


5. https://www.satscan.org/ 
6. https://portalsinan.saude.gov.br/ 
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Figure 5.5. Typologie et répartition des gîtes larvaires du moustique-tigre à La Réunion, 2021- 
2022. 


(a) Nombre de gîtes potentiels observés par quartier. (b) Répartition selon les types de gîtes observés. 
Source : (a) données ARS La Réunion. 


Données environnementales 


La mosaïque Pléiades 2020 de La Réunion est un produit généré à partir de douze 
images diffusées dans le cadre de Kalideos’. Les caractéristiques de ces images ont été 
précisées précédemment. 


Les données sur le cumul annuel de pluies et la moyenne annuelle de température 
obtenues par interpolation des données du réseau de mesure Météo-France/Cirad ont 
été téléchargées depuis le portail de l'Atlas web agricole pour la recherche, 


Méthodes 


Calcul des indices texturaux 


Après différents tests sur des extraits d'images, l'algorithme FOTOTEX a été appliqué 
selon les étapes décrites précédemment sur les zones bâties de la mosaïque Pléiades 
(figure 5.6a) avec une taille de fenêtre carrée de 67 pixels de côté (soit 33,5 m) et par la 
méthode de découpage par bloc (figure 5.6b). 


7. https://www.kalideos.fr 
8. https://aware.cirad.fr/ 


80 


Texture des images satellite et caractérisation des milieux urbains favorables aux moustiques vecteurs 


Calculs des indices spectraux 


Différents indices ont été calculés à partir de la mosaïque Pléiades : des indices de 
végétation (figure 5.6c) et d'eau à différence normalisée (ces indices sont décrits dans 
le chapitre 2), ainsi que l'indice de brillance (IB) calculé à partir des bandes rouge et 
proche-infrarouge et qui permet de mettre en évidence les surfaces réfléchissantes 
comme les zones de bâti, routes, sols nus, etc. (figure 5.6d). 


Analyses statistiques 


Les corrélations entre le nombre de gites larvaires potentiels observés (variable a 
expliquer) et la moyenne pour chaque quartier des valeurs des indices texturaux, 
spectraux, des précipitations et températures annuelles (variables explicatives) ont été 
étudiées par une analyse univariée (en testant les variables une à une), puis par une 
analyse multivariée, combinant les variables les plus significatives, à l’aide d'un modèle 
linéaire généralisé (GLM). 


Nb gîtes estimé 
L 0-50 
E 51-100 


Élevé Élevé 


- Faible - Faible 


Figure 5.6. Illustration d'indices texturaux et spectraux calculés à partir d'une mosaïque 
Pléiades, île de La Réunion, 2020, et cartographie prédictive du nombre de gîtes larvaires 
potentiels du moustique-tigre à l'échelle des zones de surveillance entomologique. 


(a) Mosaïque Pléiades 2020. (b) Composition colorée RVB en sortie de FOTOTEX. (c) Indice de végétation 
(NDVI). (d) Indice de brillance (IB). (e) Nombre de gîtes larvaires potentiels estimé par ZSE. 
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Résultats 


Comme pour les sites d'étude au Brésil, la visualisation des indices texturaux issus de 
l'algorithme FOTOTEX met en évidence une typologie de quartiers (figure 5.6b) : 

— quartiers avec une densité de bâtiments très élevée et peu de végétation; 

— quartiers avec des bâtiments plus distants, séparés par de la végétation; 

— quartiers avec une densité de bâtiments très faible avec une forte présence de 
végétation; 

— zones industrielles caractérisées par la présence de grands bâtiments, de larges 
surfaces imperméabilisées et peu de végétation. 


Ces indices, combinés aux autres variables environnementales (indices spectraux, 
précipitations, température) ont été identifiés comme significativement corrélés aux 
données de terrain sur les gîtes larvaires d’Ae. albopictus (Hoarau, 2021). La relation 
obtenue permet ensuite de prédire, pour tous les quartiers de l’île, le nombre de gîtes 
larvaires potentiels, en fonction des variables environnementales (figure 5.6e). 


» Conclusion 


Les résultats présentés dans ce chapitre, basés sur l’utilisation de FOTOTEX, illustrent 
l'apport de l'indication de texture pour la caractérisation des zones urbaines. La 
texture fournit une information complémentaire à celle issue des indices spectraux 
ou de l'occupation du sol. En particulier, il est possible de décrire l’urbain à différentes 
échelles et notamment d'identifier des unités urbaines caractérisées par des typologies 
spécifiques (structures spatiales hétérogènes au niveau du bâti et de la végétation). Ces 
informations se révèlent particulièrement pertinentes pour l'étude des relations entre 
les paysages urbains et la dynamique spatio-temporelle de la dengue à partir d’image- 
ries à très haute résolution spatiale (Pléiades, SPOT 6/7), comme nous l'avons montré 
au Brésil et à La Réunion. 


Dans ces deux zones, des analyses statistiques ont permis de mettre en évidence que 
le signal de texture des unités urbaines est fortement corrélé à la distribution spatiale 
des cas de dengue (Brésil) et représente un bon indicateur du nombre de gîtes larvaires 
potentiels dans les zones urbaines (La Réunion). La mise en évidence de ces corréla- 
tions constitue un apport important dans l'étude des dynamiques de populations de 
moustiques et du risque de transmission de la dengue en contexte urbain. En effet, la 
texture contenue dans les images satellite à très haute résolution spatiale, en prédisant 
assez bien le nombre de gîtes larvaires potentiels, permet de s'affranchir du manque 
de données collectées in situ. 


Ces jeux de données, dont l'apport est précieux pour la calibration et la validation 
des sorties de modèles de cartographie prédictive des densités de moustiques comme 
Alborun (Tran et al., 2020) ou Arbocarto (voir chapitre 8), sont rares, ce qui représente 
souvent une limite forte des applications en santé. Il serait donc possible d'utiliser l’infor- 
mation de texture comme une variable d'entrée des modèles pour des applications sur 
des sites où les données de gîtes larvaires ne sont pas disponibles. De même, les liens 
entre texture et dynamique spatio-temporelle des cas ouvrent des perspectives dans 
l'étude du risque de transmission de la dengue à partir de l'imagerie satellite. Il s’agit 
cependant d'un mécanisme très complexe faisant intervenir de nombreuses variables, 
rendant indispensable l’utilisation de données complémentaires à l'imagerie satellite. 


82 


Texture des images satellite et caractérisation des milieux urbains favorables aux moustiques vecteurs 


» References 


Brynolfsson D. Nilsson D., Torheim T., Asklund T., Karlsson C.T., Trygg J., Nyholm T., Garpebring A., 
2017. Haralick texture features from apparent diffusion coefficient (ADC) MRI images depend on 
imaging and pre-processing parameters. Scientific Reports, 7 (1), 4041. 


Busch A., Boles W., Sridharan S., 2004. Logarithmic quantisation of wavelet coefficients for improved 
texture classification performance, in IEEE International Conference on Acoustics, Speech, and Signal 
Processing, 17-21 May 2004, Montreal, QC, Canada, IEEE, iii-569. 


Caloz R., Collet C., 2001. Précis de télédétection — Volume 3 : Traitements numériques d'images de télé- 
détection, Québec, Presses de l’Université du Québec — Agence Universitaire de la Francophonie, 386 p. 
(coll. Universités francophones). 


Couteron P, Barbier N., Gautier D., 2006. Textural Ordination Based on Fourier Spectral Decomposition: 
A Method to Analyze and Compare Landscape Patterns. Landscape Ecology, 21, 555-567. 


Haralick R., 1979. Statistical and Structural Approaches to Texture. Proceedings of the IEEE, 67 (5), 
786-804. 


Hoarau O., 2021. Télédétection et modélisation de populations de moustiques vecteurs : Extraction 
d'indicateurs paysagers à partir d'images d’Observation de la Terre pour lestimation de la distribution 
des gîtes larvaires en milieu urbain, mémoire de master 2, spécialité Sciences de la terre, des planètes et 
de l'environnement, Université de La Réunion, Saint-Denis, 38 p. 


Lang M, Alleaume S., Luque S., Baghdadi N., Féret J.-B., 2018. Monitoring and Characterizing 
Heterogeneous Mediterranean Landscapes with Continuous Textural Indices Based on VHSR Imagery. 
Remote Sensing, 10 (6), 868. 


Lefebvre A., 2011. Contribution de la texture pour l'analyse d'images à très haute résolution spatiale : 
application à la détection de changement en milieu périurbain, thèse de doctorat, Université Rennes 2, 
Université Européenne de Bretagne, 284 p. 


Maillard P., 2003. Comparing Texture Analysis Methods through Classification. Photogrammetric 
Engineering & Remote Sensing, 4 (11), 357-367. 


Pacifici F, Chini M., Emery W., 2009. A neural network approach using multi-scale textural metrics 
from very high-resolution panchromatic imagery for urban land-use classification. Remote Sensing of 
Environment, 113 (6), 1276-1292. 


Pandey G., Sharma V., Chaudhary P., Chowdary V., 2021. Integration of Texture and Spectral Response 
with AI Techniques for Buildings Footprint Identification Using High-Resolution Satellite Images. 
Journal of the Indian Society of Remote Sensing, 49 (6), 1439-1452. 


Proisy C., Couteron P., Fromard F., 2007. Predicting and mapping mangrove biomass from canopy grain 
analysis using Fourier-based textural ordination of IKONOS images. Remote Sensing of Environment, 
109 (3), 379-392. 


Puissant A., Hirsch J., Weber C., 2006. The utility of texture analysis to improve per-pixel classification 
for high to very high spatial resolution imagery. International Journal of Remote Sensing, 26 (4), 733-745. 


Regniers O., 2014. Méthodes danalyse de texture pour la cartographie doccupations du sol par télédé- 
tection très haute résolution : Application à la forêt, la vigne et les parcs ostréicoles, thèse de doctorat, 
Université de Bordeaux, 165 p. https://hal.archives-ouvertes.fr/tel-01143739. 


Ruiz Hernandez I, Shi W., 2018. A Random Forests classification method for urban land-use mapping 
integrating spatial metrics and texture analysis. International Journal of Remote Sensing, 39 (4), 
1175-1198. 


Teillet C., Pillot B., Catry T., Demagistri L., Lyszczarz D., Lang M, Couteron P., Barbier N., Adou 
Kouassi A., Gunther Q., Dessay N., 2021. Fast Unsupervised Multi-Scale Characterization of Urban 
Landscapes Based on Earth Observation Data. Remote Sensing, 13 (12), 2398. 

Tran A., Mangeas M., Demarchi M., Roux E., Degenne P., Haramboure M., Le Goff G., Damiens D., 
Gouagna L.C., Herbreteau V., Dehecq J.-S., 2020. Complementarity of empirical and process-based 
approaches to modelling mosquito population dynamics with Aedes albopictus as an example - 


Application to the development of an operational mapping tool of vector populations. PLoS One, 15 
(1), e0227407. 


83 


Partie 2 


Analyser et prédire l'effet 
de variables environnementales 
sur la distribution et la dynamique 
des moustiques vecteurs 


Un modèle est la représentation simplifiée d'un système réel et comporte donc 
toujours un grand nombre de simplifications afin de répondre à une question donnée. 
C'est pourquoi les modélisateurs citent souvent la phrase du statisticien George Box : 
« All models are wrong but some are useful» (tous les modèles sont faux, mais certains 
sont utiles). Concernant les maladies à transmission vectorielle, et en particulier celles 
impliquant un moustique, les approches de modélisation sont hautement complémen- 
taires des approches expérimentales et observationnelles. La modélisation peut être 
mobilisée dans deux perspectives. Il s’agit d'une part de mieux comprendre la distri- 
bution des vecteurs, leur dynamique, leurs liens avec des variables environnementales 
et climatiques (comme celles décrites dans la première partie de cet ouvrage), ainsi 
que le cycle épidémiologique de transmission des agents pathogènes vectorisés. Elle 
permet, d'autre part, de prédire, dans l’espace et le temps, la présence des moustiques 
vecteurs, leur abondance ou le risque de transmission, et les variations de ces éléments 
selon différents scénarios. Les approches de modélisation sont variées, le choix de 
l'approche dépendant de la question que l’on souhaite adresser, des connaissances que 
Yon a du système, des données disponibles, etc. 


Dans cette deuxième partie, différentes approches de modélisation sont présentées, ainsi 
que les outils opérationnels pour la surveillance entomologique et épidémiologique et 
la lutte antivectorielle développés à partir de ces approches, avec des exemples dans 
différents contextes géographiques. Le premier chapitre de cette partie (chapitre 6) 
concerne les modèles fondés sur les données (modèles de distribution d'espèces ou 
de niche écologique). Les chapitres suivants présentent des approches basées sur 
les connaissances (chapitre 7) et les processus à l'échelle des populations (modèles 
mécanistes : chapitre 8) ou des individus (modèles comportementaux : chapitre 9). 
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Chapitre 6 


Modèles basés sur les données : 
cartographier la distribution spatiale 
des vecteurs 


Yr MouA, EMMANUEL ROUX 


La cartographie d'une espèce quelconque (animale ou végétale), et a fortiori d'une espèce 
de moustique vectrice de maladie(s), peut s'appuyer sur les connaissances expertes et faire 
appel aux méthodes d'analyse multicritère qui seront présentées dans le chapitre suivant 
(chapitre 7) ou encore résulter de l'implémentation de modèles dits «mécanistes », c'est- 
à-dire explicitement basés sur les connaissances des processus bioécologiques régissant 
le cycle de vie des individus de l'espèce (cf. chapitres 8 et 9). Cependant, il se peut que 
les lacunes dans les connaissances relatives à l'espèce qui nous intéresse excluent de 
pouvoir construire de tels modèles. Une autre solution consiste alors à se fonder sur 
les observations disponibles et relatives à des sites spécifiques afin d’induire un modèle 
des habitats de l'espèce qui permettra de prédire la qualité de ces habitats (proxy de 
la probabilité de présence de l'espèce) sur l'ensemble de la zone étudiée. Ce chapitre 
traite de tels modèles, appelés «modèles de distribution d'espèces» ou «modèles de 
niche écologique ». Il décrit plus spécifiquement un des modèles les plus utilisés, appelé 
Maxent, ainsi que son application à la cartographie de l’habitat du principal moustique 
vecteur du paludisme en Amazonie, Anopheles darlingi. Les innovations relatives à la 
minimisation des impacts des biais d’échantillonnage sur les résultats de tels modèles, 
proposées par les auteurs de ce chapitre, sont également décrites. 


» Les modèles de distribution d'espèces 


Les modèles dits de distribution d'espèces ou de niche écologique sont construits à 
partir de données d'observation et visent à prédire spatialement la qualité d’habitat 
d’une ou de plusieurs espèces (animales ou végétales). Ces modèles relient les localisa- 
tions des occurrences connues des espèces considérées aux données environnementales 
décrivant le milieu dans lequel elles ont été observées (Guisan etThuiller, 2005). En 
sortie, une carte de qualité d’habitat est produite, pouvant être convertie en une carte 
de probabilité de présence de l'espèce ou des espèces étudiée(s). 


Depuis les années 1990, l’utilisation de ces modèles n’a pas cessé de s’amplifier. Ils sont 
utilisés pour différents objectifs : 
— mieux connaître la caractérisation des niches environnementales des espèces; 
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— prédire la distribution et les aires d'expansion des espèces invasives; 

— évaluer l'impact du changement climatique, de l’utilisation et de l'occupation du sol 
sur la distribution d'espèces; 

— prédire la distribution des espèces rares, en voie de disparition, notamment en 
soutien aux plans de conservation et de réintroduction d'espèces en danger. 


En ce qui concerne plus spécifiquement la prédiction des distributions de mous- 
tiques vecteurs, les Boosted Regression Trees, basés sur le principe de la modélisation 
d'ensemble en utilisant des arbres de régression, ont été utilisés pour cartographier la 
distribution des principaux vecteurs du paludisme sur le continent américain (Sinka 
et al., 2010) puis à l'échelle du globe (Sinka et al., 2012). 


Les variables environnementales utilisées en entrée informent sur les conditions 
environnementales de la zone étudiée. Elles sont appelées «variables explicatives », 
«prédicteurs », « covariables » ou «entrées». Elles sont souvent issues d'images satel- 
lite, de photos aériennes ou de données extrapolées a partir de données recueillies 
par des observations de terrain. Les conditions environnementales qu’elles caracté- 
risent peuvent exercer un effet direct ou indirect, et favorable ou non, sur la qualité de 
Vhabitat (cf. chapitre 1). 


De nombreuses applications des modèles de distribution d'espèces exploitent un 
ensemble de variables dites bioclimatiques, dérivées de variables météorologiques 
et censées mieux caractériser les conditions biologiques associées aux habitats des 
espèces. Une source de données largement utilisée est la base de données WorldClim!. 
Un tel jeu de variables est utile aux études à larges échelles ainsi que pour prévoir 
les effets à moyen et à long terme du changement climatique sur la distribution de 
l'espèce étudiée, mais ne rend pas compte explicitement de l'occupation et de l'usage 
du sol. Les données d'occupation et d’usage du sol à haute résolution spatiale et les 
variables paysagères qui en découlent permettent, à une échelle locale — faute de 
données à haute résolution spatiale disponibles à des échelles plus larges —, de mettre 
en évidence des relations complexes entre présence de l'espèce et paysage, ainsi que de 
considérer explicitement les activités anthropiques et leurs impacts. En revanche, elles 
permettent plus difficilement de réaliser des projections selon les différents scénarios 
de changements globaux. 


Concernant les données d'observation d'espèces, deux grandes approches se 
distinguent : la modélisation à partir des informations de présence et d'absence de 
l'espèce, pouvant s'appuyer sur toute approche statistique de discrimination ou 
d'apprentissage supervisée ; et la modélisation à partir des seules données de présence 
de l'espèce, supposant des approches spécifiques afin de pallier l'absence de « contre- 
exemples» dans les observations. La deuxième approche s’avére particulièrement 
intéressante, les informations d'absence étant souvent indisponibles ou peu fiables. 
En effet, l'absence d’une espèce en un lieu peut être due: 

— à la non-détection de l'espèce alors que l'espèce est effectivement présente; 

— à l'absence réelle de l'espèce dans un habitat favorable, pour des raisons liées à la 
dispersion: l'espèce n’a pas encore colonisé le milieu, mais y parviendra tôt ou tard; des 
barrières écologiques empêchent l'espèce d'accéder à certaines parties du territoire; 
— à l'absence réelle de l'espèce du fait de conditions écologiques non favorables. 


1. https://www.worldclim.org/data/index.html 
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Les données de présence peuvent quant à elles provenir de différentes sources : collec- 
tions et bases de données historiques d’universités, d’instituts, de musées; bases de 
données en ligne consacrées à la biodiversité (p. ex. Global Biodiversity Information 
Facility — GBIF”) ou à des espèces spécifiques, comme la base VectorBase? consacrée 
aux invertébrés vecteurs de pathogènes pour les humains; observations de terrain 
réalisées spécifiquement pour l'étude. 


Un des modèles les plus utilisés et qui utilise uniquement les données de présence est 
le modèle Maxent. 


» Le modèle Maxent 


Le modèle Maxent se fonde sur le principe du maximum d’entropie (Elith et al., 2011). 
En science de l'information, l’entropie est une mesure de quantité d’information et 
d'incertitude. Une entropie d’information élevée est associée à un contenu informa- 
tionnel varié et riche, plus à même de discriminer les différentes situations ou les 
différents objets sur lesquels portent les informations considérées. Maxent consiste 
à trouver la distribution de probabilité liée à la probabilité de présence de l'espèce 
considérée, sur l'ensemble du domaine spatial étudié (l'ensemble des pixels de la zone 
d'étude), dont lentropie est maximale et qui satisfait les contraintes imposées par 
les observations. Ces contraintes tiennent compte des conditions environnementales 
dans lesquelles l'espèce a effectivement été observée. Ces contraintes se traduisent 
mathématiquement, pour chaque variable environnementale, par l'égalité entre 
l'espérance mathématique de la variable — associée à la distribution de probabilité de 
présence de l'espèce — et la moyenne empirique des valeurs prises par la variable sur 
les sites où l'espèce a effectivement été observée (sites de présence). Il a été démontré 
que ce modèle fournissait des résultats robustes même avec une faible quantité 
d'observations (Hernandez et al., 2006). 


Lalgorithme Maxent est disponible dans le package demo (pour Species Distribution 
Modeling; Hijmans et al., 2017) de l'environnement de programmation Bi. 


» Biais d’échantillonnage et minimisation de leurs impacts 
sur la modélisation 


Toute analyse ou modélisation basée sur des observations sera biaisée si les obser- 
vations ont elles-mêmes été réalisées de façon biaisée et si aucune procédure nest 
mise en œuvre pour tenter, a minima, d'identifier les biais et, au mieux, de minimiser 
leur impact. 


La paramétrisation par défaut de Maxent fait intervenir un tirage aléatoire uniforme 
d'un grand nombre de sites dits de «background», supposés représenter l’ensemble 
des conditions environnementales de la zone d'étude (dans notre exemple, le territoire 
de la Guyane française) et leur fréquence relative. Une telle configuration suppose 
que l'effort d’échantillonnage (recherche de l'espèce sur le terrain) ait été réalisé de 


2. https://www.gbif.org/ 

3. https://vectorbase.org/vectorbase/app/ 

4. R Core Team, 2021. R: A language and environment for statistical computing. R Foundation for Statistical 
Computing, Vienna, Austria, https://www.R-project.org/. 
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façon équilibrée, c’est-à-dire que les différentes conditions environnementales exis- 
tantes aient bénéficié d’un effort d'échantillonnage comparable relativement à leur 
fréquence sur l’ensemble du territoire. Or une telle condition est rarement remplie, 
les relevés s’effectuant dans des lieux accessibles, le plus souvent à proximité des voies 
de communication, et en fonction d’autres critères, explicites ou non (lorsque l'espèce 
étudiée est un moustique vecteur, l’échantillonnage s'effectuera souvent en priorité 
autour des cas détectés de la maladie). 


Nous avons ainsi proposé une méthode de correction des biais d’échantillonnage dans 
les modèles exploitant les données de présence vs background, dont fait partie Maxent 
(Moua et al., 2020). La méthode repose, dans un premier temps, sur l'estimation de 
l'effort d’échantillonnage (dans notre cas en considérant les captures de toutes les 
espèces de moustiques qui ont été réalisées) dans l’espace des caractéristiques envi- 
ronnementales” ; puis sur la sélection aléatoire des sites de background dans ce même 
espace, pondérée par l'effort d’échantillonnage estimé précédemment. La méthode a 
été évaluée sur des données synthétiques et s'avère supérieure aux autres approches de 
débiaisage dans certains contextes d'application, notamment pour de faibles nombres 
de sites de présence. Lalgorithme d'estimation de l'effort d’échantillonnage est décrit 
en pseudo-code et évalué par Yi et al. (2020). 


» Application au principal vecteur du paludisme 
en Guyane francaise 


Le paludisme dans le monde et en Guyane francaise 


Le paludisme est une maladie causée par un parasite du genre Plasmodium, qui est 
transmis par certaines espèces de moustiques du genre Anopheles. En 2020, le nombre 
de cas de paludisme dans les 85 pays où la maladie est endémique (parmi lesquels la 
France — Guyane française) a été estimé à 241 millions. Le nombre de décès dans le 
monde a quant à lui été estimé à 627 000. Les bons résultats obtenus dans la lutte contre 
le paludisme entre 2000 et 2015 ont amené les Nations unies à inscrire l'élimination de 
la maladie dans au moins 35 pays d'ici 2030 (objectif 3.3 des Objectifs du développe- 
ment durable et Stratégie mondiale de lutte contre le paludisme). Cependant, depuis 
2015, la situation tend à se détériorer et le dernier rapport mondial sur la maladie fait 
état d’une augmentation de 12% des décès entre 2019 et 2020 (WHO, 2021). 


La Guyane française comptait environ 3500 cas par an en moyenne dans les 
années 2000. Ce nombre a drastiquement diminué à partir du début des années 2010 
pour se stabiliser, après 2013, autour de quelques centaines de cas par an. Toutefois, 
une recrudescence significative a été observée en 2017 et 2018 (Mosnier et al., 2020). 
Depuis 2019, le nombre de cas est de nouveau à la baisse et l'objectif d'élimination de 
la maladie sur le territoire national français a été réaffirmé, notamment en 2019 par 
Agnès Buzyn, alors ministre de la Santé et des Solidarités, lors de la 6° conférence 
de reconstitution des ressources du Fonds mondial contre les épidémies de VIH, de 
tuberculose et de paludisme. 


5. En d’autres termes, il s’agit de l'estimation de la densité de sites de collecte (ayant débouché ou non sur 
l'observation de l'espèce) dans l’espace défini par l'ensemble des variables décrivant les conditions environ- 
nementales de la zone d'étude, et non dans l’espace géographique (latitude, longitude). 


90 


Modèles basés sur les données : cartographier la distribution spatiale des vecteurs 


La stratégie mondiale de lutte contre le paludisme inclut notamment les actions de 
lutte antivectorielle et l’estimation du risque pour mieux cibler les actions de lutte 
et/ou anticiper d'éventuelles réintroductions dans les zones où la maladie a préalable- 
ment été éliminée. Cela repose en partie sur une meilleure connaissance des habitats 
favorables aux moustiques, de leur distribution spatiale et de la probabilité de présence 
de l'espèce ou des espèces vectrices. 


Anopheles darlingi, principal vecteur du paludisme en Guyane francaise 


Le moustique Anopheles darlingi se rencontre en milieu rural guyanais. Il pique de 
préférence l'être humain (comportement dit anthropophile), à l'extérieur ou à l’inté- 
rieur des maisons (exo-endophage), et préfère des lieux de repos situés à l'extérieur 
et végétalisés. Il privilégie des gîtes larvaires ensoleillés, mais suffisamment ombragés 
pour conserver une température entre 20 et 28°C (Hiwat et Bretas, 2011). Ses larves 
ont été retrouvées dans des étendues d’eau douce propres, avec peu ou pas de courant 
et de la végétation, telles que les berges de cours d’eau, les criques, les retenues d’eau 
formées à proximité des cours d’eau après inondation, les marécages, les savanes et 
les forêts inondées ou inondables (Hiwat et al., 2010; Rozendaal, 1987). Les gîtes 
larvaires se situent rarement en forêt dense, à cause de l'acidité de l’eau et du manque 
densoleillement sous le couvert forestier. 


An. darlingi est présent dans des environnements variés, de grande étendue, et en faible 
densité. Il est également extrêmement difficile à élever en laboratoire. Tout cela rend 
complexes sa capture, la cartographie de sa distribution spatiale et l’étude de sa bioéco- 
logie. Dans ce contexte, la modélisation de distribution d'espèces — et un modèle tel que 
Maxent en particulier — s'avère particulièrement utile. Le travail décrit ci-après a fait 
l'objet d’une publication dans la revue Journal of Medical Entomology (Moua et al., 2017). 


Données de présence d’Anopheles darlingi 


Une base de données historiques des observations de présence des moustiques 
Anopheles en Guyane a été construite à partir des archives de l’Institut Pasteur d’Algérie 
puis de l’Institut Pasteur de la Guyane, des rapports d'activités de la Direction de la 
démoustication et des actions sanitaires de la région Guyane puis de la collectivité 
territoriale de Guyane, du Service de santé des armées, des publications de l'Office 
de la recherche scientifique et technique outre-mer (Orstom, devenu Institut de 
recherche pour le développement, IRD, en 1998). Elle est disponible via la plateforme 
du Global Biodiversity Information Facility (Moua et al., 2019). 


Parmi ces données, seuls les sites de présence d'An. darlingi identifiés depuis 
l'année 2000 ont été pris en compte pour la modélisation, évitant ainsi de considérer 
des sites qui auraient pu subir des transformations significatives dans le temps et pour 
lesquels l’information de présence de l'espèce pourrait être remise en doute. Ces sites 
sont au nombre de 48, mais certains d’entre eux appartenant au même pixel de 1 km’, 
seuls 39 pixels constituent les données de présence en entrée de Maxent. 


La quasi-totalité des sites de présence avérée dAn. darlingi en Guyane française se 
situent sur la bande littorale et le long des deux fleuves frontaliers, faciles d’accés 
et où réside la grande majorité de la population. Cela suggère un important biais 
déchantillonnage et justifie l'application de la méthode décrite plus haut, visant à 
minimiser l’impact de tels biais sur les résultats de modélisation. 
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Données environnementales 


Le choix des données environnementales s’est fondé sur les connaissances issues de 
la littérature et sur l'expertise des entomologistes de l’Institut Pasteur de la Guyane. Il 
s'est avéré que trois grands types de milieux devaient être considérés : 

— les milieux naturels, pour lesquels la présence d'An. darlingi dépend de la valeur ou 
de la classe de la variable environnementale associée (MIL_NAT); 

— les milieux associés aux activités anthropiques qui altérent de maniére non perma- 
nente et localisée l'environnement naturel, influençant positivement la présence 
d'An, darlingi (ANTHROP_NON_PERM); 

— les milieux fortement et durablement anthropisés, correspondant à la présence et 
aux activités humaines qui altèrent de manière permanente l’environnement naturel 
sur de larges surfaces et qui limitent la présence d'An. darlingi (ANTHROP). 


Chacun de ces milieux a été caractérisé grâce à un jeu de variables (cf. tableau 6.1) dont 
la plupart sont issues de données de télédétection et dont les résolutions spatiales ont 
été homogénéisées à 1 km. 


Tableau 6.1. Données environnementales sélectionnées pour la modélisation des habitats 
favorables à Anopheles darlingi en Guyane française. Source : Moua et al. (2017). 


Variable Donnée source Type Type de milieu auquel 
et données de télédétection de variable se réfère principalement 
exploitées la variable 

Altitude Shuttle Radar Topography Numérique MIL_NAT 
Mission (SRTM), NASA continue 

(quantitative) 

Paysages géomor- Paysages géomorphologiques  Catégorielle MIL_NAT 

phologiques (Guitet et al., 2013) (qualitative) 

Unités géomor- dérivés du SRTM Catégorielle MIL_NAT 

phologiques (qualitative) 

Occupation Gond et al. (2011), Catégorielle MIL_NAT 

du sol exploitant les données (qualitative) 

SPOT4/VÉGÉTATION 

Présence et Empreinte humaine Catégorielle ANTHROP_NON_PERM 

activités humaines (Thoisy de er al., 2010) ordonnée 

qui influent (semi- 

de maniére non quantitative) 

permanente et 

locale sur le milieu 

(HUMACT) 

Densité des BD TOPO”, IGN Numérique ANTHROP_ NON DERM 

routes et pistes continue 

(ROUTES) (quantitative) 

Pourcentage Empreinte humaine Numérique ANTHROP 

d'urbanisation (Thoisy de et al., 2010) continue 

dans le voisinage (quantitative) 

(URBAN) 


L'ensemble de la méthode d’implémentation de Maxent pour la cartographie des habi- 
tats d’An. darlingi en Guyane française est schématisé figure 6.1. 
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Résultats 


La figure 6.2 représente la qualité d'habitat pour An. darlingi. La probabilité de 
présence de l'espèce est liée de façon monotone à la qualité d'habitat (plus la qualité 
est «bonne», plus il y a de chances que l'espèce soit présente. 


Les résultats montrent qu'une qualité d’habitat élevée est significativement liée à 
la présence et aux activités humaines, caractérisées par les variables HUMACT et 
ROUTES, en particulier sur le littoral, au centre et à l’ouest du territoire, même si dans 
ces deux dernières zones les voies de communication sont essentiellement les rivières 
et les pistes non asphaltées. 


La contribution significative de la longueur de routes et de pistes (ROUTES) à la 
qualité d'habitat apparaît cohérente avec les connaissances sur l'impact des modi- 
fications environnementales et de l’activité humaine : la construction des routes 
et des pistes est accompagnée de trouées dans la forét, générant potentiellement 
de nouveaux gites larvaires (fossés p. ex.) [Singer et Castro de, 2001] et apportant 
Yensoleillement nécessaire au développement des larves ainsi que la présence d'hôtes 
humains. Il est cependant à noter qu’au-dela d'un certain seuil, autour de 7km de 
routes/pistes par km’, la qualité d’habitat n’augmente plus puis diminue, signifiant 
que la densification du réseau routier (caractéristique d'une urbanisation de la 
zone) provoque des changements des milieux défavorables à An. darlingi (absence 
de végétation pour le repos des femelles gravides, pollution des eaux de ruissel- 
lement, etc.). Pour des raisons similaires, le phénomène d'urbanisation (caractérisé 


Saint-Laurent- 
du-Maroni , 


` BRÉSIL 


Indice de qualité d’habitat 
E Très élevé 

E Elevé 

B Moyen 

E Faible 

O Très faible 


25 50 75 100 km 


RGFG95/UTM22 
Unité : métre 


Figure 6.2. Qualité d’habitat d’Anopheles darlingi en Guyane française. 
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par la variable URBAN) n'est pas favorable à An. darlingi (Stefani et al., 2013). Cela 
sobserve particulièrement dans les zones de Cayenne et de Kourou. En revanche, 
Saint-Laurent-du-Maroni, une des trois plus grandes zones urbaines de Guyane 
par la taille et la densité d’urbanisation, avec Cayenne et Kourou, apparait avec une 
qualité d’habitat très élevée du fait notamment d’une densité de bâti estimée plus 
faible. Cependant, Saint-Laurent-du-Maroni connaissant une croissance démo- 
graphique soutenue (+2,3% par an entre 2013 et 2019 d’après l’Institut national 
de la statistique et des études économiques, INSEES), une actualisation de la carte 
de distribution des habitats favorables au moustique vecteur du paludisme tenant 
compte des récentes expansion et densification du bâti s'avère nécessaire. 


Les activités humaines (variable HUMACT) liées à des habitats particulièrement 
favorables à An. darlingi dans le centre du territoire sont essentiellement associées à 
l'activité minière. En 2014, l’orpaillage illégal en Guyane française représentait environ 
700 sites où travaillaient 10000 à 15000 personnes (Douine et al., 2016). Il contribue 
de différentes manières au maintien du paludisme dans la région. En effet, l’activité 
dorpaillage favorise la prolifération d'An. darlingi en créant des trouées en forêt dense 
et des retenues d'eau mi-ombragées qui sont autant de gîtes larvaires potentiels. La 
fréquentation des sites d’orpaillage par une population humaine nombreuse et forte- 
ment exposée aux piqûres de moustiques, par ses activités en extérieur et des conditions 
d'habitation précaires, génère des situations à très haut risque de transmission locale 
du paludisme. De plus, cette population, particulièrement mobile et peu suivie médi- 
calement, contribue fortement à la circulation du parasite dans les populations de la 
région, favorisant l'émergence et le maintien de foyers de transmission, de résistances 
aux traitements contre le paludisme, etc. (Thoisy de et al., 2021; Douine et al., 2021). 


Il est également à noter que An. darlingi rest pas le seul vecteur du paludisme présent 
sur les sites d’orpaillage. En particulier, An. marajoara semble contribuer significative- 
ment à la transmission dans des sites forestiers très anthropisés comme les sites miniers 
(Pommier de Santi et al., 2016). L'ensemble de ces éléments désignent lorpaillage, et 
plus spécifiquement lorpaillage illégal, comme un des principaux facteurs permettant 
le maintien d’une transmission autochtone du paludisme et comme un des principaux 
freins à l'élimination de la maladie dans la région. 


Certaines classes d’occupation du sol sont associées 4 une qualité d’habitat élevée. 
C'est le cas de la classe « Savane boisée / forêt sèche», qui correspond à des zones 
sèches mais régulièrement inondées et donc associées à la présence de potentiels gîtes 
larvaires (Gond et al., 2011; Rosa-Freitas et al., 2007). Un tel résultat est cohérent 
avec des captures effectuées dans les savanes le long du littoral (Dusfour et al., 2013; 
Vezenegho et al., 2015) ayant confirmé la présence d'An. darlingi dans ce type de 
milieu, en abondance parfois élevée. 


» Conclusion 


Dans ce chapitre, nous avons vu qu'il est possible de construire, à partir des seules 
informations de présence de l'espèce et des conditions environnementales sur les sites 
d'observation, des modèles permettant de prédire la qualité d'habitat et par exten- 
sion la probabilité de présence de l'espèce en tout point de la zone d'étude Parmi ces 


6. https://www.insee.fr/fr/statistiques/6012651 
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modèles, Maxent est un des plus populaires, notamment par sa relative robustesse 
lorsque le jeu de données de présence est de taille réduite. Comme tout modèle, 
il est cependant sensible aux biais d’échantillonnage, a fortiori lorsque le nombre de 
sites d'observation est faible, mais des approches variées, dont une proposée par les 
auteurs de ce chapitre, permettent de minimiser l'impact de ces biais. La tendance en 
modélisation de distribution d'espèces est également à la modélisation ensembliste, 
consistant à appliquer plusieurs modèles de distribution d'espèces puis à combiner les 
sorties des modèles donnant les meilleurs résultats. Une librairie R a été développée 
pour ce type d'approche ` biomod?2 (Thuiller et al., 2022). Cependant, toute démarche 
de modélisation requiert de comprendre la paramétrisation du modèle et de l'adapter 
au contexte et aux objectifs de l'étude afin d'assurer l'obtention de résultats fiables et 
pertinents. L'application de plusieurs modèles peut requérir un effort conséquent. 


In fine, la cartographie de la distribution spatiale de l'espèce vectrice est censée 
guider les actions de lutte antivectorielle. Elle est cependant peu exploitée dans 
ce cadre, faute d’une résolution spatiale suffisamment fine et aussi parce qu'elle ne 
prédit pas l'abondance des vecteurs. Des modèles de prédiction de l'abondance de 
moustiques construits à partir de données ont été proposés pour pallier ce manque. 
Tran et al. (2020) présentent par exemple un modèle basé sur l'approche d’appren- 
tissage par machines à vecteurs de support (SVM) à partir de relevés des densités 
de larves de moustiques Aedes dans différentes localités de l’île de La Réunion. Ce 
modèle produit des résultats prometteurs et complémentaires à ceux obtenus par 
modélisation mécaniste (cf. chapitre 8). Adde et al. (2016) ont, eux, modélisé avec 
succès, par régression logistique ordinale à effets mixtes, la dynamique de la densité 
dAn darlingi dans une localité de Guyane française. Cependant, de tels modèles 
sont difficilement applicables sur de vastes zones. Ils sont en effet particulièrement 
gourmands en données et les données d’abondance de moustiques sont particuliè- 
rement difficiles et coûteuses (en temps, personnel, argent) à obtenir, donc rares. De 
plus, il est difficilement envisageable de les récolter sur de vastes territoires et à des 
fréquences temporelles élevées (à noter cependant que les méthodes de piégeage de 
moustiques se diversifient et deviennent de plus en plus efficaces et bon marché, ce 
qui permet d'envisager de développer des modèles qui soient de plus en plus robustes 
et généraux). Enfin, une des raisons de la faible appropriation des modèles de distri- 
bution d'espèces pour la lutte contre les maladies vectorielles est que la distribution 
des vecteurs ne peut s’interpréter directement en termes de risque de transmission 
et encore moins de risque épidémique. Ces derniers font intervenir de nombreux 
autres facteurs liés au vecteur (abondance et âge des populations), mais aussi liés à 
la présence et aux activités humaines, au système de santé, à des facteurs individuels 
humains (comportement, susceptibilité génétique, immunité, etc.). 


Pour une appropriation effective de la modélisation de la distribution d'espèces dans la 
lutte antivectorielle et plus généralement pour la planification des actions de contrôle 
des maladies vectorielles, il conviendrait donc de mieux spécifier les apports de 
chaque approche de modélisation. On pourrait ainsi aller vers des approches plus inté- 
grées, systémiques, permettant de considérer différentes composantes du risque en 
exploitant différentes approches et formalismes de modélisation (à partir des connais- 
sances, mécanistes, à partir des données) et leurs avantages respectifs, en fonction des 
données et des connaissances disponibles. 
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Chapitre 7 


Modèles fondés sur les connaissances : 
exemple d’un outil d'évaluation multicritère 
pour la santé publique 


FANJASOA RAKOTOMANANA, 
HOBINIAINA ANTHONIO RAKOTOARISON 


Pour réduire l'impact d'une maladie pour les individus, la prise en charge adéquate des 
cas est fondamentale. Cependant, les actions de santé publique humaine et animale 
visant à prévenir les risques de diffusion des maladies sont indispensables pour réduire 
leur impact à l'échelle de la population. Les maladies infectieuses affectant les humains 
et les animaux surviennent à la suite de leur exposition à des agents infectieux 
présents dans leur environnement direct ou indirect. Ces agents infectieux, ainsi que 
leurs arthropodes vecteurs (organismes vivants transmettant l'agent infectieux), sont 
eux-mêmes sous l'influence de déterminants environnementaux dont les paramètres 
varient en fonction des conditions physiques, chimiques, climatiques et des pressions 
anthropiques exercées sur leur habitat naturel. La compréhension de ces déterminants 
dans l'épidémiologie de ces maladies et la surveillance de leur évolution sont donc 
indispensables pour détecter des changements dans l’environnement et appréhender 
les risques d’émergence d'une épidémie. 


Lors de la survenue d’épidémies ou de pandémies, surtout avec l'émergence d'un 
agent pathogène inconnu, le recours à la modélisation est devenu une pratique 
courante en santé publique. La modélisation descriptive d'un phénomène de santé 
est habituelle (elle permet notamment d'estimer l'incidence et la prévalence de la 
maladie), mais la modélisation peut également viser à comprendre les mécanismes 
de transmission ou à prédire l'ampleur d'une épidémie ou l'impact de mesures de 
lutte (Dubois, 2005). Les modèles apportent des informations permettant d'estimer 
les moyens techniques et financiers nécessaires pour faire face à une éventuelle situa- 
tion critique en santé publique. Ils jouent un rôle essentiel dans l’aide à la décision. 
L'introduction de la dimension spatiale dans un modèle est aussi essentielle pour 
la compréhension du phénomène d'’initialisation d’une épidémie. La plus ancienne 
cartographie thématique de la distribution des cas de choléra de Londres en 1854 a 
montré que l'intégration des informations spatiales dans une démarche d'épidémio- 
logie explicative a permis de comprendre le mécanisme de diffusion de la maladie et 
de contrôler l'épidémie (Snow, 1855). 
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Un grand nombre de modèles reposent sur l'analyse de données (cf. chapitre 6). D’ail- 
leurs, des initiatives et des plateformes de collecte et de gestion de données existent à 
l'échelle mondiale ou continentale, avec création de tableaux de bord de surveillance 
des maladies à potentiel épidémique. Cependant, de nombreuses contraintes pèsent 
sur l'acquisition de données robustes, complètes, temporellement et spatialement en 
adéquation avec la problématique. Par exemple, depuis 2019 le ministère de la Santé 
publique de Madagascar s'appuie sur le Système de gestion d’information sanitaire 
(SGIS) proposé par la plateforme DHIS2!, utilisée par plus de 100 pays dans le monde. 
En 2021, les taux de complétude et de promptitude de rapportage des données restent 
inférieurs à 80% pour le système national de surveillance intégrée opérationnel du 
paludisme, ce qui constitue une des faiblesses dans la gestion et mise en œuvre des 
réponses aux épidémies de paludisme (PNLP, 2021). De plus, les données concernant 
les facteurs de risques environnementaux font également souvent défaut dans les pays 
à ressources limitées comme Madagascar, mais aussi globalement lors de l'émergence 
de nouveaux pathogènes. Pourtant ces données sont indispensables dans les stratégies 
de lutte contre les maladies à potentiel épidémique comme le paludisme (WHO, 2001). 


Ce chapitre traite des modèles fondés sur les connaissances, qui peuvent être mis en 
œuvre en l'absence de données épidémiologiques, avec l'exemple de l'analyse multi- 
critère spatialisée. Cette méthode aide à contourner les besoins en importants jeux 
de données épidémiologiques, nécessaires pour développer des modèles basés sur 
les données (cf. chapitre 6). Ce chapitre montre un exemple de modèle fondé sur les 
connaissances utilisé dans la détermination des zones à risques prioritaires pour la 
campagne de lutte antivectorielle contre le paludisme à Madagascar. Enfin, ce chapitre 
illustre l'outil Full MCE for public Health développé à l'Institut Pasteur de Madagascar 
(IPM), qui permet d'appliquer cette méthode par une extension (plug-in) pour le 
logiciel de système d’information géographique (SIG) QGIS. 


» L'approche multicritère spatialisée, une approche fondée 
sur les connaissances 


Les logiciels de SIG sont adaptés pour gérer des données environnementales multi- 
sources et offrir un outil d’aide a la décision. Lapproche d’analyse multicritére (AMC) 
combinée au SIG est de plus en plus adoptée dans différents domaines. Elle est fondée 
sur les connaissances des experts. Quelques exemples ont montré que ‘AMC combinée 
au SIG, également appelée analyse multicritère spatialisée (AMCS), est très puissante 
pour répondre aux problèmes liés à la santé humaine ou animale. Par exemple, Bell 
et al. (2007) ont construit des indices de privation, qui sont des indicateurs liés à la 
santé des populations, en proposant aux décideurs des gradients socio-économiques 
basés sur une approche AMCS. En l'absence de données précises sur les foyers de 
peste porcine africaine (PPA), lAMCS a permis de cartographier les risques de trans- 
mission de la PPA en Afrique et à Madagascar, ce qui permet de compenser le manque 
d'études quantitatives sur l’épidémiologie de la maladie en Afrique (Glanville de et al., 
2014). La modélisation fondée sur les connaissances a également été utilisée pour la 
détermination des zones à risques d'amplification et de propagation de la fièvre de 
la vallée du Rift (FVR) dans quatre pays de l'Afrique de l'Est : le Kenya, la Tanzanie, 
l'Ouganda et l’Éthiopie (Tran et al., 2016). 


1. https://dhis2.org/fr/about/ 
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LAMCS est constituée de cing principales étapes, considérées comme des éléments 
constitutifs des procédures spatiales d'aide à la décision : 


1. La définition du risque à cartographier (risque d'introduction, d'amplification, de 
diffusion, etc.) et du territoire concerné selon l'objectif nécessitant un (des) choix pour 
une prise de décision. La bonne connaissance des risques aide à identifier les facteurs 
à prendre en compte dans l'étape suivante. 


2. L'identification des critères pour l'analyse : les données qui vont servir comme 
facteurs sont identifiées avec lavis des experts et collectées à partir de différentes 
sources (instituts nationaux de géographie ou autres portails d’information géogra- 
phique). Les critères peuvent se décliner en deux types : les contraintes (un masque qui 
limite les zones à prendre ou ne pas prendre en considération) et les facteurs (critères 
définissant un certain degré d’aptitude pour une zone donnée) [Rakotoarison et al., 
2020; Rakotomanana et al., 2007]. Les facteurs (p. ex. la température, qui impacte la 
dynamique de population des moustiques et le cycle biologique du plasmodium chez 
les moustiques) et les contraintes (p. ex. les zones habitées, cibles des interventions 
de lutte antivectorielle : il n’est pas utile de prédire le risque en dehors de ces zones) 
sont donc traités et combinés de manière différente. Les experts comme les épidémio- 
logistes, entomologistes, climatologues, écologues ou spécialistes en santé humaine et 
animale interviennent à ce stade d'identification des critères. Leurs avis sont appuyés 
par la revue de littérature résultant d'une recherche bibliographique préalable. Les 
données de télédétection constituent une source particulièrement utile pour acquérir 
des informations spatialisées sur des facteurs ou contraintes climatiques et environ- 
nementales. L'utilisation des données satellitaires pour appréhender les différents types 
d'occupation du sol et leur évolution dans le temps est devenue l'application la plus utile 
dans un SIG pour disposer de données récentes sur une zone d'étude (cf. chapitre 1). 


3. La standardisation des facteurs : elle consiste à transformer les facteurs d'évaluation 
originaux en unités comparables afin d'obtenir des aptitudes à un événement (p. ex. 
aptitude à la survenue d'une maladie) sur une échelle continue (0 à 1 ou 0 à 255 selon 
le choix). Une valeur de zéro correspond à une aptitude nulle tandis qu’une valeur 
de 1 ou 255 (selon le choix de l'échelle utilisée) représente une aptitude maximale à 
l'événement étudié. La standardisation peut se faire en utilisant des seuils (au-delà 
ou en deçà d’un seuil, le risque est présent ou absent) ou en s'appuyant sur la logique 
floue, approche la plus commune en analyse multicritère pour l’aide à la décision. 
Avec cette technique, la transition entre l'aptitude ou non d’un facteur à engendrer un 
risque se fait de manière progressive d’une zone à une autre, selon la relation entre le 
facteur et le risque (croissante, décroissante, symétrique, etc.). La figure 7.1 montre 
deux exemples de fonctions de relation entre les facteurs et les risques à évaluer. 


4. La pondération des facteurs : cette étape traite de l'attribution d’un poids (entre 0 
et 1) à chaque facteur selon son importance par rapport aux autres facteurs consi- 
dérés. Plus le poids est proche de 1, plus la contribution du facteur à l'événement 
est importante. Plusieurs méthodes, dont la comparaison par paires, ont été suggé- 
rées pour pondérer les facteurs selon les contextes. La méthode de comparaison par 
paires a été développée par Saaty dans le cadre d’un processus d'analyse hiérarchique 
(AHP pour Analytic Hierarchy Process) [Saaty, 1977]. Elle repose sur la construction 
d'une matrice où les facteurs sont comparés par paires en attribuant des notes à partir 
d’une échelle de comparaison (tableau 7.1). À partir de cette matrice, la méthode AHP 
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a Abondance b Abondance 
des moustiques des moustiques 
b,c,d 
1 1 b c 
0 > 0 > e 
a Surface en eau a d Température 


Figure 7.1. Exemples de fonctions d'appartenance pour la standardisation des facteurs : (a) 
fonction croissante; (b) fonction symétrique avec plateau. 


permet de calculer le poids de chaque facteur. La comparaison par paires peut, néan- 
moins, présenter des relations incohérentes. Ces incohérences peuvent être mesurées 
par le calcul du ratio de cohérence RC. Des comparaisons par paires jugées cohérentes 
ont un RC < 0,1. Une valeur de RC = 0,1 indique des jugements incohérents et une 
comparaison qui devra être réévaluée. 


Tableau 7.1. Exemple d'une échelle de comparaison de deux facteurs de risque. 


Importance du facteur A par rapport au facteur B 


5 z z S 
2 2 t S g S = 2 2 
£ n Ÿ $ £ g £ 2 n l £ 
vo v E E oO E oO E © E vo 
£ ZG © g 5 us 3 T =g £ 
D Ka E = © S Ka E Aes D 
z S £ -o a £ S š 
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A 
Moins important Plus important 
1/9 1/7 1/5 1/3 1 3 5 7 9 


5. La dernière étape consiste en la combinaison des indices de risque (facteurs standar- 
disés et pondérés) et en leur intégration dans un système d’information géographique 
produisant une carte de risque. La carte obtenue traduit le niveau de risque par unité 
géographique élémentaire. Avec un choix d'échelle continue allant de 0 à 255, le risque 
faible est traduit par la valeur 0 tandis que la valeur 255 est interprétée comme un risque 
élevé. La variation des poids selon l'estimation et lavis des experts permet d'estimer 
l'incertitude au niveau d’une unité géographique donnée. Une carte d'incertitude peut 
être associée à la carte de risque. Le modèle est considéré comme stable s’il ne varie 
pas significativement, quelle que soit la variation de poids attribuée aux facteurs par 
les différents experts. Une validation de la carte de risque peut être effectuée selon la 
disponibilité des données épidémiologiques. 


» Exemple d'application sur le risque de paludisme 
dans les Hautes Terres centrales de Madagascar 


Contexte 


La lutte contre le paludisme demeure un défi mondial pour la santé publique, mais 
également pour l'économie. L'objectif de la lutte contre le paludisme à l'échéance 2030 
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est de réduire de 90% l'incidence et la mortalité liées au paludisme dans le monde, par 
rapport aux données de références de 2015. En 2019, le nombre de cas de paludisme a 
été de 227 millions dans 85 pays où la maladie est endémique; 14 millions de plus ont 
été estimés en 2020, soit 241 millions (WHO, 2021). Pour les pays à ressources limi- 
tées, la lutte contre ce fléau dépend principalement des partenaires financiers à travers 
le fonds mondial destiné aux grandes épidémies, telles que le paludisme, le sida et la 
tuberculose (PNLP, 2017; WHO/GMP, 2021). À Madagascar, le paludisme constitue 
la quatrième cause de morbidité dans les centres de santé et de mortalité hospitalière 
(PNLP, 2017). Le paludisme est caractérisé par quatre faciès selon la durée et l'intensité 
de la transmission : Est, Ouest, subdésertique et Hautes Terres centrales (HTC). 


Application de la méthode AMCS à l'aide de l'outil Full MCE for Public Health 


Rakotoarison et al. (2020) ont appliqué la méthode AMCS sur les HTC afin de déter- 
miner les zones à risques prioritaires pour les campagnes d’aspersion intra-domiciliaire 
dinsecticide (AID). Le modèle a été développé avec l'outil Full MCE for Public Health. 


En effet, pour des utilisateurs non spécialistes en géomatique, l’'IPM a créé une exten- 
sion au logiciel libre QGIS pour permettre aux acteurs de santé de s'approprier un 
outil interactif et semi-automatique. L'extension Full MCE for public Health a été 
développée pour QGIS version 2.x avec le langage Python 2 et une interface convi- 
viale? (figure 7.2). Des formations à son utilisation ont été dispensées par l'équipe 
de l’'IPM à l'échelle nationale et internationale pour faciliter la prise en main par les 
personnels de santé humaine ou animale. 


Les étapes de la méthode AMCS décrites dans le paragraphe précédent sont détaillées 
ci-après, avec des illustrations de l'interface de l'outil Full MCE for Public Health. 


if Full MCE for Public Health v1.0 ? xX 


Plugin d'Analyse Multicirtère pour la santé publique Version Desktop 
Fonctions Développeur 
Redassification des contraintes Mampionona RASAMIMALALA 
Normalisation des facteurs 


Pondération par la méthode de comparaison par paires 


Agrégation par la méthode de combinaison linéaire pondérée 


Concepteurs miciaire (CAID) 


domaine de 


Mampionona RASAMIMALALA (rasamimalala @gmail.com) 
Anthonio RAKOTOARISON (anthonio@pasteur mg) 
Fanjasoa RAKOTOMANANA (fanja@pasteur.mg) 


Patrice PIOLA (ppiola @pasteur.mg) 


[Précédent | | suivent > Executer Annuler 


Figure 7.2. Interface d'accueil de l'extension Full MCE for public health dans QGIS 2.x. 


1. Définition des risques et des zones prioritaires 


En premier lieu, les HTC ont été identifiées comme zone prioritaire, car il s’agit d’une 
zone à paludisme instable, à risque épidémique. La stratégie de lutte antivectorielle 
contre le paludisme y est basée sur les campagnes d’AID. Le choix des zones d'inter- 
vention dépend de la disponibilité des moyens de lutte et la cartographie des zones les 
plus à risque doit permettre de mieux cibler les zones prioritaires. 


2. https://github.com/SaGEOTeam/FullMCE 
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2. Identification des critères 


Les différents critères qui interviennent dans le risque de transmission du paludisme 
dans les HTC ont été identifiés dans le cadre d’une collaboration avec des spécialistes 
en entomologie médicale, des acteurs de santé au niveau du programme national de 
lutte contre le paludisme (PNLP) et d’autres experts locaux. Cette étape repose sur les 
connaissances sur le paludisme en général et les HTC de Madagascar en particulier. 


Le chapitre 1 a traité du rôle de l’environnement dans la transmission des mala- 
dies dues aux moustiques et de l'intérêt de l’utilisation de la télédétection. En ce 
qui concerne le paludisme sur les HTC, les rizières sont les gîtes larvaires de prédi- 
lection du moustique Anopheles funestus, vecteur principal dans cette région. Pour 
différentes raisons (variations des précipitations, possibilité d'irrigation, insécurité), 
les surfaces cultivées peuvent varier d’une année à l’autre. Pendant une période de 
repos d’une rizière, elle peut rester humide ou sèche selon le système d'irrigation. 
Parfois les parcelles changent complètement de vocation en devenant des lieux de 
prélèvement de terre pour la fabrication de briques. Les trous humides laissés par 
ces activités anthropiques constituent également des gîtes larvaires idéaux pour les 
moustiques vecteurs du paludisme. Aussi, quand les surfaces rizicoles diminuent, 
on assiste au développement, dans les espaces naturels, d’autres espèces vectrices, 
telles ou An. gambiae. Les zones humides, constituant de potentiels gîtes larvaires 
des anophèles vecteurs du paludisme, sont donc un des principaux critères retenus 
pour l'AMCS. La télédétection (images Landsat) permet de cartographier la super- 
ficie rizicole ainsi que celle des zones humides naturelles à une fréquence annuelle 
(Rakotoarison et al., 2020). 


Laltitude, la température, la précipitation, la densité de la population humaine (qui 
peut avoir un impact sur la distribution des gites larvaires), les zones habitées (qui 
sont les cibles des AID) ont également été retenues comme critéres (figure 7.3) 
(Rakotoarison et al., 2020). 


Données environnementales Données démographiques Données climatiques 
Images Densité Oe z EE 
SRTM Sadir del R [Zones habitées | | Température | | Précipitation | 
Extraction Classification 
Attitude Occupation 
du sol 


Standardisation 


Pondération 
Agrégation 


Carte de risque = 


Figure 7.3. Données en entrée et flux de traitement des données pour l'obtention des cartes de 
risque de paludisme sur les Hautes Terres centrales de Madagascar. Adapté de Rakotoarison et 
al. (2020). 
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Des prétraitements (p. ex. classification d'Images satellite) sont nécessaires pour 
préparer les données correspondantes pour l'analyse (figure 7.3). L'interface du plug-in 
Full MCE for Public Health permet d'intégrer ces différents critères comme contraintes 
ou facteurs, au format raster (figure 7.4). 


a # Full MCE for Public Health v1.0 ? X 

Données d'entrée: CONTRAINTES 

Nombre de contraintes | 2 gl 

Noms Chemins Prêtes 
1|DensPop D:/Traitement/DensPop15Htc.rst 
2|ZoneHab D:/Traitement/ZHabHtcBool.rst Ar v 
< Précédent Suivant > Executer Annuler 

b # Full MCE for Public Health v1.0 ? x 


Données d'entrée: FACTEURS 


Nombre de facteurs [ 3h 


Noms Chemins Normalisés 
1/Altitude D:/Traitement/AltHtc.rst 
2|DensPop D:/Traitement/DensPop15Htc.rst 
3 [Precipitation D:/Traitement/PrecipHtc.rst 
< Précédent Suivant > Executer Annuler 


Figure 7.4. Illustration de l'interface du plug-in Full MCE for Public Health correspondant à la 
prise en compte (a) des contraintes et (b) des facteurs à considérer dans l'analyse multicritère 
spatialisée. 


3. Standardisation 


L'interface du plug-in permet ensuite de standardiser les facteurs. Selon les choix des 
experts, la forme des fonctions de standardisation (linéaire, sigmoide, etc.) et leur sens 
de variation (croissant, décroissant, symétrique) peuvent être paramétrés, ainsi que les 
seuils (figure 7.5a). 


4. Pondération des facteurs 


Dans l'étape de pondération, l'utilisateur doit comparer les différents facteurs deux à 
deux. La saisie des valeurs de la matrice de comparaison par paires (tableau 7.1) est 
facilitée par l'interface du plug-in (figure 7.5b). Le poids attribué à chaque facteur est 
calculé automatiquement selon la méthode AHP et résumé dans l'interface de calcul 
final (figure 7.6). 
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Dans l'étude de Rakotoarison et al. (2020), la densité de population a été considérée 
comme le facteur le plus important par les experts, avec une valeur de pondération 
calculée selon la méthode AHP la plus élevée (poids : 0,50), suivi de la distance par 
rapport aux rizières et zones humides (plus la distance est importante, moins le risque 
est élevé; poids : 0,18), suivi de la température (poids : 0,17), de l'altitude (poids : 0,09) 
et de la précipitation (poids : 0,06). 


a 


# Full MCE for Public Health v1.0 ? x 
Paramétres de normalisation des facteurs 
Nombre de facteurs 3 + 
Noms Fonctions Sens A B SS D 

1 [Altitude S-Shaped (Sigmt > (Croissant | 500 1500 

2 |DensPop J-Shaped > |Décroisse ~ 120 300 

3 [Precipitation Linéaire ~ Symétriq + 120 140 160 200 

< Précédent Suivant > Executer Annuler 

b # Full MCE for Public Health v1.0 ? x 


Pondération des facteurs - Méthode de comparaison par paires 


Matrice de jugement 


Altitude DensPop Precipitation | 


Altitude 1 0.14286 1 
DensPop 7 1 3 Tester 
Precipitation 1 0.33333 1 
Enregistrer 
Importer 
RC = 0.07 RC < 0.1 => Matrice de jugement cohérente et ac 
< Précédent Suivant > Executer Annuler 


Figure 7.5. Illustration de l'interface du plug-in Full MCE for Public Health correspondant aux 
étapes (a) de standardisation et (b) de création de la matrice de comparaison par paires pour la 
pondération des facteurs. 


5. Combinaison des indices de risque 


Les facteurs (densité de population, distance aux zones humides, température, altitude, 
précipitation), standardisés et pondérés, ainsi que les contraintes (zones habitées, alti- 
tude, densité de population) sont finalement combinés pour obtenir les cartes des zones 
à risques prioritaires pour AID pour différentes années (2014, 2015 et 2016) par la 
mise à jour annuelle des cartes d'occupation du sol, de température et de précipitations. 
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Les zones à risques identifiées ont été classées selon un gradient allant de très faible à 
très élevé et mettent en évidence des variations spatiales et interannuelles de risque de 
transmission importantes (figure 7.6). 


Gradient de risque e? Faible E Élevé 0 25 50km 
= 
Limite des districts [J Très faible L Modéré ` J Tres élevé 


Figure 7.6. Cartes de risque de paludisme dans les Hautes Terres centrales malgaches. Adapté 
de Rakotoarison et al. (2020). 


» Conclusion 


L'influence de la localisation géographique et le rôle de l'environnement sur la distri- 
bution des maladies sont des concepts très anciens dans le domaine de la santé. 
Lutilisation des données satellites pour appréhender les différents types d'occupation 
du sol et leur évolution dans le temps est cruciale pour disposer de données récentes 
et mettre à jour les cartes de risque. Le recours aux modèles basés sur les connais- 
sances favorise l'approche participative et l'implication des experts dans l'appréciation 
des risques. L'outil développé à partir de cette méthode est libre d'utilisation. Il peut 
être utilisé dans d’autres domaines, tels l'agriculture, la protection de l’environnement 
proprement dite ou l'aménagement du territoire. 
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Chapitre 8 


Arbocarto : un modèle mécaniste fondé 
sur le cycle de vie des moustiques Aedes 


RENAUD MARTI, MARIE DEMARCHI, 
MATHIEU CASTETS, ANNELISE TRAN 


Par la nature de leur construction, les modèles basés sur les processus, désignés 
également comme «modèles mécanistes», sont des modèles explicatifs, centrés 
sur la causalité des relations entre les entrées et les sorties (Craver, 2006). Ce type 
d'approche nécessite d’expliciter ces relations (p. ex. l'effet de la température sur le 
développement des moustiques) sous la forme de liens de causalité énoncés dans un 
espace simplifié, le plus souvent sous une forme mathématique (p. ex. une équation 
exprimant le taux de développement comme une fonction de la température), en 
se basant sur des connaissances déjà établies (par des études observationnelles ou 
expérimentales) sur le système modélisé. 


Selon les objectifs de modélisation, cette approche implique aussi une nécessaire et 
parfois difficile démarche de simplification. Cette démarche, à l'initiative du modé- 
lisateur et mobilisant ses connaissances sur l’importance des processus à prendre en 
compte, lui permet de contrôler la parcimonie du modèle. D'autre part, l'explicita- 
tion des processus permet, par l'approche mécaniste, de simuler in silico des scénarios 
(p. ex. l'impact de traitements insecticides selon différents protocoles pour contrôler 
les populations de moustiques) et, par leur analyse, d'identifier les points de contrôle 
du système étudié. Des études expérimentales seraient comparativement plus difficiles 
ou coûteuses à mettre en place. 


Depuis l'entrée dans l'ère des gros volumes de données (big data ou données massives), 
l'intérêt de telles approches associées à la compréhension, l'exploitation et la valida- 
tion de mécanismes explicatifs est parfois questionné au regard de l’hégémonie des 
approches pilotées par les données (cf. chapitre 2) et du pouvoir prédictif, parfois 
spectaculaire, des méthodes de machine learning (Baker et al., 2018). Pourtant, les 
deux types d’approches apportent des informations complémentaires (Tran et al., 
2020) et pourraient être avantageusement combinés (Baker et al., 2018). 


Nous présentons dans ce chapitre l'exemple d'un modèle mécaniste devenu outil appli- 
catif dans le cadre de la surveillance entomologique et de la lutte antivectorielle contre 
les moustiques Aedes et qui s'appuie sur la description explicite du cycle de vie du 
moustique. Cet outil est développé dans le cadre du projet Arbocarto!, qui intervient 


1. https:// www.arbocarto.fr 
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dans un contexte de consolidation des stratégies de surveillance et de contrôle mises en 
œuvre par les opérateurs de la lutte antivectorielle dans différentes régions françaises. 
Le modèle du même nom est un outil de cartographie prédictive des densités de popu- 
lations des moustiques Aedes, conçu à une échelle spatiale adaptée à l’organisation des 
actions de surveillance et de contrôle, et s'appuyant sur des données météorologiques 
(températures et pluviométrie quotidienne) et des produits de télédétection à haute et 
très haute résolution spatiale. 


» Un modèle générique construit sur le cycle de vie du moustique 


Tous les moustiques ont en commun un cycle de vie qui commence par un stade aqua- 
tique immature avec trois grandes étapes de développement : œuf, larve, puis nymphe. 
La rapide phase d’émergence aboutissant à la forme adulte marque la transition vers 
le stade aérien, où les moustiques vont alors alterner des comportements de repro- 
duction, d'alimentation, de repos et de dispersion jusqu’à leur mort (cf. chapitre 1). Ce 
développement qui prend la forme d'une métamorphose complète (insecte holomé- 
tabole) permet d'interpréter le cycle biologique des moustiques comme un modèle à 
compartiments (Cailly et al., 2012). A chaque étape du cycle est associée une représen- 
tation de type «boîte» (ou «compartiment ») dans laquelle sont regroupés l'ensemble 
des moustiques partageant le même état morphologique et des comportements 
similaires (p. ex. «Nymphes»). Les boîtes sont reliées entre elles par des flèches qui 
marquent les transitions entre les différents états (figure 8.1). 


Les variations démographiques associées aux compartiments peuvent être décrites par 
un système d'équations différentielles ordinaires (EDO). Dans le modèle générique de 
population de moustiques développé par Cailly et al. (2012), un système EDO permet 
ainsi d'exprimer chacun des 10 états populationnels au cours du temps, sous la forme 
d'autant de bilans démographiques élémentaires : 3 états aquatiques (E œufs, L larves, 
P pupes), l'état émergent (Aem), 3 états nullipares (A1h, Alg, Alo) et 3 états pares 
(A2h, A2g, A20), les adultes femelles étant divisées en fonction de leur comportement 
durant le cycle gonotrophique (4 : en recherche d’hôte, g : transition entre gorgées et 
gravides, o : à la recherche d'un site d’oviposition). La variation du nombre de mous- 
tiques d’un état donné (p. ex. «Larves») sur une période définie (ici, 1 journée) est 
ainsi exprimée comme l'addition d’une fraction de la population issue du comparti- 
ment précédent (p. ex. « Œufs») et les soustractions d’une fraction de la population 
de l'état donné évoluant vers le compartiment suivant (p. ex. «Nymphes») et d'un 
terme de mortalité (figure 8.1). Pour quantifier ces flux d'individus d'un compartiment 
à l’autre, ce modèle utilise des paramètres considérés comme constants dans le temps 
(p. ex. le nombre d'œufs pondus) et des fonctions qui, elles, varient selon les conditions 
météorologiques comme la pluie ou la température (p. ex. les taux de mortalité ou les 
taux de transition d’un stade à l’autre) [Cailly et al., 2012]. Ce formalisme (figure 8.1) 
permet ainsi de synthétiser les connaissances de l'impact des variables météoro- 
logiques et environnementales sur la dynamique d’une population de moustiques et 
d'adapter une structure de modèle générique à une espèce ou une zone géographique 
donnée (Ezanno et al., 2015). 


Comme évoqué dans la section précédente, des choix de simplification sont effectués 
dans ce modèle, par rapport à l'intérêt de formaliser certains mécanismes biologiques 
au détriment d’autres, jugés secondaires. En raison de leur pouvoir de nuisance 
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Mortalité additionnelle 


Mortalité d en déplacement 


Ath Nullipares A2h Pares en 
en recherche d'hôte recherche d'hôte 


Aem Femelles | 
adultes LO A2g Pares gorgées 


émergentes 


Ato Nullipares A2o Pares 
en recherche en recherche 
de sites de ponte de sites de ponte 


= à 1 1 
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Mortalité d & Mortalité d & Mortalité 
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js nm ge de l’environnement 
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Figure 8.1. Représentation schématique du cycle de vie des moustiques selon une approche de 
modèle à compartiments. Adapté de Cailly et al. (2012). 


et/ou de leur importance épidémiologique liés à la prise d’un repas sanguin (piqûre), 
dans ce modèle, l'accent a été mis sur la modélisation de l'abondance de moustiques 
(femelles) adultes susceptibles de prendre un repas sanguin (adultes femelles en 
recherche d'hôtes, somme de A1h et A2h) et de transmettre un agent pathogène après 
une première piqûre infectante (A2h). D'autre part, les principaux états du cycle de 
vie du moustique ont été conservés, afin de pouvoir simuler l'effet d’une action de 
contrôle sur un stade en particulier. 


En revanche, certains processus ne sont pas représentés : 

— les compartiments des moustiques mâles (en raison de l'absence de piqûres 
de la part du mâle); ce compartiment devra cependant être ajouté pour tester 
l'impact d'actions de contrôle comme la technique de l'insecte stérile (TIS), qui 
consiste à lâcher des moustiques mâles stériles afin de contrôler les populations 
(Haramboure et al., 2020); 

— le déplacement des moustiques (qui reste limité à l'expression d'un terme de 
surmortalité pour les moustiques femelles en recherche d’hôtes ou de gîtes de ponte, 
traduisant la prise de risque associée au vol, à la piqûre et à la ponte); ce modèle ne 
permet donc pas de simuler la dispersion des moustiques. 
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» Adaptation du modèle aux espèces Aedes albopictus 
et Aedes aegypti et spatialisation 


La modélisation de l'abondance des moustiques Ae. albopictus et Ae. aegypti présente 
un intérêt majeur dans le cadre de la lutte antivectorielle. À l'échelle mondiale, ces deux 
espèces sont les principaux vecteurs des virus de la dengue, du chikungunya et Zika. 
Sur le territoire français, Ae. albopictus est notamment présent dans l'océan Indien (La 
Réunion et Mayotte) et dans de nombreux départements de France métropolitaine, 
tandis ou Ae, aegypti est observé aux Antilles, en Guyane et dans l'océan Indien. Si 
un utilisateur souhaite simuler les dynamiques d’abondance de l’une ou l’autre de ces 
espèces dans un contexte biogéographique spécifique, il est alors nécessaire d'adapter 
le modèle présenté dans la section précédente, fondé sur le cycle de vie générique 
du moustique. Cette phase d'adaptation s'effectue au travers du choix des valeurs à 
affecter aux différents paramètres du modèle et des fonctions de transition entre les 
différents stades de vie (Ezanno ef al., 2015) ainsi que de la spatialisation du modèle. 


Choix des paramètres et fonctions du modèle de dynamique 
de population de moustique 


Cette étape nécessite une revue de la littérature ciblée sur l'espèce et le contexte envi- 
ronnemental (Ezanno et al., 2015) afin d'identifier les publications citant explicitement 
les grandeurs recherchées (p. ex. Delatte et al. 2009 ou Lacour et al. 2010 pour Aedes). 
Les fonctions de transition entre stades peuvent généralement être établies à partir 
de résultats d'études expérimentales in vitro, avec un contrôle strict des paramètres 
influents (température, pression, humidité) et en faisant varier une seule grandeur à 
la fois (Delatte et al., 2009; Lacour et al., 2010). Ceci permet d'établir une relation 
biostatistique, caractérisant la dynamique de la transition. Si les valeurs de certains 
paramètres ou fonctions ne sont pas disponibles dans cette littérature, les échanges 
avec les entomologistes médicaux et les observations de terrain des acteurs de la lutte 
antivectorielle constituent alors une contribution déterminante et irremplaçable pour 
les estimer en prenant en compte le contexte environnemental régional (p. ex. la date 
de sortie de diapause pour Ae. albopictus en zone tempérée). 


L'ensemble des paramètres (constantes du modèle) et des fonctions utilisés pour modé- 
liser la dynamique de population d’Ae. albopictus (en milieu tempéré et tropical) et 
d'Ae, aegypti sont détaillés dans les publications de Focks et al. (1993) et de Tran et al. 
(2013, 2020). Pour une même espèce (p. ex. Ae, albopictus en milieu tropical), ces para- 
mètres et fonctions pourront être réutilisés dans un contexte géographique différent. 
C'est ainsi que les données expérimentales obtenues sur l'espèce Ae. albopictus à Vile 
de La Réunion (Delatte et al., 2009) ont permis de modéliser la dynamique de popu- 
lation du moustique-tigre sur cette ile (Tran et al., 2020), mais aussi sur l’île Maurice 
voisine (Iyaloo et al., 2021). 


Spatialisation du modèle : La notion de capacité de charge 
de l'environnement 


Au-delà de la dynamique temporelle pilotée par les variations météorologiques, les varia- 
tions spatiales d’abondance des populations de moustiques entre différents territoires, 
ou au sein d'un même territoire, sont contraintes par la distribution des ressources 
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disponibles, en particulier la densité de gîtes potentiels pour la ponte et le dévelop- 
pement des stades larvaires. Ces gîtes sont de faibles dimensions pour Ae. aegypti et 
Ae. albopictus (coupelles, pneus, cavités, etc.) et doivent garantir une mise en eau à 
court, à moyen, voire à long terme en raison de la forte résistance des œufs d’Aedes à 
la dessiccation (cf. chapitre 5). Cette mise en eau, qui peut être due à une action natu- 
relle (précipitations) ou artificielle (arrosage), déclenche l’éclosion des œufs arrivés à 
maturation et le passage à l'état larvaire (cf. chapitre 1). Si l'on considère un espace de 
quelques hectares, la densité de gîtes larvaires peut être agrégée sous la forme d'une 
moyenne spatiale en une seule valeur, représentative d'une «capacité de charge de 
l'environnement ». Cette définition tire son origine de la constante K en écologie fonda- 
mentale, associée au modèle d’équation logistique (Verhulst, 1845), et exprime un terme 
de mortalité densité-dépendant, concernant ici le stade larvaire et de l'émergence de 
l'adulte. Cette capacité de charge KI est un paramètre important et représente le nombre 
maximum de larves que l'environnement peut supporter sans affecter la population elle- 
même. Elle dépend donc directement du nombre de gîtes larvaires disponibles pour une 
zone donnée et de leur production moyenne associée (p. ex. 10 larves par gîte). 


La typologie et la densité de ces gîtes peuvent fortement varier d’une région à l’autre, en 
fonction de la nature de l'environnement (plus ou moins végétalisé, arboré, etc.) et des 
usages et comportements des populations humaines (stockage d’eau, arrosage de plantes, 
stockage de déchets, etc.). Le modèle doit donc permettre à l'utilisateur de paramétrer la 
capacité de charge de l’environnement en fonction d'observations de terrain (Tran et al., 
2020) ou d’estimations basées sur des informations extraites d'images de télédétection 
comme l'occupation du sol, la végétation ou la structure du bâti (cf. chapitre 5), et de 
données externes comme la localisation de maisons abandonnées, de cimetières, etc. En 
raison de la typologie des gîtes des moustiques Ae. albopictus et Ae. aegypti, la capacité 
de charge de l’environnement peut être décomposée en deux termes Kifix et Klvar, selon 
l'origine de la mise en eau des gîtes larvaires : K/fix, anthropique due à l'arrosage (consi- 
dérée comme « fixe», c'est-à-dire invariant au cours du temps) ou Kivar, naturelle due à 
la pluie (considérée comme naturellement variable). 


>» Implémentation, initialisation et simulation des abondances 
de moustiques Aedes 


Implémentation à l'aide du langage de modélisation Ocelet 


Afin de simuler les abondances de moustiques sur un territoire donné, associé a 
la présence de l’une ou l’autre des espèces d’Aedes décrites précédemment, nous 
implémentons le modèle dans un environnement informatique. Limplémentation 
(dérivé de l'anglais to implement : réaliser, mettre en œuvre) consiste à définir et 
écrire le code et les variables informatiques qui vont permettre d’automatiser toutes 
les tâches de calcul visant à simuler la dynamique des populations de moustiques 
pour chaque parcelle d'étude. Dans la modélisation spatiale, nous bénéficions d'un 
cadre conceptuel et pratique particulièrement adapté avec la plateforme Ocelet 
(Degenne et Lo Seen, 2016). Cette plateforme propose à la fois un langage métier, 
précis et générique, pour décrire des phénomènes dynamiques dans un paysage en 
s'appuyant sur la notion de graphes d'interaction, et un ensemble d'outils permettant 
de renseigner les relations entre des entités de nature géographique, compatibles 
avec les systèmes d’information géographique (SIG). 
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Dans Ocelet, un système est représenté par la description d’«entités» qui peuvent 
avoir différentes représentations spatiales (polygones, points, etc.). La dynamique du 
système est simulée à l’aide d'un «scénario » dans lequel l'état des entités va évoluer 
au cours du temps à l’aide de fonctions : ces fonctions peuvent être propres à chaque 
type d'entité («services») ou dépendre d’autres entités avec lesquelles elles sont en 
interaction («relations ») [Degenne et Lo Seen, 2016]. 


Les entrées du modèle sont (figure 8.2a) : 

— la localisation des stations météorologiques (format géographique shapefile) ; 

— les valeurs des températures minimales et maximales et des précipitations journa- 
lières (format texte). Comme précisé plus haut (figure 8.1), ces variables vont impacter 
les taux de transition d’un stade à l’autre ainsi que les taux de mortalité; 

— un «fichier environnemental» (format shapefile) comprenant la localisation des 
parcelles pour lesquelles les densités de chaque stade vont être estimées, caractérisées 
par leur altitude et la capacité de charge de l'environnement (notée KI, voir section 
précédente). 


a) Données d'entrée 


Localisation Je 


des stations @ << TENA 
et données 
météo journalières 


e) Interface utilisateur 


Module de Simulation 


Fichier 

environnemental 

c) Actions E 

de prevention “J < 

ou \ / EU, cu 

de contrôle ocelet 

d) Nouvelles données f) Sorties 

environnementales Fichier géographique ET? 

(p. ex. : occupation des densités de moustiques | 4 E 

du sol, végétation, e AS 

cimetières...) Pour chaque « parcelle », nombre ATN $ 
Wis prédit de moustiques par état = 


Figure 8.2. Schéma de fonctionnement de l'outil Arbocarto. 


(a) Données d'entrée requises pour lancer une simulation, (b) modèle biologique intégré (figure 8.1), 
(c, d) données d'entrée optionnelles, (e) interface (figure 8.4), et (f) sorties générées. 


Chaque parcelle est considérée comme indépendante des autres et, à chaque pas de 
temps : 

— les valeurs de température et de précipitations de chaque station météorologique 
sont lues et attribuées à chaque parcelle par correspondance spatiale (la station la plus 
proche); 

— les fonctions du modèle sont mises à jour pour chaque parcelle; 

— le système EDO (figure 8.2b) est résolu pour chaque parcelle permettant d’estimer 
le nombre d'individus dans chaque stade (E, L, P Aem, A1h, Alg, Alo, A2h, A2g, A20). 


Les sorties du modèle (figure 8.2f) sont, pour chaque parcelle du «fichier environ- 
nemental », le nombre d'individus dans chaque stade à chaque pas de temps. À noter 
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que le modèle prédit les densités de moustiques du jour correspondant à la dernière 
date des relevés météorologiques, et pas les densités futures. Cependant, la prédiction 
des densités des stades aquatiques (P et L) permet d'anticiper de quelques jours les 
densités adultes. 


Initialisation et exemple de simulations 


La phase d’initialisation est réalisée de la manière suivante : affectation arbitraire d'un 
grand nombre d'œufs (p. ex. 1000 000) pour chaque parcelle élémentaire du territoire 
et d'effectifs nuls pour tous les autres stades, puis observation d’une période d’embargo 
de simulation d’au moins un an (les résultats de simulation de la première année ne 
sont pas exploités, mais servent à initialiser de manière plus réaliste la dynamique 
de population des moustiques). Cette phase rappelle une considération importante 
évoquée lors des choix de simplification dans la section 2 : le modèle est conçu pour 
représenter les phénomènes de dynamiques d’une population de moustiques dûment 
établie, une fois stabilisé le processus de colonisation de l’environnement du site étudié. 


La simulation consiste à définir un pas de temps souhaité pour le calcul des abondances 
(minimum : une journée, généralement une semaine), une période de simulation 
cohérente avec la disponibilité des données météorologiques et l'application (ou non) 
d'actions de contrôle des populations de moustiques. L'intérêt de l'approche mécaniste 
se révèle ici dans la facilité d’implémenter la plupart des actions de lutte antivectorielle : 
— la destruction de gîtes se modélise à travers la diminution de la capacité de charge 
de l’environnement pour la ou les parcelle(s) concernée(s) ; 

— l’utilisation de larvicides ou d’adulticides (fumigation) équivaut à augmenter tempo- 
rairement le taux de mortalité du stade ciblé ` larves ou adultes. 


Le modèle permet ainsi deffectuer plusieurs simulations successives (avec ou sans 
actions de contrôle) et de tester in silico l'impact potentiel d'une de ces actions ou de 
la lutte intégrée (combinant plusieurs actions) [figure 8.3]. 


Validation du modèle 


La confrontation entre les données de simulation (pilotées par les valeurs journalières 
de température et de précipitations collectées au niveau des stations météorolo- 
giques correspondantes) et les données entomologiques de terrain a permis de 
valider les sorties du modèle sur plusieurs sites tests en climats tempérés et tropicaux 
(Iyaloo et al., 2021; Tran et al., 2013, 2020). 


Cette étape permet de vérifier la validité des sorties du modèle, avant d’utiliser le modèle 
à des fins prédictives ou pour évaluer l'impact de différents scénarios de contrôle. 


» Arbocarto : une interface adaptée aux actions de la lutte 
antivectorielle 


Le modèle présenté, fondé sur le cycle de vie du moustique, prend en compte à la fois les 
variations des données météorologiques (température et précipitations) pour piloter la 
dynamique temporelle du modèle et un niveau de favorabilité de l’environnement en ce 
qui concerne la densité de gîtes larvaires pour permettre une spatialisation cohérente à 
l'échelle du paysage (le plus souvent des surfaces de l'ordre de quelques hectares). 
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Limplémentation sous la plateforme Ocelet permet de lancer des simulations rapide- 
ment dans un contexte scientifique. Cependant, afin de viser un public d'utilisateurs le 
plus large possible, la prise en main du modèle doit être accessible aux utilisateurs sans 
connaissances informatiques particulières. L'outil Ocelet permet d'exporter un fichier 
exécutable (format jar) dans un environnement d'exécution Java (JRE pour Java Runtime 
Environment). Il reste cependant à pouvoir accéder aux choix des grandes variables du 
modèle : choix de l'espèce, de la période et du lieu de simulation, prise en compte de points 
noirs observés sur le terrain (forte densité de gîtes larvaires) ou d'actions de lutte anti- 
vectorielle. À cette fin, une interface développée en JavaFx? permet de piloter l'ensemble 
du modèle, phases d’initialisation et de simulation, en quelques clics (figure 8.4). 


L'interface permet également à l'utilisateur de modifier le fichier environnemental 
de départ en fonction de la typologie de l’occupation du sol ou de la végétation, ou 
en ajoutant une ou plusieurs couches d’information géographique complémentaire 
comprenant la localisation de sites producteurs de gîtes larvaires (cimetières, garages, 
jardins familiaux, terrasses à plot, ventes de plantes, maisons abandonnées, etc.; 
figure 8.2f); l'ajout de cette information modifie la capacité de charge de l’environne- 
ment de chaque parcelle concernée. 


L'interface est organisée en trois blocs principaux (figure 8.4) : à gauche le menu avec les 
différents modules à activer; au centre, les paramètres à préciser pour chaque module; 
en bas, un espace d'avertissement ou d’information du bon déroulé du traitement. 


Chaque sortie produite par l'application est accompagnée d'un fichier texte retraçant 
les paramètres choisis par l'utilisateur pour la génération de cette sortie. 
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Figure 8.4. Capture de l'interface Arbocarto permettant de piloter le modèle de dynamique de 
population des moustiques Aedes albopictus et Aedes aegypti. 


Le pilotage se fait à travers des modules spécifiques d’initialisation et de simulation, permettant la prise en 
compte d'actions de lutte antivectorielle et d'observations entomologiques sur le terrain. 


Afin de faciliter la prise en main de l'outil par les utilisateurs, des formations ont été 
dispensées. Loutil Arbocarto est donc entré dans sa dernière phase : son déploiement 


2. https://openjfx.io/ 
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opérationnelet son utilisation en routine au sein des structures en charge de la surveillance 
et de la lutte antivectorielle, que sont les EID (ententes inter-départementales pour la 
démoustication), les ARS (agences régionales de santé) et les collectivités départemen- 
tales. Ces structures utiliseront cet outil en complément de leur travail de prospection, 
d'information aux populations et de contrôle en cas d'épidémie. Le projet dispose d'un 
site web? d'information sur les actions en cours et sur les modalités d’accés à l'application 
sous licence libre CeCILL-C. 
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Chapitre 9 


Simulation spatiale du risque de propagation 
de la dengue à partir de modèles 
comportementaux vecteurs et hôtes 


ÉRIC DAUDÉ, SÉBASTIEN REY-COYREHOURCQ, 
ALEXANDRE CEBEILLAC 


Les maladies arbovirales telles que la dengue, le chikungunya et celle transmise par le 
virus Zika constituent un enjeu majeur de santé publique lorsqu'elles atteignent des 
niveaux épidémiques. À l'échelle mondiale, les épidémies sont en lien étroit avec les 
paramètres climatiques (température, précipitations) et les niveaux d'urbanisation 
des territoires, mais la compréhension de leurs dynamiques de propagation reste un 
défi majeur à l'échelle suburbaine. En effet, à cette échelle, l’hétérogénéité environ- 
nementale, les variations spatio-temporelles des densités de vecteurs et la mobilité 
des populations humaines complexifient la géographie de la maladie et de sa diffu- 
sion. La lutte contre ces virus cible principalement leurs vecteurs, essentiellement 
le moustique urbain Aedes aegypti. Cependant, cette lutte est compromise par la 
superficie des zones à traiter, la difficulté d'atteindre les lieux de présence potentielle 
des moustiques ou de leurs gîtes de ponte, le coût et l'effort humain requis. Si les 
stratégies de lutte développées dans les années 1960, qui passaient notamment par 
l'usage massif d’insecticides comme le DDT, sont aujourd’hui en partie inefficaces, 
la compréhension du système pathogène et les avancées technologiques permettent 
d'envisager le développement de stratégies de lutte ciblées, centrées sur le moustique 
vecteur, notamment en l'absence de vaccin (Hoffmann et al., 2011; Seixas et al., 
2019). Cependant, l'identification des zones à cibler pour la lutte antivectorielle dans 
les grandes métropoles reste un défi majeur. 


En épidémiologie spatiale, une des hypothèses les plus fréquentes postule que 
certains lieux spécifiques ont le potentiel d'être des points chauds de contamination 
hyperlocalisés, responsables de la multiplication des cas et à l’origine de processus 
d'amplification virale. Depuis ces lieux, et par la voie des mobilités humaines, 
d’autres lieux contribuent à la propagation virale à plus large échelle, agissant comme 
des nœuds de dispersion qui facilitent la propagation des pathogènes dans la ville 
(Daudé et al., 2015). Comment identifier ces lieux? Ont-ils des spécificités socio- 
environnementales ? Comment évoluent-ils dans le temps et avec quels impacts sur 
la dynamique épidémiologique ? 
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À l'échelle d'une grande métropole de plusieurs millions d'habitants, l'identification 
de ces lieux devrait permettre un contrôle plus ciblé et moins intrusif, réduisant ainsi 
les foyers épidémiques, soulageant les infrastructures de soins en période épidémique 
et diminuant la charge économique pour les individus et la collectivité. Explorer 
le poids de ces facteurs — noyau de contamination locale vs dissémination via les 
mobilités humaines — dans la dynamique globale d'une épidémie n'est cependant pas 
trivial du fait des relations non linéaires qui existent entre tous ces composants et des 
difficultés d'accès ou de collecte des données pertinentes pour évaluer ces différents 
mécanismes. Pour explorer ces deux axes de recherche qui portent l’un sur le poids 
différencié des lieux dans l'exposition au risque de transmission et, l'autre, sur le rôle 
des mobilités individuelles dans la propagation du virus, la modélisation spatiale est 
donc particulièrement adaptée. 


» Modèles individus-centrés et spatialisés 


Différentes approches peuvent être utilisées pour modéliser et formaliser un objet ou 
domaine d'étude Le choix de l’une ou de l’autre est fonction des objectifs assignés au 
modèle, des connaissances associées au phénomène étudié et des données disponibles 
pour le décrire. 


Différents types de modélisation 


Un premier type de modèle s'appuie sur la disponibilité de données quantitatives pour 
en faire un résumé, une description, il s’agit des modèles statistiques. Des modèles 
de régression logistique multinomiale et des statistiques spatiales telles que lauto- 
corrélation spatio-temporelle peuvent être utilisés pour rechercher des liens entre 
des clusters de cas d'une maladie, les caractéristiques socio-économiques et envi- 
ronnementales des quartiers et les modes d'occupation du sol. L'étude des résidus de 
ces modèles permet alors de rechercher des explications locales à ces écarts, pour 
comprendre par exemple pourquoi certains quartiers semblent davantage exposés à la 
maladie que la tendance générale (Zellweger et al., 2017). 


Un second type de modèle s'appuie sur la formalisation de processus pour estimer 
notamment des valeurs qui ne sont pas connues localement ou prédire l'évolution 
d'un système. Un modèle gravitaire va permettre par exemple de calculer les flux de 
navetteurs entre différentes villes en fonction de la masse de ces villes (exprimée par 
le volume de leur population) et la distance qui les séparent. Ce type de modèle peut 
alors être utilisé pour estimer le risque de propagation d’un virus au sein d’un système 
urbain, voire entre pays via les vols internationaux (Salami et al., 2020). 


Enfin, une troisième approche repose sur la modélisation de comportements, c’est- 
à-dire la manière de fonctionner ou de réaliser des actions d’une entité dans son 
environnement. Un modèle de déplacement peut ainsi décrire les mobilités quoti- 
diennes d’un individu, contraintes par les activités qui devront être réalisées. La 
simulation de milliers de ces déplacements permet alors de reproduire les dynamiques 
collectives de mobilité au sein d’un territoire et d'y associer un mécanisme d'inter- 
action interpersonnelle pour simuler la transmission d’un virus (Karl et al., 2014). La 
modélisation par système multi-agent (Treuil et al., 2008) est particulièrement bien 
adaptée pour ce type de problématique. 
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Utilisation des modèles multi-agents en épidémiologie 


Les systèmes multi-agents (SMA) relèvent des méthodes dites «individus-centrées ». 
Ils permettent d’implémenter informatiquement des modèles qui sont ensuite utilisés 
comme des « mondes artificiels » afin de mener des expérimentations proches de celles 
d’un laboratoire. Il est en effet possible, avec ces modèles, de contrôler les paramètres, 
de les modifier et de répliquer des séries d'expériences, par la simulation. L'observation 
et la mesure des évolutions et des changements d'état des agents permettent alors den 
inférer ce qui se passerait dans le monde réel dans des conditions similaires. 


Fruit des recherches en intelligence artificielle distribuée et en robotique, les agents 
d'un SMA interagissent entre eux et avec leur environnement — le plus souvent 
selon des modes de coopération, de concurrence ou de coexistence — dans le but 
de résoudre des problèmes dont la complexité ne permet pas une résolution indivi- 
duelle. En offrant une formalisation explicite de l’espace, du temps, des interactions, 
des comportements et des niveaux d'organisation, ces méthodes trouvent un terrain 
favorable aux questionnements géographiques, notamment dans l'étude exploratoire 
de phénomènes qui se propagent dans l’espace. Il s’agit alors de proposer une forma- 
lisation explicite des mécanismes fondamentaux qui concourent à la dynamique du 
phénomène, ceux-ci étant ensuite implémentés dans un système informatique afin 
d'en éprouver la validité par la simulation. 


Dans les zones endémiques, notre première hypothèse est que les villes se caractérisent 
par un nombre limité de lieux qui présentent des conditions socio-environnementales 
favorables au maintien de populations de moustiques pendant les périodes inter- 
épidémiques. Cette présence peut être suffisante pour assurer une circulation 
locale, continue et de faible intensité des virus. Lors de changements saisonniers 
infra-annuels, tels que la hausse des températures et/ou les pluies de mousson, les 
populations de vecteurs explosent et augmentent le risque global de propagation 
virale. Notre deuxième hypothèse est alors que l’organisation socio-spatiale et les 
discontinuités au sein des villes, qui structurent les mobilités quotidiennes, sont 
alors à l’origine de la propagation des pathogènes qui, à partir de ces quelques foyers 
initiaux, se diffusent ensuite dans la zone urbaine et entraînent des épidémies de plus 
ou moins grande envergure. Il est donc essentiel d'identifier les réservoirs potentiels 
de diffusion des pathogènes, entendus ici comme un ensemble de lieux hyperlocalisés, 
en lien avec les facteurs environnementaux et anthropiques, pour agir efficacement 
contre la diffusion des épidémies. 


L'évaluation du poids épidémiologique de ces zones est cependant un défi. En effet, les 
systèmes de surveillance enregistrent les cas de contamination au lieu de résidence des 
personnes malades, qui n'est potentiellement pas le lieu de contamination. Et tous les 
cas ne sont pas recensés. Ainsi, les études épidémiologiques réalisées à partir de ces 
seules données ont tendance a surestimer le risque local (contagion par proximité) par 
rapport au risque global (diffusion par sauts). C'est dans ce contexte que la modélisa- 
tion et les enquêtes de terrain ciblées peuvent révéler le rôle des micro-environnements 
locaux et de la mobilité humaine dans la dispersion des pathogènes en milieu urbain. Un 
obstacle important dans la modélisation et l’acquisition de données intervient lorsqu'il 
s'agit de considérer lhétérogénéité à haute résolution des environnements urbains et des 
conditions socio-économiques pour étudier leurs effets sur la dynamique épidémique. 
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» Application à la dengue à Bangkok : MO, méthodes et données 


Nous avons développé un ensemble de modèles intégrés au sein du simulateur MO? 
(Daudé et al., 2015). Celui-ci permet de construire un modèle de simulation composé 
de différentes briques (figure 9.1) qui intègrent les connaissances biologiques du 
vecteur et les mobilités humaines quotidiennes, les dynamiques de ces deux modèles 
étant liées aux hétérogénéités socio-environnementales à haute résolution. Le terrain 
d'application de ce modèle est Bangkok, Thaïlande. 


i Remplissage S 
i E | 


Figure 9.1. Les principaux composants du modèle MO?. 
MOMA: le modèle des moustiques. MOMOS : le modèle des humains. MODE: le modèle de l’environnement. 


Nous présentons dans les sections qui suivent la modélisation de l’environnement 
(MODE), du moustique vecteur (MOMA), des hôtes (MOMOS) et les mécanismes 
de transmission d'agents infectieux entre vecteurs et hôtes. La dernière section du 
chapitre est consacrée à l'illustration d'une simulation avec MO? et aux méthodes 
d'exploration de ce type de modèle. 


Modélisation de l’environnement 


Les dynamiques de l’environnement sont simulées par le modèle MODE (Misslin et 
Daudé, 2016). Les agents de ce modèle sont des cellules qui représentent des portions 
du territoire sur lequel le modèle est appliqué. Ces cellules sont décrites par diffé- 
rents paramètres environnementaux qui influencent le comportement d'Ae, aegypti : 
les modes d'occupation du sol (végétation, eau, types de bâti, etc.), le nombre de gîtes 
larvaires potentiels (récipients pouvant contenir de l’eau claire), les températures de 
Yair et les précipitations. Les données qui permettent de quantifier ces attributs sont 
issues de traitements d'images satellite, de relevés de stations météorologiques et 
de savoirs experts. Chaque agent de l’environnement (cellule) dispose également de 
fonctions qui décrivent l'évolution de la quantité d’eau disponible dans les gîtes de 
ponte potentiels (remplissage et vidange) et la variation des températures. Ces varia- 
tions sont calculées à partir des données des images thermiques MODIS (satellites 
Aqua et Terra). Lorsque la couverture nuageuse est trop épaisse, certains pixels sont 
susceptibles de présenter des valeurs nulles. Si cette situation se présente, MODE 
attribue aux zones concernées une valeur de température globale issue d’une station 
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météorologique à Bangkok!. Ce sont ces mêmes stations qui fournissent quotidien- 
nement les précipitations. La résolution spatiale d'une cellule MO? est de 30 x 30m, 
elle correspond à celle des images Landsat 8 (OLI). Cette résolution est compatible 
avec la superficie moyenne de l'aire de dispersion du moustique observée par des 
études de type marquage-libération-recapture (Sheppard et al., 1969). 


MODE permet ainsi de produire de nombreuses données spatialisées et dynamiques 
(figure 9.2) qui sont utilisées par les modéles comportementaux des moustiques et des 
humains. Nous calculons ainsi des valeurs de température de l'air qui permettent de 
repérer des différentiels importants à l'échelle de Bangkok (figure 9.2a; cf. chapitre 3), 
un indice de végétation (figure 9.2b) et des estimations de la distribution des popu- 
lations humaines (figure 9.2c; cf. chapitre 4). Le nombre de gites de ponte pour le 
moustique est également estimé (figure 9.2d). Une approche fondée sur la connais- 
sance terrain et à dire d’experts a été ici privilégiée, car il est coûteux de collecter cette 
donnée sur le terrain. Les densités de gîtes dans l’espace urbain étant très largement 


T° minimale de l’air Indice de végétation 


(au 22/01/2009) (NDVI modifié) 
| 
19°C 26 °C 0 0,82 


Cellules non habitées Cellules sans GLP* 
EZ? Zones en eau E Zone en eau 
ES Zones végétalisées 

20km N : i 
CH Routes be À *GLP = Gîtes larvaires potentiels 


Figure 9.2. Facteurs environnementaux constituant la niche écologique du vecteur. 
(a) La température de l'air estimée. (b) L'indice de végétation. (c) La répartition de la population humaine. 
(d) La répartition des gîtes larvaires potentiels. 


1. https://aqicn.org/city/thailand/bangkok/thai-meteorological-department-bangna/vn/wxqa.com 


123 


Télédétection et modélisation spatiale 


liées à la présence des hommes et à leurs activités (Ooi et al., 2006), nous avons 
construit un indice d’abondance de gîtes de ponte potentiels basé sur les densités de 
ménages présentes localement (Misslin et al., 2016). 


Une partie de ces données spatialisées servent alors d'entrées pour le modèle du mous- 
tique (Misslin et Daudé, 2017). 


Modélisation du moustique vecteur 


Le modèle MOMA est un modèle comportemental intégré au simulateur MO. 
Il permet de simuler le développement de populations vectorielles en lien avec les 
caractéristiques de l'environnement (Maneerat et Daudé, 2016). La figure 9.3 présente 
la conception générale des transitions entre les états d’Aedes (phases aquatique et 
aérienne), la figure 9.4 représente les interactions entre le moustique adulte femelle et 
l'environnement dans le modèle. 


Pour des simulations à l'échelle d’une ville, les phases de développement aquatique sont 
décrites par des modèles de stocks à une résolution spatiale de 30 x 30m, alors que la 
phase aérienne est décrite par un modèle individu-centré (figure 9.3). La femelle Aedes 
pond généralement dans des récipients contenant de l’eau claire. Le développement de 
l'œuf commence au moment où il entre en contact avec l’eau. Cette phase de dévelop- 
pement dure 2 à 3 jours selon les conditions de température de l’eau. Le moustique passe 
au stade larvaire au moment de l'éclosion. Les larves sont capables de se déplacer dans 
l'eau pour respirer et pour se nourrir. Au bout de 4 à 6 jours, la larve passe au stade de 
nymphe. Cette phase de développement dure entre 2 et7 jours en fonction de la tempéra- 
ture. Le développement du moustique est jusque-là totalement contraint par la présence 
d'eau dans les gîtes larvaires ce qui explique le lien fort qui existe entre l'augmentation des 
densités de vecteurs et les précipitations. Une fois la phase aquatique terminée, le mous- 
tique passe à son stade aérien. Ici le stock de nymphes arrivées à maturité déclenche la 
création d'agents moustiques, seules sont simulées les femelles. Le prélèvement de sang 
leur permet d'obtenir l'énergie dont elles auront besoin lors de la phase d’oviposition (O). 
Les femelles entament leur cycle gonotrophique (Gn) au cours duquel elles sont amenées 
à piquer des mammifères (essentiellement des humains). 


Nymphe 


U 
U 
1 
= Cycle 
Larve (>) Œuf gonotrophique Oviposition 
i e ý Tel 


Stocks aquatiques 


États adultes 


| 


Figure 9.3. Diagramme états-transitions du cycle de vie d'un agent Aedes, selon son stade de 
développement (immature ou aquatique, adulte ou aérien). 


Les fléches représentent les possibilités de transition d’un état vers un autre état. 
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MOMA, le modèle individu-centré du moustique, modélise donc la femelle adulte 
Ae. aegypti responsable de la transmission de virus lors de ses repas sanguins. Les acti- 
vités qui modélisent le comportement du moustique sont (figure 9.4) : s’alimenter, se 
reproduire, se reposer, se déplacer et survivre. Ces activités sont liées à des besoins 
qui peuvent être assouvis selon la présence de ressources en sang et en nectar pour 
l'alimentation, de gîtes larvaires pour la reproduction, de zones ombragées pour le 
repos et de zones de circulation pour les déplacements. La survie quant à elle dépend 
pour l'essentiel des conditions de température, d'humidité et d'un paramètre stochas- 
tique qui factorise les dangers auxquels le moustique est exposé au quotidien (prédation 
notamment). Aedes agit donc selon son état biologique, défini par ses variables internes, 
et selon les ressources qui sont présentes dans son environnement immédiat. 


Pondre des œufs dans 


Canaux de communication 
entre les deux agents 
J + 
B1-s’alimenter E … E éi 
, + — nectar , 
B2 - se reproduire 4 1 C3- gîte larvaire alimente 
B3-sereposer + | C4- ombre 
FA B4 - se déplacer & 
TA B5 - survivre pL 7 C5-barrière 
gë E =+ C6- température / 
humidité 
+ être actif z jè 
ii 4 + C7 - lumière Stocks 
— m immatures 
a 


C Création d'un nouvel agent femelle moustique Aedes ee e DE 


Figure 9.4. Modélisation des interactions entre le moustique femelle et l'environnement. 


Modélisation des hôtes et de leurs mobilités 


Le modèle MOMOS décrit le composant hôte, correspondant pour l'essentiel au 
comportement de mobilité d'un humain (Cebeillac et al., 2017). Un agent hôte est 
localisé dans l’environnement dans lequel il se déplace pour réaliser ses activités. 
Les déplacements et changements de position dans l’espace exposent l'humain à 
différents risques (figure 9.5). Dans le modèle, les déplacements s'effectuent entre 
différents espaces d'activités (maison, école, lieu de travail, de loisir, etc.) sur la base 
d'un agenda de mobilité construit pour chacun des agents. Différentes sources de 
données sont mobilisées pour construire ces agendas : les réseaux sociaux (Twitter, 
Facebook) et les feuilles de temps issues d'enquêtes (Cebeillac et al., 2018). La géné- 
ration des agents est réalisée à partir d’un générateur de population synthétique (voir 
chapitre 4), chaque agent étant caractérisé par des attributs géographiques, tels qu'un 
lieu de domicile et un agenda d'activités journalier. Les activités réalisées par un indi- 
vidu, qui déterminent les lieux fréquentés, dépendent dans le modèle de son âge, 
de son lieu de résidence et de son statut épidémiologique (un agent présentant des 
symptômes de dengue ne se déplace pas). 


Mécanismes de transmission du virus 


Les agents de type moustique (m) et hôte (h) sont décrits selon leur exposition au virus 
(figure 9.6), comme respectivement susceptibles (Sm, Sh), exposés (Em, Eh) et infectés 
(Im, Ih). Seul l'agent hôte peut être immunisé (Rh), l’agent moustique restant, jusqu'à 


125 


Télédétection et modélisation spatiale 
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Figure 9.5. Représentation de l'espace d'activité défini par Hägerstrand (1970) en contexte de 
santé. Dans cet exemple, quatre personnes d'un même foyer fréquentent différents lieux au 
cours d’une journée. Chacun de ces lieux est associé à un risque sur la santé selon l’activité qui 
y est réalisée et l'heure de la journée. D’après Cebeillac et Daudé (2022). 


Schéma de transition 


I D ‘as 
de l'humain 


Figure 9.6. Schéma de transition entre les différents états pour un agent moustique (m) et un 
agent hôte (h). 

L'agent-hôte, dans l’état susceptible (Sh), peut être contaminé par le virus de la dengue à la suite d'une 
piqûre de femelle moustique infectée (Im) avec un risque de transmission (a). Une fois contaminé, l'hôte 
développe le virus durant une période (T1) appelée incubation intrinsèque puis entre dans sa phase de 
virémie (T2), phase durant laquelle il peut contaminer un moustique sain qui le piquerait (Sm) avec une 
probabilité (y). Les moustiques Aedes contaminés entrent dans une période d’incubation extrinsèque (EIm) 
avant de devenir contagieux (Im). 
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sa mort, infecté donc contagieux (Im). La transmission du virus d’un hôte infecté vers 
un moustique sain ou d'un moustique infecté vers un hôte sain se fait au cours d'un 
repas de sang pris par le moustique. 


Ces différents états peuvent ou non modifier le comportement des agents. Par exemple il 
est possible de limiter, voire de stopper, les mobilités des agents-hôtes en cas de contami- 
nation afin de prendre en compte l'apparition de symptômes plus ou moins incapacitants. 


» Simulations du modèle MO? 


Les simulations du modèle MO? sont multi-temporelles, un pas de temps (itération) 
est fixé à une minute. Le temps d’une simulation se décompose ainsi en trois échelles 
selon les processus : la journée, l'heure et la minute. La mise à jour des variables 
environnementales et le développement aquatique du moustique se font en fonction 
des précipitations, de l'évaporation et des températures moyennes journalières. Les 
mobilités des agents humains sont simulées par tranche horaire en fonction de leurs 
agendas et les activités des moustiques adultes en fonction des températures, de la 
luminosité et de leurs besoins à la minute. L'enregistrement de différents paramètres à 
chaque itération permet de reconstruire et d'analyser l'évolution de nombreuses dyna- 
miques. Il est ainsi possible de suivre l’évolution des stocks de moustiques en phases 
aquatiques et de calculer des densités de moustiques adultes à la résolution spatiale de 
la cellule. Il est également possible d'identifier les différents lieux d'interactions entre 
les hôtes et les vecteurs, et de calculer des profils environnementaux à risque lorsque 
ces interactions aboutissent à des boucles de contamination. L'interface de simula- 
tion du modèle permet de visualiser l'évolution de différents indicateurs, sous formes 
graphiques et cartographiques (figure 9.7). 


{È Simulation:mo3_modetMO3madel - MO3/madels/mo3.gaml - Gama - oa 
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Figure 9.7. Interface de simulation du modèle MO. 
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Des dynamiques de vie des moustiques 


Au cours d’une itération, chaque minute, les agents moustiques explorent leur 
environnement et peuvent sélectionner une cible qui répond à leurs besoins et 
se déplacer. Si aucune cible nest identifiée, car aucune ne peut satisfaire l’un des 
besoins exprimés, le moustique se déplace au hasard. Si une cible est sélectionnée, 
le moustique se déplace dans sa direction. Dans le mécanisme de déplacement, le 
coefficient de porosité associé à l'environnement joue un rôle important, car il va 
limiter, par exemple, le passage entre un objet spatial de type route et une maison 
(effet barrière), alors qu'il va faciliter le passage d'une route à un espace ouvert 
(effet corridor) [Maneerat et Daudé, 2016]. Nous avons par exemple simulé les 
dynamiques de populations de moustiques et leurs dispersions spatiales à partir de 
différentes configurations géographiques. L'objectif était de mesurer les effets des 
catégories d'usage du sol et des densités humaines sur la dispersion de moustiques 
depuis différents lieux d’émergence (figure 9.8). 


Ces simulations ont permis d’une part de vérifier la cohérence des distances 
moyennes parcourues par des cohortes de moustiques virtuels en zone urbaine au 
regard des données de terrain et de la littérature et, d'autre part, de mettre en relation 
ces distances relativement faibles (quelques dizaines de mètres et jusqu’à plusieurs 
centaines de mètres pour les plus élevées) avec les configurations locales d'usage du 
sol et les densités humaines. Ainsi, la zone très peuplée et au bâti très dense du bidon- 
ville (figure 9.8, au centre de la zone) connaît une faible dispersion des cohortes de 
moustiques, au contraire des zones plus périphériques, moins denses et caractéri- 
sées par des successions de résidences et de jardins dans lesquels les moustiques, de 
proche en proche et en fonction de leurs besoins, peuvent parcourir de plus grandes 
distances au cours de leur vie. 


0,125 0,25 


Figure 9.8. Motifs de déplacement d’agents-moustiques dans un quartier de Delhi (Inde). 


Chaque tache de couleur correspond à la trace des déplacements d'un moustique au cours de sa vie. 
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Plus de 50 paramètres et variables sont utilisés dans le modèle MO? développé dans le 
langage de programmation multi-paradigme Scala?. Ce niveau de détail rend la cali- 
bration et la validation particulièrement complexe. En effet les combinaisons possibles 
de paramètres rendent quasiment impossible l'exploration systématique de leurs 
domaines de variation. La validation du modèle moustique a par exemple été réalisée 
selon deux critères : le savoir expert et les données d'enquêtes (Maneerat et Daudé, 
2017). Les recherches se poursuivent pour étudier les comportements de ce modèle, le 
calibrer et explorer différents scénarios de lutte antivectorielle. 


Exploration du comportement du modèle 


Les données de calibration utilisées sont des données issues de la surveillance vecto- 
rielle et des données épidémiologiques, les unes devant permettre de calibrer et de 
valider le modèle moustique, les autres d’éprouver l'hypothèse du rôle des mobilités 
humaines dans la dispersion de la maladie à l'échelle d’une ville. La reconstruction d’un 
monde virtuel à partir d'un certain nombre de proxys et de faits stylisés observés dans 
la réalité produit des résultats qui nécessitent malgré tout une évaluation pour pouvoir 
être rapportés à la réalité. En effet, un modèle à lui seul ne saurait capturer l’ensemble 
des phénomènes complexes à l'œuvre dans la réalité. Une validation à l'échelle micro et 
macro des processus est nécessaire pour s'assurer de la cohérence et de la robustesse des 
résultats produits par ces données issues de la simulation générée en vue de répondre, 
non pas à toutes les questions, mais bien à cet ensemble de questions spécifiques pour 
lequel les mécanismes et paramètres ont été sélectionnés. Les modèles dans MO? mobi- 
lisent chacun un grand nombre de paramètres et de mécanismes en interaction, ce qui 
génère en sortie une variété et un volume de comportements et de données impossible à 
énumérer et donc à évaluer sans passer par des méthodes plus ciblées. 


Pour faire face à cette combinatoire gigantesque — sur des modèles dont une exécu- 
tion dure potentiellement plusieurs minutes voire plusieurs heures —, nous faisons 
appel à des algorithmes d'exploration (algorithmes évolutionnaires) et à la plateforme 
OpenMOLE (Reuillon et al., 2013) capable d'exploiter le parallélisme offert par les 
infrastructures HPC (High Performance Computing). Les modèles MO? bénéficient 
ainsi à la fois des infrastructures de la grille de calcul européenne (4000 cœurs) et des 
moyens de calcul intensif haute-performance du mésocentre régional de Normandie, 
le CRIANN (10000 cœurs), et cela pendant plusieurs semaines, afin d'exécuter et 
d'évaluer la cohérence et les résultats de simulations sur ces questions. 


» Conclusion 


MO? est un simulateur basé sur des modèles comportementaux à haute résolution 
spatiale et temporelle. Il décrit un système pathogène complexe pour étudier les 
facteurs environnementaux et anthropiques de la diffusion de la dengue. En phase de 
développement, MO? est utilisé tout au long du projet du même nom? pour guider 
la collecte de données et les expérimentations sur le terrain et, à terme, pourrait être 
employé pour identifier les stratégies de contrôle vectoriel les plus efficaces pour une 
évaluation in situ de celles-ci à partir des algorithmes d'optimisation. 


2. https://www.scala-lang.org/ 
3. www.MO3.cnrs.fr 
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Les maladies vectorielles dont sont responsables les moustiques, comme le palu- 
disme, constituent toujours un fardeau pour les systèmes de santé et le développement 
de nombreux pays tropicaux, notamment en Afrique subsaharienne. En outre, ces 
derniéres décennies, de nombreuses maladies vectorielles dues aux moustiques, comme 
la dengue, le chikungunya, le virus Zika, la fièvre du Nil occidental et la fièvre de la vallée 
du Rift ont émergé ou réémergé de par le monde, gagnant de nouveaux territoires et 
affectant aussi bien la santé animale que la santé humaine. La persistance ou l'émergence 
de ces maladies vectorielles est conditionnée par de nombreux facteurs climatiques, 
environnementaux, démographiques, propres à chaque territoire, influençant a diffé- 
rentes échelles leur fonctionnement. Depuis deux décennies, les approches géomatiques 
sont utilisées pour mieux comprendre les processus épidémiologiques des maladies 
vectorielles et mieux prédire leurs évolutions dans un monde en plein bouleversement. 


Cet ouvrage dresse un panorama original et actualisé des applications potentielles de 
télédétection et de modélisation spatiale pour mieux surveiller et contrôler les mala- 
dies infectieuses dont l'agent pathogène est transmis par des moustiques. La présence, 
l'abondance et le comportement des moustiques vecteurs et, de fait, l’épidémiologie des 
maladies qu'ils transmettent sont des phénomènes complexes dépendant de différents 
facteurs géographiques, climatiques et environnementaux, y compris de facteurs liés 
à la démographie et aux activités humaines, qui conditionnent localement l'expression 
de ces maladies dans l’espace et le temps. 


Les techniques de télédétection et de SIG couplées à différentes approches de modéli- 
sation appliquées aux maladies vectorielles présentent de multiples atouts : 

— cartographier la distribution des gîtes larvaires et les densités de populations de 
moustiques adultes; 

— évaluer à différentes échelles spatiales le risque vectoriel, pour des maladies déjà 
présentes ou risquant d'être introduites dans un territoire; 

— déterminer des lieux et des saisons à haut risque de transmission (hot spots) afin 
d'optimiser la surveillance et le contrôle; 

— déterminer les facteurs environnementaux, démographiques et socioculturels qui 
exposent davantage certaines communautés aux risques d'infection afin de cibler la 
prévention; 
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— prendre en compte les changements humains et les modifications environnemen- 
tales induites, notamment les transformations des terres agricoles et de la couverture 
végétale, ’étalement urbain et l’industrialisation, afin de prévoir à court et moyen 
termes l’évolution du risque vectoriel; 

— évaluer l'efficacité et aider à optimiser différentes méthodes de lutte antivectorielle, 


prises isolément ou en combinaison dans le cadre d’une stratégie de lutte intégrée. 


Les informations dérivées des données satellites utiles pour mieux appréhender 
le risque vectoriel et décrites ici dans la partie 1 sont disponibles facilement et à 
moindre coût pour les chercheurs et les gestionnaires de la santé. Ces informations 
sont intégrées avec d’autres données climatiques, entomologiques, démographiques 
et épidémiologiques dans des plateformes SIG permettant d'appréhender le fonction- 
nement complexe d'une maladie vectorielle tenant compte de tous ses déterminants et 
de leur dimension spatiale. 


Différentes approches de modélisation présentées dans la partie 2 de cet ouvrage sont 
alors utilisées pour mieux comprendre le fonctionnement d'une maladie vectorielle 
dans un contexte donné et mieux prédire son évolution dans l'espace et le temps, tenant 
compte des interventions de lutte antivectorielle ou des changements environnemen- 
taux et climatiques. Ces approches diffèrent selon les objectifs, les connaissances 
et les informations disponibles sur un vecteur et une maladie vectorielle donnée. 
La partie 2 présente ainsi ces approches diverses mais souvent complémentaires 
(modèles de distribution d'espèces, modèles basés sur les connaissances, modèles 
mécanistes, modèles comportementaux) déployées sur des maladies vectorielles 
variées (couples Anopheles/paludisme, Aedes/dengue) dans des contextes différents 
(Guyane, Madagascar, La Réunion, Thaïlande) et permettant de développer des outils 
opérationnels utiles aux décideurs et aux gestionnaires de santé pour optimiser la 
surveillance entomologique et la lutte antivectorielle. 


À l'échelle d’un territoire insulaire tropical soumis à différents risques vectoriels 
comme l’île de La Réunion, des réussites récentes illustrent tout le potentiel de ces 
approches de modélisation appliquées comme une boîte à outils. La caractérisation 
spatiale et saisonnière du risque lié au typhus murin ou aux «bavites»! permet de 
mettre en place des campagnes de sensibilisation ciblées et adaptées auprès des profes- 
sionnels de santé ou des éleveurs (Grimaud et al., 2021, 2019; Tran et al., 2021). Des 
outils cartographiques de prédiction des densités vectorielles et du risque de transmis- 
sion de la dengue sont appropriés et utilisés depuis plusieurs années par l'opérateur de 
lutte antivectorielle (ARS La Réunion) pour optimiser la surveillance entomologique 
et les actions de prévention et de contrôle (Tran ef al., 2020). Des outils de modéli- 
sation sont développés et transférés pour évaluer et optimiser les méthodes de lutte 
antivectorielle seule ou en combinaison dans un contexte de stratégie de lutte intégrée 
(Douchet et al., 2021). Ces réussites sont favorisées à La Réunion par la disponibilité et 
la richesse des connaissances et données expérimentales et observationnelles ainsi que 
la proximité locale de chercheurs issus de différentes disciplines et travaillant étroite- 
ment avec les bénéficiaires de ces travaux de recherche comme ARS La Réunion au 
sein d’un réseau comme le Dispositif de recherche et d'enseignement en partenariat 


1. La fièvre catarrhale ovine et la maladie hémorragique épizootique des cervidés sont deux maladies 
animales virales touchant régulièrement l’île de La Réunion. Ces «bavites», surnom local tiré des symp- 
tômes de salivation excessive qu’elles induisent, sont responsables de pertes économiques importantes. 
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One Health Océan Indien’. La flexibilité et l'adaptabilité de ces approches permettent 
néanmoins de travailler dans d’autres contextes moins documentés comme cela est 
montré dans cet ouvrage à travers les travaux menés dans la grande île voisine de 
Madagascar (Rakotoarison et al., 2020). 


Comme nous l'avons vu avec l'exemple de l’île de La Réunion, la convergence de 
facteurs tels que la disponibilité de données satellites multi-temporelles et de données 
entomologiques et épidémiologiques géoréférencées, la collaboration entre les scien- 
tifiques de la télédétection, les modélisateurs et les biologistes, et la disponibilité de 
systèmes statistiques et d’algorithmes de traitement d'images adaptés créent un envi- 
ronnement de recherche fertile. L'association des bénéficiaires en amont des travaux 
de recherche et une politique réelle de transfert des résultats issus de ces travaux à 
travers la coconstruction d'outils opérationnels de gestion facilitent leur appropriation 
et leur utilisation effective dans les programmes de surveillance et de lutte antivec- 
torielle. Ces approches de télédétection et de modélisation spatiale relativement 
génériques et transversales se déclinent aussi bien sur des maladies vectorielles zoono- 
tiques (fièvre de la vallée du Rift, fièvre du Nil occidental, typhus murin), humaines 
(paludisme, dengue) qu’animales («bavites») et associent naturellement les sciences 
environnementales et humaines. On retrouve ici les piliers de l'approche «une seule 
santé», y compris la multidisciplinarité et l'intersectorialité. 


Le développement de la géomatique appliquée aux maladies vectorielles a connu un 
engouement certain ces dernières décennies. Un large choix de puissants logiciels de 
traitement d'images, de SIG, d'outils statistiques et d’approches de modélisation, comme 
illustré dans cet ouvrage, est désormais disponible dans un environnement informatique 
de bureau abordable, ce qui permet aux épidémiologistes et aux biologistes d’expé- 
rimenter de nouvelles techniques d'analyse spatiale. Bien que les études examinées 
dans ce document démontrent l'efficacité de la télédétection et d’autres technologies 
géospatiales dans la surveillance et le contrôle des maladies vectorielles, plusieurs para- 
mètres doivent être pris en compte pour que ces technologies soient adoptées de façon 
routinière dans la gestion de la santé publique, notamment dans les pays du Sud. Il 
s’agit entre autres de la disponibilité des ressources pour la collecte, le traitement et 
la modélisation des données géospatiales, de la formation du personnel à l'acquisition 
et l'interprétation correcte des résultats, de la rentabilité de ces techniques de surveil- 
lance et de la disponibilité continue des données de télédétection en temps opportun. 
Il convient également de souligner que l'affectation de ressources à ces nouvelles 
techniques de surveillance ne doit pas se faire au détriment des activités de base de 
prévention et de gestion des maladies au niveau communautaire. 


Pour renforcer les capacités dans le domaine de la géomatique appliquée aux mala- 
dies vectorielles et construire des collaborations interdisciplinaires et intersectorielles, 
il convient de souligner l’importance des formations : les formations diplômantes 
(masters, théses) de jeunes étudiants permettent de développer la masse critique des 
compétences sur ces thématiques, notamment dans les pays du Sud. Lintégration de 
modules spécifiques à l'utilisation de ces approches dans la formation initiale des futurs 


2. https://www.onehealth-oi.org/ 

3. L'initiative One Health (une seule santé) est portée par l'OMS, l'Organisation mondiale de la santé animale 
(OMSA), l'Organisation des Nations unies pour l’agriculture et l'alimentation (FAO) et le Programme des 
Nations unies pour l’environnement (PNUE). 
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acteurs/décideurs de santé et l’organisation de formations professionnelles permettent 
de partager l'avancée des travaux de recherche sur l'apport des données d'observation 
de la Terre et des modèles dans les domaines de la santé humaine, de la santé animale 
et des écosystèmes, en soulignant leurs potentiels et leurs limites. Le développement 
récent des formations en ligne tend à favoriser cette appropriation par les utilisateurs 
des produits et méthodes de télédétection et de modélisation. 


Les perspectives de l’utilisation de l'information spatiale, et en particulier des images 
d'observation de la Terre, pour mieux comprendre, prédire et prévenir les risques de 
transmission de maladies vectorielles, sont nombreuses. Tout d’abord, et malgré leur 
développement exponentiel, l'application de ces techniques de télédétection et de modé- 
lisation spatiale reste à faire pour beaucoup de maladies vectorielles dans de nombreux 
territoires. L'exactitude des cartes prédictives de ces maladies doit aussi être vérifiée sur 
le terrain : on aura toujours besoin de spécialistes en entomologie et en épidémiologie. 
Outre les approches biologiques, les modèles sociaux et comportementaux tels que le 
temps passé par les individus à l'extérieur — qui augmente le risque d'exposition aux 
vecteurs anthropophiles — les types de construction des maisons, l’utilisation de mous- 
tiquaires et autres répulsifs, ainsi que la disponibilité d'installations sanitaires de base et 
de soins de santé primaires — qui sont liés aux conditions socio-économiques —, sont 
importants pour la prévention et le contrôle des maladies vectorielles. Ainsi, il a été 
démontré que les caractéristiques de ces maladies sont étroitement liées à la pauvreté 
et aux inégalités sociales. Ces facteurs socio-anthropologiques et entomologiques ne 
peuvent être déduits des seules techniques de télédétection. 


Une des pistes intéressantes qu'ouvre la modélisation est la possibilité de tester des 
scénarios futurs de changements et d'étudier leurs potentiels impacts sur les risques 
sanitaires. Par exemple, l'impact des changements climatiques sur la distribution 
et l'abondance des populations de vecteurs, ainsi que sur le risque de transmission 
d'agents pathogènes, peut être étudié par la modélisation (Guis et al., 2012; Kraemer 
et al., 2019). En effet, la plupart des modèles présentés dans la partie 2 intègrent les 
variables température et précipitations : utiliser comme entrée de ces modèles les 
projections de ces deux variables (augmentation des températures, modification 
des régimes de pluies) selon les scénarios du Groupe d'experts intergouvernemental 
sur l'évolution du climat (GIEC) permet de tester in silico leurs effets. D’autre part, 
l'impact d'actions de lutte contre les vecteurs peut également être testé par les modèles 
(Douchet et al., 2021; Haramboure et al., 2020) et discuté avec les acteurs de la LAV 
pour optimiser leurs actions, notamment dans le cadre de la mise en place de stratégies 
de lutte intégrée. C’est d’ailleurs un des principaux besoins d’évolutions exprimés par 
les utilisateurs de l'outil Arbocarto que sont les opérateurs de la lutte antivectorielle en 
France (cf. chapitre 8). Enfin, d’autres types de scénarios peuvent étre testés, comme 
Yimpact des changements démographiques (populations, mobilités), mais aussi des 
changements environnementaux. Dans ce dernier exemple, l'intérêt de la modélisa- 
tion est de permettre aux décideurs la visualisation sous forme de cartes de l'impact 
des mesures d’aménagement ou d’atténuation, afin d’éclairer la prise de décision en 
intégrant des enjeux de santé et notamment la prévention des risques épidémiques. 


Enfin, si les données spatiales et les méthodes présentées dans cet ouvrage portent 
sur les maladies dont sont responsables les moustiques, les approches présentées 
(extraction d’informations à partir de l’imagerie satellite, modélisation de distribution 
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d'espèces, de dynamiques de population, de mobilités) ont un caractère générique. 
Elles peuvent donc être adaptées à d’autres questions et enjeux, comme les enjeux de 
biodiversité (et l'étude de leurs liens avec la santé), mais aussi les enjeux de sécurité 
alimentaire (avec des questions communes liées à l’épidémiologie végétale ainsi qu’à 
la distribution et à la dynamique des ravageurs de culture). Ces questions rejoignent 
l'approche «une seule santé» qui, au-delà de la santé animale et humaine, doit 
également intégrer la santé des plantes et des écosystèmes. 
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Agent pathogène : agent (virus, bactérie, protozoaire, nématode, etc.) capable dengendrer une 
lésion ou de causer une maladie chez les animaux (y compris l'humain) ou chez les plantes. 


Arbovirus : contraction des termes anglais Arthropode Born Virus, virus transmis par des 
arthropodes. On appelle alors arbovirose la maladie causée par cet arbovirus. 


Arthropode : groupe d'animaux à corps segmenté et recouvert d'un exosquelette, tels que les 
tiques et les moustiques. 


Autocorrélation spatio-temporelle : mesure statistique qui répond à la question suivante : 
les valeurs d’entités spatiales proches évoluent-elles dans le temps de la même manière que les 
valeurs d’entités lointaines ? 


Biais : Le biais statistique d’un estimateur est la différence entre la valeur de l'espérance de cet 
estimateur et la valeur qu'il est censé estimer. Le terme «biais» recouvre aussi l'ensemble des 
démarches ou procédés qui engendrent des biais statistiques, c’est-à-dire des erreurs dans les 
résultats d'une analyse ou d'une modélisation. 


Biotope : espace écologique de taille variable qui offre des conditions viables, généralement 
assez homogènes et constantes, à l’ensemble des êtres vivants qui y coexistent. 


Cas index : en épidémiologie, ce terme désigne la première personne d'une épidémie (ou d'un 
cluster épidémiologique) à avoir été contaminée par un agent pathogène. 


Cluster (agrégat) : regroupement d'objets reliés dans le temps et dans l'espace. En épidémio- 
logie, on parle de foyer de contagion pour définir un événement sanitaire qui lie entre eux des 
individus. 


Compétence vectorielle : estimation quantitative de l'efficacité du vecteur à transmettre un 
agent pathogène en conditions naturelles. 


Diapause : état physiologique de repos, génétiquement déterminé, pour un organisme, qui 
diminue l'intensité de ses activités métaboliques. Souvent saisonnier, il peut être induit par 
exemple par la photopériode (Duvallet et al., 2017, voir références du chapitre 1). 


Dormance : vie ralentie sans croissance ni développement, minimisant les risques en relation 
avec des conditions environnementales défavorables et dont le déterminisme est environ- 
nemental (Duvallet et al., 2017, voir références du chapitre 1) 


Espérance mathématique : l'espérance d'une variable aléatoire est une valeur numérique 
représentant le résultat moyen d’une expérience aléatoire. C’est l'équivalent en probabilité de la 
moyenne d'une série statistique en statistiques. 


Estivation : période pendant laquelle certains arthropodes réduisent leur métabolisme et 
puisent dans leurs réserves de graisse du fait de conditions environnementales (température, 
hygrométrie) défavorables. Selon la saison concernée, on peut aussi parler d’hivernation. 


Eurygame : caractérise une espèce dont l'accouplement a lieu au cours d’un vol dans un grand 
espace; l'élevage en cage est donc difficile, voire impossible. Par opposition, sténogame désigne 
une espèce dont l’accouplement a lieu dans un espace restreint. 
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Gîte larvaire : espace dans lequel se développent les stades aquatiques des moustiques, égale- 
ment associé au lieu de ponte des femelles. Les gîtes larvaires varient considérablement selon les 
espèces et leurs besoins écologiques aux stades aquatiques : eau douce ou saumâtre, exposition 
au soleil ou à l’abri, permanents ou temporaires, chargés de débris végétaux ou composés d’eau 
claire, etc. 


Modèle gravitaire : en géographie, les modèles d'interaction spatiale cherchent à décrire les 
relations entre les lieux. Le modèle gravitaire permet de caractériser l'intensité de ces relations 
entre unités géographiques en tenant compte de leur potentiel (population, richesse, etc.) et de 
leur distance. 


Nullipare : désigne une femelle qui n’a pas encore pondu, par opposition à «pare ». 


Oviposition : (du lat. ovum, œuf), action de ponte des œufs à un endroit choisi par la 
femelle gravide (en lien avec l'écologie larvaire), pour permettre leur développement dans 
l'environnement. 


Période d’incubation extrinsèque : temps nécessaire aux étapes de développement du virus 
ingéré pour rendre le vecteur infectant lors du repas de sang suivant. 


Raster : un raster est une grille, aussi appelée matrice, organisée en lignes et en colonnes. 
Chaque cellule de cette grille est un pixel unique associé à des valeurs. C'est la distribution de 
tous les pixels, dans la grille (ligne i, colonne j), qui crée une image (en général au format «.tif»). 


Résidus : Dans un modèle statistique, les résidus correspondent à la part des données 
d'observation qui nest pas expliquée par le modèle. Dans un modèle de régression, les résidus 
correspondent à l'écart entre les données estimées par le modèle et les données observées. 


Vecteur (biologie) : se dit d’un organisme qui transmet un agent infectieux. En particulier, il 
peut s'agir d’un arthropode hématophage qui assure la transmission biologique (ou mécanique) 
active d’un agent pathogène d'un vertébré à un autre vertébré. 


Vecteur (géomatique) : La donnée vecteur est utilisée pour représenter les entités du monde 
réel dans un SIG. Une entité vectorielle peut avoir une géométrie de type point, polyligne ou 
polygone. Chaque entité vecteur a des données attributaires qui la décrivent. La géométrie de 
l'entité est décrite en termes de sommets!. 


1. https://docs.qgis.org/2.8/fr/docs/gentle_gis_introduction/vector_data.html#:~:text=La%20 
donn%C3%A9e%20vecteur%20est%20utilis4%C3%A9e,d%XC3%A9IcriteX20en%20termes%20de%20sommets 
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ACP 
AHP 
AID 
AMC 
AMCS 
ARS 
BI 
CES 
Cnes 
Cirad 
DDT 
DGS 
DLR 


EDO 
ECOSTRESS 
EID 

ESA 

ESRI 

ETM+ 

EVI 

FAO 


FFT 
FOTO 


GBIF 


GHSL 
GIEC 
GLCM 
GLM 
GPM 


Liste des sigles 


Analyse en composantes principales 

Analytic Hierarchy Process / Processus d’analyse hiérarchique 
Aspersion intra-domiciliaire d’insecticide 

Analyse multicritére 

Analyse multicritére spatialisée 

Agence régionale de santé 

Brightness Index / Indice de brillance 

Centre d'expertise scientifique (du pôle Theia) 

Centre national d’études spatiales 

Centre de coopération internationale en recherche agronomique pour le développement 
Dichlorodiphényltrichloroéthane (insecticide) 

Direction générale de la santé 


Deutsches Zentrum für Luft- und Raumfahrt / Centre allemand de l’aéronautique 
et astronautique (agence spatiale allemande) 


Equations différentielles ordinaires 

Ecosystem Spaceborne Thermal Radiometer Experiment on Space Station 
Entente interdépartementale pour la démoustication 

European Spatial Agency / Agence spatiale européenne 

Environmental Systems Research Institute (entreprise privée) 

Enhanced Thematic Mapper Plus 

Enhanced Vegetation Index / Indice de végétation amélioré 


Food and Agriculture Organization / Organisation des Nations unies 
pour l'agriculture et l'alimentation 


Fast Fourier Transform | Transformée de Fourier rapide 


Fourier-based Textural Ordination / Ordination texturale basée sur la transformée 
de Fourier 


Global Biodiversity Information Facility | Centre d’information mondial 
sur la biodiversité 


Global Human Settlement Layer | Couche mondiale des établissements humains 
Groupe d'experts intergouvernemental sur l'évolution du climat 

Grey Level Co-occurrence Matrix / Matrice de cooccurrence des niveaux de gris 
Generalized Linear Model / Modèle linéaire généralisé 


Global Precipitation Measurement | Mesure des précipitations mondiales 
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GSMaP 


HPC 
HR 
HTC 
IB 
IMERG 
INRAE 
IPF 
IPM 
IRD 
IGN 
Insee 
Iris 


JRC 


JRE 
LAV 
LMI 
LST 
MAE 
MIR 
MNE 
MNDWI 


MODIS 


MTV 
NASA 


NDBI 
NDVI 


NDWI 
OLI 
OMS 
OSM 
Peps 
PNLP 
Pnue 
PPA 
PIR 
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Global Satellite Mapping of Precipitation / Cartographie mondiale des précipitations 
par satellite 


High Performance Computing / Calcul de haute performance 
Haute résolution (spatiale) 

Hautes Terres centrales (de Madagascar) 

Indice de brillance 

Integrated Multi-satellite Retrievals for GPM 

Institut national de recherche pour l’agriculture, l'alimentation et l'environnement 
Iterative Proportional Fitting / Ajustement proportionnel itératif 
Institut Pasteur de Madagascar 

Institut de recherche pour le développement 

Institut national de l'information géographique et forestière 
Institut national de la statistique et des études économiques 
Îlots regroupés pour l'information statistique 


Joint Research Center / Centre commun de recherche (service scientifique interne 
de la Commission européenne) 


Java Runtime Environment | Environnement d’exécution Java 
Lutte antivectorielle 

Laboratoire mixte international (de IRD) 

Land Surface Temperature | Température de surface au sol 
Mean Absolute Error / Erreur moyenne absolue 

Moyen infrarouge 

Modèle numérique d’élévation 


Modified Normalized Difference Water Index / Indice d’eau par différence 
normalisée modifié 


Moderate-Resolution Imaging Spectroradiometer | Radiométre spectral pour imagerie 
de résolution moyenne 


Maladie à transmission vectorielle 


National Aeronautics and Space Administration | Administration nationale 
de l’aéronautique et de l’espace (des Etats-Unis) 


Normalized Difference Built-up Index / Indice de bâti par différence normalisée 


Normalized Difference Vegetation Index / Indice de végétation 
par différence normalisée 


Normalized Difference Water Index / Indice d’eau par différence normalisée 
Operational Land Imager 

Organisation mondiale de la santé 

OpenStreetMap 

Plateforme d’exploitation des produits Sentinel 

Programme national de lutte contre le paludisme 

Programme des Nations unies pour l'environnement 

Peste porcine africaine 


Proche infrarouge 


RC 

RVB 
SAVI 
SCO 
SEAS-OI 


SGIS 
SIG 
SMA 
SRTM 
SVM 
SWIR 
TIRS 
TIS 
THRS 
TRMM 
TVX 
UMR 
USGS 
WHO 
WMS 
WNE 
WNV 
ZSE 


Liste des sigles 


Ratio de cohérence 

Rouge, vert, bleu 

Soil Adjusted Vegetation Index / Indice de végétation ajusté au sol 
Space Climate Observatory / Observatoire spatial du climat 


Surveillance de l’environnement assisté par satellite — Océan Indien (pôle de recherche 
et de valorisation en télédétection et géomatique) 


Système de gestion d’information sanitaire 

Système d’information géographique 

Système multi-agent 

Shuttle Radar Topography Mission / Mission de topographie radar par navette spatiale 
Support Vector Machine / Machine à vecteurs de support 

Short Wave Infrared | Infrarouge à ondes courtes 

Thermal Infrared Sensor / Capteur infrarouge thermique 

Technique de l'insecte stérile 

Très haute résolution (spatiale) 

Tropical Rainfall Measuring Mission / Mission pour la mesure des pluies tropicales 
Temperature-Vegetation Index / Indice température-végétation 

Unité mixte de recherche 

United States Geological Survey / Service géologique des États-Unis 

World Health Organization | Organisation mondiale de la santé 

Web Map Service | Service de carte en ligne 

West Nile fever / Fièvre du Nil occidental 

West Nile virus / Virus du Nil occidental 


Zone de surveillance entomologique 
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Les moustiques sont vecteurs de nombreux agents pathogènes responsables 
de maladies, telles que le paludisme, la dengue, le chikungunya ou la fièvre 
jaune. Selon l'Organisation mondiale de la santé, ils provoquent plusieurs 
centaines de milliers de décès chaque année. Ils sont aussi à l'origine de 
zoonoses, comme la fièvre de la vallée du Rift et la fièvre du Nil occidental. 


Dans ce contexte, les besoins en outils opérationnels permettant d'orienter 
les actions de surveillance et de contrôle sont importants, à la fois dans les 
pays du Sud — les zones tropicales et subtropicales sont les plus touchées 
par les maladies causées par les moustiques —, mais également dans les 
pays du Nord, où l'installation de nouvelles espèces comme le moustique- 
tigre augmente le risque d'émergence de maladies. Pour répondre à ces 
besoins, les images d'observation de la Terre présentent un fort intérêt : la 
distribution dans l'espace et la dynamique temporelle des moustiques sont 
influencées par des variables climatiques (températures, précipitations, 
humidité] et environnementales (disponibilité de zones en eau, végétation), 
dont les indicateurs peuvent être dérivés d'images satellite. 


De nombreuses études récentes ont permis de développer des méthodes 
innovantes couplant télédétection et modélisation spatiale pour prédire la 
dynamique spatiale et temporelle des moustiques vecteurs et des maladies 
associées. Au-delà de l'étude de faisabilité, certaines de ces méthodes ont 
abouti à des outils et à des chaînes de traitement aujourd'hui opérationnels, 
utilisés par les acteurs de santé publique et les opérateurs chargés de la 
lutte antivectorielle. 


Cet ouvrage, destiné aux étudiants et aux chercheurs comme aux acteurs de 
santé publique, présente une synthèse de ces travaux de recherche et de ces 
outils. 
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